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This publication gives details of the 
computational subroutines available in the 
PL/I Library. These subroutines are used 
by the PL/I (F) compiler in the implementa- 
tion of PL/I built-in functions and of the 
operators used in the evaluation of PL/I 
expressions. Not all PL/I built-in func- 
tions and expression operators are support- 
ed by the PI/I Library; the compiler gener- 
ates in-line code for a small number of 
them. The details provided include timing 
figures, summaries of the mathematical 
methods used, and (where appropriate) fig- 
ures for range and accuracy. This informa- 
tion is intended to be of interest chiefly 
to those programmers concerned with the 
performance of computational subprograms. 








PREFACE 


This publication provides the PL/I pro- 
grammer with detailed information about the 
computational subroutines which are part of 
the OS/360 PL/I Library. 


The reader is assumed to be a programmer 
with a particular concern for performance 
information associated with individual 
modules. The numerical analyst is provided 
with a description of the algorithms, and a 
specification of accuracy and range, where 
these are considered to be significant. 


Useful background reading is provided in 
the following IBM publications: 


IBM System/360 Principles of Operation, 
Form A22-6821 


IBM  System/360 Operating System PL/I: 
Language Specifications, Form C28-6571 





IBM System/360 Operating System: PL/I (F) 
Programmer's Guide, Form C28-6594 


IBM System/360 Operating System: Assem- 
pler Language, Form C28-6514 


Copies of this and other IBM publications can be obtained through IBM 


Branch Offices. 


A form for reader's comments appears at the back of this publication. 
It may be mailed directly to IBM, Address any additional comments 
concerning this publication to the IBM Corporation, Department D39,1271 


Avenue of the Americas, New York, N.Y. 10020. 


© International Business Machines Corporation, 1966 
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INTRODUCTION 


The PL/I Library computational subrou- 
tines provide support for the operators and 
functions of the PL/I language in four 
major categories: 

1. Bit and Character Strings 
2. Arithmetic 

3. Mathematical 

4. Arrays 

These subroutines have been designed to 
allow their use in a multi-tasking environ- 
ment. 


This publication gives detailed informa- 
tion in each of the four sections mentioned 


COMPUTATIONAL SUBROUTINES 


above with respect to performance, accura- 
cy, choice of algorithm, and range of 
values handled (where appropriate). 


A number of exceptional conditions may 
arise in the execution of the library 
subroutines. Many of these are not direct- 
ly related to PL/I ON conditions. The 
method of treatment in these cases is to 
write a diagnostic message and raise the 
ERROR condition. This allows the user the 
opportunity of investigating the error by 
use of the ONCODE built-in function in his 
ON ERROR unit and of making a choice on the 
action which he wishes to take. Full 
details of the diagnostic messages printed 
at object time and of the ONCODE values 
associated with them are specified in IBM 


System/360 Operating System: PL/I (F) 
Programmer's Guide, Form C28-6594. 


CHAPTER 1: BIT AND CHARACTER STRING OPERATIONS AND FUNCTIONS 


The Library string package contains 
modules for handling bit and character 
strings. Generally, a string function or 


operator is supported by only one module, 
but in the interests of efficiency some of 
the kit string operators are provided with 
additional modules to deal with byte- 
aligned input data. 


Information, 


Form A22-6825. They are 
intended to be simple enough to give a good 
general guide to performance while averag- 
ing out many logical variations, the effect 
of which is comparatively small. 


A complete list of the modules provided 
in the Library string package is given in 


Execution times are based on information Figure 1. 
in IBM System/360 Instruction Timing 
v eU LA | Si ee e ee ee dr 
| PL/I | PL/I | Bit String [Character | 
| Operation | Function [---------- T------------ ] String | 
| | General  |Byte-aligned| | 
p----------------- }---------- ł----------}------------}------ =- 1 
| "And" (8) | - | Use BOOL | IHEBSA | - | 
| 'or' (|) | - | Use BOOL | IHEBSO | - | 
| "Not" (4) | - | Use BOOL | IHEBSN | - | 
| Concatenate (||)| REPEAT |  IHEBSK | - | IHECSK | 
| Compare | x | IHEBSD | IHEBSC | IHECSC | 
| Assign | - |Use IHEBSK | IHEBSM | IHECSM | 
| Fill | - | IHEBSM | - | IHECSM | 
| - | HIGH/LOW | - | - | IHECSM | 
| - | SUBSTR |  IHEBSS | - | IHECSS | 
| - | INDEX | IHEBSI | - | IHECSI | 
| - | BOOL | IHEBSF | - | - | 
A iaskie a AAA DEAR EEE APPS Le aa J 


Figure 1. 


BIT STRING OPERATIONS 


The 'And' Operator (6) 


Module Name: IHEBSA 


Entry Point: IHEBSAO 


Function: 


To implement the 'and' operator between 
two byte-aligned bit strings, placing the 
result in a byte-aligned target field. 


Bit and Character String Operations and Functions 


Method: 


The current length of the 
is set equal either to the maximum of 
those of the operands, or to the maximum 
length of the target field (when trunca- 
tion is necessary to avoid exceeding the 
length of this field). The strings are 
'and'ed together for a length equal to 
the minimum of the lengths of the oper- 
ands, and the result is extended with 
zeros, if necessary, up to the current 
length calculated for the target field. 


target string 


Implementation: 


M 


E 


F 


e Module size: 296 bytes 


e Execution time: 


Let L4 - the length of the shorter 

string 

Le = the difference between the 
string lengths 

Mi = FLOOR((Ly - 1)/2048) 

Ni = MOD(L;, 2048) - 1 

Pi = FLOOR(N;/8) 

S = SIGN(L2) 

i =1,2 


Then the approximate execution times in 
microseconds for the System/360 models 
given below are obtained from the fol- 
lowing formula: 


a + b*M4 + c¥P, + S*(d + e*Ma + f*P2) 


eeepc era pene poto TS 1 
| 30 | 40 | 50 | 65 | 75 | 
~--}------ ł------- ł------- ł------- ł------ 1 
a | 2184 | 693 | 296 | 86 | 57 | 
b | 2569 |1543 | 757 | 243 1202 | 
c | 9 | 5.6 | 2.8 | 0.9 | 0.7 | 
d | 363 | 117 | 48 | 15 | 9 | 
e | 1171 | 696 | 316 | 107 | 95 | 
£ | 4 | 2.5 | I. | 0.4 | 0.3 | 
| | j | | | 
a'| 1437 | 466 | 200 | 60 | 43 | 
— — ee dae RENEE dica sii J 
ote: If the two strings are equal in 
length and end on a byte boundary, 


then a is replaced by a‘. 


The 'Or' Operator (|) 


odule Name: IHEBSO 
ntry Point: IHEBSOO 
unction: 


To implement the ‘or’ operation between 
two byte-aligned bit strings, placing the 
result in a byte-aligned target field. 


Chapter 1. 


Method: 


The current length of the target string 
is set equal to either the maximum of 
those of the operands or to the maximum 
length of the target field (when trunca- 
tion is necessary to avoid exceeding the 
length of this field). The strings are 
'or'ed together for a length equal to the 
minimum of the lengths of the operands 
and the remainder of the longer string is 
moved into the target field up to the 


current length calculated for it; the 
remainder of the target field is left 
unchanged. 
Implementation: 
e Module size: 312 bytes 
e Execution time: 
Let L4, = the length of the shorter 
string 
L2 = the difference between the 


Note: 


string lengths 

Mj; = FLOOR((Lj - 1)/2048) 
MOD(L4,2048) - 1 
Pj = FLOOR(N,/8) 
SIGN(La) 
1,2 

Then the approximate execution times in 
microseconds for the System/360 models 
given below are obtained from the fol- 
lowing formula: 

a + b*M, + c*P, 


+ S*(d + e*M, + f*P3) 


T queam quee par PESAS 1 

| 30 | 40 | 50 | 65 | 75 | 
-------- ł------- ł------- ł------- ł------ 4 
a | 2294 | 724 | 316 | 97 | 65 | 
b | 2569 |1543 | 757 | 243 | 202 | 
c | 9 | 5.6 | 2.8 | 0.9 | 0.7 | 
d | 364 | 117 | 53 | 18 | 15 | 
e | 1146 | 688 | 311 | 103 | 92 | 
£ | 4 | 2.5 | 1.1 | 0.4 | 0.3 | 
| | | | | | 
a'| 1437 | 466 | 200 | 60 | 43 | 
E ccc dl ccc AE dust Loros J 
If the two strings are equal in 
length and end on a byte boundary, 


then a is replaced by a'. 
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The 'Not' Operator (1) 


Module Name: IHEBSN 

Entry Point: IHEBSNO 

Function: 
To implement the "not' operator for a 
byte-aligned bit string, placing the 


result in a byte-aligned target field. 


Method: 


The current length of the target string 
is set equal to either the current length 
of the operand or to the maximum length 
of the target field (when truncation is 
necessary to avoid exceeding the length 
of this field). The target field is set 
to a string of 1's for a length equal to 
its calculated current’ length and the 
result is obtained by an ‘exclusive or' 
with the operand. The remainder of the 
target field beyond the calculated cur- 
rent length is left unchanged. 


Implementation: 
e Module size: 192 bytes 


e Execution time: 


Let L4, = the string length 
M4 = FLOOR((L4 - 1)/2048) 
N, = MOD(L,,2048) - 1 
Pa = FLOOR(N,/8) 


Then the approximate execution times in 
microseconds for the System/360 models 
given below are obtained from the fol- 
lowing formula: 

a + b*M, + C*P, 


a’ + b*M, + c*P, 


| iero Neri E A e casa e Tomte A 1 
| | 30 | 40 | 50 | 65 | 75 | 
}---4------ j------- +------- ł------- ł------ { 
| a | 1599 | 530 | 245 | 72 | 52 | 
| b | 2544 |1533 | 754 | 242  |201 | 
[co | 9| 5.6| 2.8 | 0.9 | 0.7 | 
| | | | | | | 
| at] 1183 | 429 | 192 | 55 | 40 | 
SS Ses dom do Al----—— dh J 
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Bit String Concatenate/REPEAT 


Module Name: IHEBSK 


Entry Points: 


Entry 

Operation point 
Concatenate (||) IHEBSKK 
REPEAT(Bit string,n) IHEBSKR 


Function: 


IHEBSKK: to concatenate two bit strings 
into a target field. 


IHEBSKR: to concatenate n + 1 instances 
of the single source string into a 
target field. If n < 0, the result is 


the string itself. 


Method: 


The current length of the target field is 
made equal to the smaller of two values: 
the sum of the current lengths of the 
source strings, and the maximum length of 
the target field. Both entry points use 
a loop which obtains data from the source 
fields, aligns it correctly, and moves it 
to the target field. The remainder of 
the target field beyond the calculated 
current length is left unchanged. 


Implementation: 


Y 


e Module size: 328 bytes 


e Execution time: 


Let L4,,La2 = the lengths of the two 
strings 
F į = FLOOR (Lj/32) 
R =n + 1, ifn20 
= 1, if n <0 
i = 1,2 


Then the approximate execution times in 
microseconds for the System/360 models 
given below are obtained from the for- 
mula for the appropriate entry point: 


IHEBSKK: a + b* (F4 + Fo) 
IHEBSKR: c + R*(d + b*F,) 


FYT—"T UE ona UN A A E Moet sg TI Passo 1 
| | 30 | 40 | 50 | 65 | 75 | 
}---+------ ł------- ł------- ł------- ł------ 
| a | 3497 | 1311 | 445 | 126 | 74.7] 
| b | 345 | 111 | 42 | 11.8 | 8.1 | 
| c | 594 | 187 | 73 | 16.4 | 11.0 | 
| à | 1505 | 479 | 196 | 56.3 | 33.2| 
Loriano dr do d---—---4------—- J 


Bit String Byte-aligned Comparison 


Module Name: IHEBSC 


Entry Point: IHEBSCO 


Function: 


To compare two byte-aligned bit strings 
and to return a condition code as bits 2 
and 3 of a full-word target field as 
follows: 


00 if strings are equal 

01 if first string compares low at the 
first inequality 

10 if first string compares high at the 
first inequality 


The shorter string is treated as though 
extended with zeros to the length of the 
longer. 


The first kyte of the target field is 
also used to preserve the program mask in 
the PSW for the calling routine. This 
byte contains: 


Bits Contents 
0 to 1 Instruction length code 01 


2 to 3 Condition code as above 
4 to 7 Program mask (calling routine) 


Method: 


The two strings are compared up to the 
current length of the shorter string. 
The remainder of the longer string is 
compared with zeros. 


Implementation: 
e Module size: 272 bytes 
e Execution time: 


Let L, = the number of bytes compared 
up to the first inequality 


La = the number of bytes in the 
additional part of the longer 
string, compared to zeros if 
necessary 

Ma = FLOOR ((L, ~- 1)/256) 

Na = MOD (Li ~ 1,256) 

S = SIGN (La) 

Then the appropriate execution times in 
microseconds for the System/360 models 
given below are obtained from the fol- 
lowing formula: 


a + b*M4 * C*N4 * S*(d * e*La) 


po Wu E TOT TST eee ADT E NET) Wi ig gt = oe 1 
| | 30 | 40 | 50 | 65 | 75 | 
}---+------ j------- +------- j------- ł------ { 
| a | 981 | 304 | 128 | 40.3 | 27.0 | 
| b | 1477 | 794 | 289 | 114 {110 | 
[ c | 5| 2.8] 1.0] 0.4 | 0.8 | 
| di 846 | 256 | 107 | 29.2 | 18.8 | 
le] 99} 40 | 16 | 4.9 | 3.6 | 
Lllciclcloclocloo-2c- Al-2c-oc- Lua A i ss d 


Bit String General Comparison 


Module Name: IHEBSD 


Entry Point: IHEBSDO 


Function: 


To compare two bit strings and return a 
condition code as bits 2 and 3 of a 
full-word target field as follows: 


00 if strings are equal 

01 if first string compares low at the 
first inequality 

10 if first string compares high at the 
first inequality 


The shorter string is treated as though 
extended with zeros to the length of the 
longer. 
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The first byte of the target field is 
also used to preserve the program mask in 
the PSW for the calling routine. This 
byte contains: 


Bits Contents 
0 to 1 Instruction length code 01 
2 to 3 Condition code as above 
4 to 7 Program mask (calling routine) 
Method: 
The two strings are compared up to the 
current length of the shorter string. 
The remainder of the longer string is 
compared with zeros. 
Implementation: 
e Module size: 192 bytes 
e Execution time: 
the number of 32-bit words 


compared up to the first in- 
equality 


Let Fi 


F2 = the number of 32-bit words 
compared to zero if necessary 


Then the approximate execution times in 
microseconds for the System/360 models 
given below are obtained from the fol- 
lowing formula: 


a + b*F4, + c*F> 


a each ae) gcc quem Um quem m prm 1 
| | 30 | 40 | 50 | 65 | 75 | 
E 4------- ł------- ł------- ł------ 1 
|a | 1113 | 344 | 120 | 31.7 | 18.2 
| b | 624 | 194 | 76 | 21.3 | 14.9 | 
| c | 437 | 139 | 55 | 15.3 | 11.4 | 
base LS OP ica Ls Loco. > Lota J 
Bit String Assign/Fill 
Module Name: IHEBSM 
Entry Points: 
Entry 
Operation point 
Fixed-length assign IHEBSMF 
Variakle-length assign IHEBSMV 
Zero fill only IHEBSMZ 
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Function: 


IHEBSMF: to assign a byte-aligned string 
to a byte-aligned fixed-length target, 
filling out with zero bits if neces- 
sary. 


IHEBSMV: to assign a byte-aligned string 
to a byte-aligned variable-length tar- 
get. 


IHEBSMZ: to fill out the target area from 
its current length to its maximum 
length with zero bits. 


Method: 


IHEBSMF: the minimum of the source  cur- 
rent length and the target maximum 
length is calculated and the source 
string is moved to the target for a 
length equal to this length. Zero 
filling of the target is performed if 
necessary. The current length of the 
target is set equal to the maximum 
length. 


IHEBSMV: the source string is moved to 
the target field as above, but without 
zero filling. The current length of 
the target is set appropriately. 


IHEBSMZ: zeros are propagated in the 
target from the current length to the 
maximum length. The current length of 
the target is set equal to the maximum 
length. 


implementation: 


e Module size: 384 bytes 


e Execution time: 


Let B, the bit length of the field to 


be assigned 


Bz = the bit length for zero fill- 
ing 


Li = FLOOR(P4/8) 


= 
pH- 
Il 


FLOOR (Lj/256) 


a 
H- 
Il 


MOD (L4, 256) 
S = SIGN (L2) 


i = 1,2 


O 


Then the approximate execution times in 
microseconds for the System/360 models 
given 
mula for the appropriate entry point: 


IHEBSMF: f + Lb*M, + CEN, + S*(g + e*M, 
+ C*N2) 
IHEBSMV: a + b*M, + C*N4 
IHEBSMZ: d + e*M, + c*Na 
DER: ca foe ps Ao a PS > ce qmm 1 
| 30 | 40 | 50 | 65 | 75 | 
np + 4------- fono 1 
a | 1379 | 442 | 192 | 58.9 | 41.2 | 
k | 1196 | 706 | 319 | 108 | 95.8 | 
c | 4 | 2.5 | 1.1 | 0.4 | 0.3 | 
a | 1865 | 589 | 243 | 73.6 | 48.0 | 
e | 1171 | 696 | 316 | 107 | 92.9 | 
£ | 1897 | 623 | 266 | 81.1 | 56.6 | 
g | 1086 | 349 | 150 | 46.7 | 31.8 | 
| | | | | | 
a'| 1079 | 343 | 140 | 42.9 | 29.1 | 
f'| 1597 | 524 | 214 | 65.1 | 44.5 | 
Fu S rata Latas a Ls SER J 
ote: a‘ and f' replace a and f respective- 
ly if MOD (B,,8) = 0 
ther Information: 


This routine supplies assignment of byte- 
aligned bit strings of roth fixed and 
variable lengths. Non-aligned strings 
may be assigned by using the REPEAT entry 
IHEBSKR with n equal to O. Any filling 
required for fixed length strings can 


then be obtained using the IHEBSMZ entry 
described akove. 
BIT STRING FUNCTIONS 
Bit String SUBSTR 
Module Name: IHEBSS 
Entry Points: 
Entry 
Oreration point 
SUBSTR(Bit-string,i) IHEBSS2 
SUBSTR(Bit-string,i, j) IHEBSS3 
Function: 
TO produce a string dope vector 


describing the SUBSTR pseudo-variakle and 
function of a bit-string. 


Chapter 1. 


below are obtained from the for- . 


Module Name: 


Entry Point: 


Method: 
Arithmetic is performed according to the 
function definition, using the current 
length of the argument string. The 


result describes a fixed-length string. 


Implementation: 


e Module size: 192 


bytes 
e Execution times: 


Approximate execution time in micro- 
seconds for the System/360 models given 
below is shown for each entry point: 


+] 


PL + — + —— 


Bit String INDEX 


IHEBSI 


IHEBSIO 


Function: 


TO compare two bit strings to see if the 
second is identical to a substring of the 
first, and, if it is, to produce a binary 
integer (the index) which indicates the 
first bit position in the first string at 
which such a substring begins. If no 
such index is found, or if either string 
is null, the function value is zero. 


Method: 


and com- 
of the two strings in 


The index is found by shifting 
paring portions 
registers. 
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implementation: 


e Module size: 296 bytes 


e Execution time: 
Let Fj = the length (in words) actually 
processed at the ith  compari- 
son (the length up to the 
first inequality or the length 
of the second string if no 
inequality is found) 


(length of first string) - 
(length of second string) 


I = the resulting index, except 
when this is 0 with B 2 0, in 
which case I = B + 1 
J = MOD(I,32) 
i = 1,n 
Then the approximate execution times in 
microseconds for the System/360 models 


given below are obtained from the fol- 
lowing formula: 


I 
a + b*I + c*J + d* 5, (Fi - 1) 


i-i 

——T--—----— Tea Goo es po ————T T-———-— 1 

| 30 | 40 | 50 | 65 | 75 | 
e 1------- ł------- 1------- f------ 1 
a | 2026 | 646 | 246 | 63.6 | 37.4 | 
b | 422 | 124 | 48 | 13.9 | 9.2 | 
c | 104 | 41 | 18 | 8.3 | 3.1 | 
d | 543 | 180 | 73 | 18.9 | 13.2 | 
A A E E A AS deem J 


Bit String BOOL (Boolean Function) 


Module Name: 


Entry Point: 


IHEBSF 


IHEBSFO 


Function: 


To take two source strings and perform 
one of the sixteen possible logical oper- 
ations between corresponding bits. The 
particular operation performed is defined 
by inserting the bit pattern - n,4ngnan,- 
yielded by the third argument into the 
table below: 
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Method: 


The current length of the target string 
is set equal to either the maximum of the 
current lengths of the source strings or 
to the maximum length of the target field 
(when truncation is necessary to avoid 
exceeding the length of this field). The 
necessary operation is performed on the 
strings and the result stored in the 
target field. If one string is shorter 
than the other, it is regarded as being 
extended on the right with zeros up to 
the length of the longer. The field 
between the calculated current length and 
the maximum length of the target is left 
unchanged. 


Implementation: 


e Module size: 480 bytes 


e Execution time: 


Let B4, = the bit length of the shorter 
string 
Ba = the difference in bit lengths 


of the strings 
Fj = CEIL(B,/32) 
i = 1,2 
Then the approximate execution times in 
microseconds for the System/360 models 
given below are obtained from the fol- 
lowing formula: 


a + b*F, + C*Fa 


FOTO Te IS nicae d Met EO n dedans 1 
| 30 | 40 | 50 | 65 | 75 | 
E---t------ $------- ł------- ł------- ł------ 1 
a | 2275 | 767 | 291 | 75.0 | 46.5 | 
b| 728 | 230 | 95 | 30.0 | 19.4 | 
c | 501 | 162 | 68 | 21.8 | 15.2 | 
b --cee dot Larissa dou comics areas Lise >= J 


CHARACTER STRING OPERATIONS 


Character String Concatenate/REPEAT 


Module Name: IHECSK 
Entry Points: 
Entry 
Operation point 
Concatenate (| |) IHECSKK 
REPEAT IHECSKR 


(Character string,n) 


Function: 


IHECSKK: to concatenate two character 
strings into a target field. 

IHECSKR: to concatenate n + 1 instances 
of the single source string into a 


target field. If n < 0, the result is 
the string itself. 


Method: 


The current length of the target field is 
made equal to the smaller of two values: 
the sum of the current lengths of the 
source fields, and the maximum length of 
the target field. The source strings are 
then moved to the target. Characters 
beyond the range of the target current 
length remain unaltered. 


Implementation: 


e Module size: 208 bytes 


e Execution time: 


Let L4,,L2 = the lengths of the source 
strings 
M 1 = FLOOR((L; - 12/256) 
Nj = MOD(Lj - 1,256) 
R =n+1, if ne2 0 
= 1, ifn< 0 
i 2 1,2 


Then the approximate execution times in 
microseconds for the System/360 models 
given below are obtained from the fol- 
lowing formulas: 

IHECSKK: a + k*(M, + Ma) + c*(N, + N) 
IHECSKR: R*(b*M, + c*N, + d) + e 


T Tose pessoa pte quu 1 
| 30 | 40 | 50 | 65 | 75 | 
+------ }------- ł------- ł------- T 1 
a | 1407 | 491 | 209 | 63.7 | 46.4 | 
b | 1196 | 706 | 319 | 108 | 95.8 | 
c | 4 | 2.5 | 1.1 | 0.4 | 043 | 
a | 703 | 236 | 101 { 31.7 | 22.0 | 
e| 108 | 54.6 | 26.0 | 5.4 | 6.2 | 
isse Leuca Li de Los 4 
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Module Name: 


Entry Point: 


Character String Compare 


IHECSC 


IHECSCO 


Function: 


To compare two character strings and to 
return a condition code as bits 2 and 3 
of a full-word target field as follows: 


00 if strings are equal 

01 if first string compares low at the 
first inequality 

10 if the first string compares high at 
the first inequality 


The shorter string is treated as though 
extended with blanks to the length of the 
longer one. 


The first byte of the target field is 
also used to preserve the program mask in 
the PSW for the calling routine. This 
byte contains: 

Bits Contents 
0 to 1 Instruction length code 01 
2 to 3 Condition code as above 
4 to 7 Program mask (calling routine) 


Method: 


The two strings are compared in storage. 
If the strings are of different lengths 
and are identical up to the length of the 
shorter, the remainder of the longer is 
compared with blanks. 


Implementation: 


e Module size: 200 bytes 


e Execution time: 


Let L4 the length of the strings com- 
pared up to the first in- 
equality (proceeding left to 


right) 


La the length of the additional 
part of the longer string com- 


pared with blanks if necessary 


FLOOR((Lji- 1)/256) 


= 
H- 
Il 


Ni = MOD(Lj - 1,256) 
S = SIGN(La) 


i =1,2 


Strings: Character String Operations 15 


Then the approximate execution times in 
microseconds for the System/360 models 
given below are obtained from the fol- 
lowing formula: 


a t b*M4 + C*N4 + S*(d + e*M3 + CN) 


CASTA MAR e ir e | Ve oe ER T ^ 77 2 
| [| 30 | 40 | 50 | 65 | 75 | 
p---4------ ł------- ł------- ł------- i------ 1 
| a | 849 | 284 | 116 | 36.8 | 27.1 | 
| b | 1469 | 790 | 289 | 115 |110 | 
| c | S.t 298 1.0| 0.4 | 0.4 |] 
| aj 620 | 210 | 88 | 27.3 | 21.7 | 
| e | 1474 | 794 | 290 | 114 |110 | 
L2-llIllzllzl Loo Se do ose es E ee J 


Character String Assign/Fill/HIGH/LOW 


Module Name: IHECSM 


Entry Points: 


Entry 

Operation point 
Fixed-length assign IHECSMF 
Variable-length assign IHECSMV 
Blank fill only IHECSMB 
HIGH IHECSMH 
LOW IHECSML 

Function: 


IHECSMF: to assign a character string to 
a fixed-length target, filling out with 
blanks if necessary. 


IHECSMV: to assign a character string to 
a variakle-length target. 


IHECSMB: to fill out the target field 
from its current length to its maximum 
length with blanks. 


IHECSMH: to fill a target field with the 
highest character in the collating 
sequence, up to its current length. 


IHECSML: to fill the target field with 
the lowest character in the collating 
sequence, ur to its current length. 


Method: 

IHECSMF: The minimum of the source cur- 
rent length and the target maximum 
length is calculated and the source 
string is moved to the target for a 
length equal to this length. Filling 


of the target with blanks 
target maximum 


up to the 
length is performed if 


necessary. The current length of the 
target is set equal to its maximum 
length. 
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IHECSMV: moves the string as above, but 


without blank filling. The current 
length of the target is set  appropri- 
ately. 

IHECSMB:  propagates blanks and sets the 


current length of the target 
its maximum length. 


equal to 


IHECSMH,  IHECSML: uses part of the blank 
fill routine to propagate the highest 
or lowest character in the collating 
Sequence up to the current length of 
the target. 


Implementation: 
e Module size: 280 bytes 
* Execution time: 


Let Lı = either the specified length 


(for IHECSMH/L/B) or the 
length for blank filling 
(IHECSMF) 


L2 = the length of the 
the source 
(IHECSMF/V) 


shorter of 
and target fields 
M, = FLOOR((L, - 2)/256) 


Na = MOD(L, - 2,256) (where La 2 2 
for both M, and N.) 


Ma = FLOOR((L2 - 1)/256) 
Na = MOD (Là = 1, 256) 


S = SIGN(L,) 


Then the approximate execution times in 
microseconds for the System /360 models 
given below are obtained from the fol- 
lowing formulas: 


IHECSMF: h + g*M, + c*N, + S*(i + b*M, 


+ C*N4) 

IHECSMV: f + g*M, + c*N2 

IHECSMB: e + b*M, + c*N, 

IHECSMH: a + b*M, + CXN4 

IHECSML: d + b*M, + C*N, 
p[^77T*4-25-- gres ToL TSK TH Tessere y.===== 1 
| | 30 | 40 | 50 | 65 | 75 | 
}---+------ ł------- ł------- d---- I------ 4 
| a | 830 | 282 | 120 | 37.0 | 26.8 | 
| b | 1171 | 696 | 316 | 107 | 95.2 | 
| c | 4 | 245: || 1.1 | 0.4 | 0.3 | 
pa | 799 | 268 | 116 | 35.5 | 25.6 | 
| e | 892 | 300 | 130 | 40.5 | 28.0 | 
| £ | 790 | 267 | 114 | 34.6 | 24.9 | 
| g | 1196 | 706 | 319 | 108 | 95.8 | 
| h | 1298 | 446 | 190 | 57.4 | 40.6 | 
| i | 400 | 129 | 57.6 | 18.5 | 13.7 | 
lo AO eee ee AA AA Loans. 4 
CHARACTER STRING FUNCTIONS 
Character String SUBSTR 
Module Name: IHECSS 
Entry Points: 

Entry 
Operation point 
SUBSTR(Character-string,i) IHECSS2 
SUBSTR(Character-string,i, Jj) IHECSS3 


Function: 


To produce a string dope vector describ- 
ing the SUBSTR pseudo-variable and func- 
tion of a character string. 


Method: 
Arithmetic is performed according to the 
function definition, using the current 
length of the argument string. The 


result descrikes a fixed-length string. 


Chapter 1. 


Implementation: 
e Module size: 176 bytes 


e Execution times: 


Approximate execution times in micro- 
seconds for the System/360 models given 
below are as shown for each entry 
point: 
CECT FEE SS Mr es ee ria 
|Entry | 30 | 40 | 50 | 65 | 75 | 
| | Point | | | | | | 
---~---}------ t------}------}-----}------4 
|IHECSS2| 887 | 310 | 127 | 36.9| 26.2 | 
------ from apa fr nena fen fr] 
|IHECSS3| 1018 | 347 | 143 | 41.0] 29.1 | 
USS dellla a eee Sal discs J 


Character String INDEX 


Module Name: IHECSI 


Entry Point: IHECSIO 


Function: 


TO compare two character strings to see 
if the second is identical to a substring 
of the first, and, if it is,to produce a 
binary integer (the index) which indic- 
ates the first character position in the 


first string at which such a substring 
begins. If no such index is found, or if 
either string is null, the function value 
is zero. 

Method: 


The point required is located by compar- 
ing in storage the second string with a 
corresponding number of characters in the 
first string. 
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Implementation: 
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e Module size: 168 bytes 


e Execution time: 


Let Li. 


Mi 


the length processed at the 
ith comparison (the length up 
to the first inequality, or 
the length of the second 
string if no inequality is 
found) 


FLOOR (14/256) 


MOD(L4, 256) 


(length of first string) - 
(length of second string) 


the resulting index, except 
when this is 0 with B 20, in 


which case I = B + 1 


1,n 


Then the approximate execution times in 
microseconds for the System/360 models 
given below is obtained from the fol- 
lowing formula: 


I 
a + b*I + È) (c*Mj + d*N4) 
i=1 
n ee Win gee we FS SEES ære IDA DT AUR Ber e 1 
| 30 | 40 | 50 | 65 | 75 | 
pecore peces > EE teres prose 1 
| 954 | 314 | 127 | 38.2 | 24.5 | 
| 160 | 59.9 | 28.8 | 10.1 | 9.3 | 
| 1427 | 776 | 280 | 110 [107 | 
| 5| 2.8]| 1.0 | 0.4 | 0.4 | 
docce docs docs lc an an cn n ee den ci ee a ee d 


CHAPTER 2: ARITHMETIC OPERATIONS AND FUNCTIONS 


Library arithmetic modules support all 
those arithmetic generic functions and 
operators for which the compilers neither 
produce in-line code nor (as for the func- 
tions FIXED, FLOAT, BINARY and DECIMAL) use 
parts of the conversion package. 


Statistics for accuracy of floating- 
point modules are given where considered 
meaningful and helpful; an explanation of 
their use is given in the chapter on 
mathematical routines. Precise results are 
obtained from all fixed-point modules 
except complex division and complex ABS, 
where small truncation errors inevitably 
occur, and the ADD function (fixed 
decimal), in which the effect of truncation 
errors depends on the relative values of 
the scale factors of the arguments. 


Any restrictions on the admissibility of 
arguments are noted under the headings 
'Range' and 'Error and Exceptional 
Conditions'. 


Range: This states any ranges of arguments 
which a module assumes to have been exclud- 
ed prior to its being called. 


Error and Exceptional Conditions: These 
cover conditions which may result from the 


use of a routine; they are listed in four 
categories: 


P -- Programmed conditions in the module 
concerned. Programmed tests are 
made where this is not too costly 
and, if an invalid argument is 
found, a branch is taken to the 
entry point IHEERRC of the execution 


error package (EXEP). This results 
in the printing of an appropriate 
message and in the ERROR condition 
being raised. 


I -- Interrupt conditions in the module 
concerned. For those routines where 
SIZE and FIXEDOVERFLOW are detected 
by programmed tests or where hard- 
ware interruptions. may occur, the 
OVERFLOW,  UNDERFLOW, FIXEDOVERFLOW, 
SIZE and ZERODIVIDE conditions pass 
to the ON handler (IHEERR) and are 
treated in the normal way. The 
machine is assumed to be enabled for 
all interruptions except signifi- 
cance, which is masked off. 


O -- Programmed conditions in modules 
called by the module concerned. 
These occur when invalid arguments 


are detected in the module called. 


H -- As ES but the interrupt conditions 
occur in the modules called by the 
module concerned. 


Speed 


The average execution times given are 
based on the IBM System/360 Instruction 
Timing Information, Form A22-6825. These 
times include the times taken by the 
modules called. 


A summary of the Library arithmetic 
modules is given in Figures 2 and 3. 
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| ARITHMETIC OPERATION 
fon—------------------___n_no_mnoeco-o esca 
| Operation | Binary 
| | fixed 
pa re one È EA 
| Real Operations 
E a presse 
| Integer exponentiation: x**n | IHEXIB 
| General exponentiation: x**y | = 

| Shift-and-assign, Shift-and-load | = 

p-----7-].-]-.-]--------------2---------------- asias 


T 

| Multiplication/division: Za4*z2,24/22 | IHEMZU 
| Multiplication: z4*za | - 

| Division: 24/22 | = 

| Integer exponentiation: z**n | IHEXIU 
| General exponentiation: 24**Z2 | = 

ae ae III IE LIA demum eram 
Figure 2.  Arithmetic Operations 

POSTE m e n uem 1 
| ARITHMETIC FUNCTIONS | 
[pesmet TaT e D ee TTI 1 
| Function | Binary | Decimal| Short | Long | 
| | fixed | fixed | float |. float | 
}---------- br Las AA gll y 
l Real Arguments | 
}---------- T-------- q--------1-------- -------- 1 
| MAX, MIN | IHEMXB | IHEMXD | IHEMXS | IHEMXL | 
| ADD | - | IHEADD | - | - | 
}---------- Pose bosses AAA LA 4 
| Complex Arguments | 
pS poe que e RAR pS? > 7 
| ADD | - | IHEADV | - | =- | 
| MULTIPLY | IHEMPU | IHEMPV | - | - | 
| DIVIDE | IHEDVU | IHEDVV |  - | - | 
| ABS | IHEABU | IHEABV | IHEABW | IHEABZ | 
AA Le AE Lts A RE ER J 
Figure 3. Arithmetic Functions 

REAL OPERATIONS Me 


Positive Integer Exponentiation (fixed 
binary) 


Module Name: IHEXIB 


Entry Point: IHEXIBO 


Function: 


TO calculate x**n, where n is a positive 
integer. 


20 


S 


fixed | 


I 
1 
1 
I 
l 
i 
1 
I 
H 


H———4 


p————. 


thod: 


The result is set initially to the 
The final result is 


of the ar 


IHEXXS 


IHEMZW 
IHEDZW 
IHEXIW 
IHEXXW 


gument. 


Pp_—d4 FRA 


-—————4 


IHEXXL 


IHEMZZ 
IHEDZZ 
IHEXIZ 
IHEXXZ 


then obtained by repeated 


the argument. 


Range: 


0 < n < 2** 
The  precis 
precision 


implies 
such cases. 


31 
ion 


(p,q), 


rules 


this 
called only if n*(p + 1) 
that n < 32/(p + 1) 


© L 


squaring of 
this value or squaring and multiplying by 


of PL/I 
further restriction in that if 
module will be 


^ IA 


16 for all 


Implementation: 
e Module size: 88 bytes 


e Execution time: 


Let M = number of significant bits in 
the exponent 
N = number of 1 bits in the expo- 
nent 


Then the approximate execution times in 
microseconds for the System /360 models 


given beiow are obtained from the  fol- 
lowing formula: 
-a * b*M * c*N 

fob SEES ToS es quee TOS cre rane 1 
| | 30 | 40 | 50 | 65 | 75 | 
}---+------ ł------- ł------- i------ ł------ 1 
| a 238 | 9 | 5.5 | -1.46| -0.6 | 
| b | 708 | 188 | 63.8 | 15.6 | 9.8 | 
| c | 335 | 94 | 33.0 | 6.1 | 3.9 | 
A AD decesos SA AA dirne ADELS J 


Positive Integer Exponentiation (fixed 
decimal) 


Module Name: IHEXID 


Entry Point: IHEXIDO 


Function: 


To calculate x**n, where n is a positive 
integer. 


Method: 


The result is set initially to the value 
of the argument. The final result is 
then obtained by repeated squaring of 
this value or squaring and multiplying by 
the argument. 


Range: 
The precision rules of PL/I impose the 


restriction that if x has a precision 
(p,q), this module will be called only if 


n*(p + 1) - 1 < 15. This implies that 
n < 16/(p + 1) 8 for all such cases 
and, in fact, this module will operate 


only for the range 0 <n s 8. 


Implementation: 
e Module size: 136 bytes 
e Execution time: 


Let M = number of significant bits in 
the exponent 
N = number 
nent 


of 1 bits in the expo- 


Then the approximate execution times in 
microseconds for the System/360 models 
given below are obtained from the fol- 
lowing formula: 


-a + b*M + c*N 


E ie T Toce a ee Ts 1 
| | 30 | 40 | 50 | 65 | 75 | 
E---1------ persone ł------- fo i------ 1 
{al 580 | 48 | -16.3 | 37.1 | 44.5 | 
| b | 1113 | 279 | 91 | 45.5 | 40.3 | 
| c | 752 | 192 | 65 | 38.8 | 37.2 | 
RN EAS K cae ames eas discos dlllccaudu am J 


Integer Exponentiation (floating-point) 


Module Names and Entry Points: 


Module Entry 

Argument name point 

Short float  IHEXIS IHEXISO 

Long float IHEXIL IHEXILO 
Function: 


To calculate x**n, where n is an integer 
between -2**31 and 2**31 - 1 inclusive. 


Method: 


If the exponent is zero and the argument 
non-zero, the result 1 is returned 
immediately. Otherwise the result is set 
initially to the value of the argument 
and the exponent is made positive. The 
argument is raised to this positive power 
by repeated squaring of the contents of 
the result field or squaring and multi- 
plying by the argument. Then, if the 
exponent was negative, the reciprocal of 
the result is taken, otherwise it is left 
unchanged. 
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Accuracy: Then the approximate execution times in 


microseconds for the System/360 models 
given below are obtained from the fol- 
lowing formulas: 


The values given here are for the rela- 
tive error divided by the exponent for 
exponents between 2 and 1023; the argu- 
ments are uniformly distributed over the 


a + b*M + c*N for positive exponents 
full range for each exponent for which 


neither OVERFLOW nor UNDERFLOW occurs. a'* b¥M + c*N for negative exponents 
There are 2**(10 - k) arguments for each 
exponent in the range 2**k < exponent < IHEXIS 
2**(k + 1) - 1, where k has integral a T777-79--T------- Tocco T7----- 1 
values from 1 to 9 inclusive. | | 30 | 40 | 50 | 65 | 75 | 
---}------ ł------- ł------- ł------- ł------ 1 
IHEXIS | a | -104 | 23 | 29 | 10.7 | 7.6 | 
PO POR 1 | b | 701 | 176 | 56 | 14.3 | 8.7 | 
| R.M.S. relative | Maximum relative | | c | 342 | 90 | 26 | 5.7 | 3.2 | 
| error/exponent |  error/exponent | | af] 552 | 171 | 57 { 18.6 | 12.7 | 
| *105*6 | *10**6 | Cee eee Locos dc ccce Lao 5 A J 
p-_--__---=-------- += 
| 0.00871 | 0.692 | IHEXIL 
Medeae eei emen SSeS x repo pos nea qorana pesce 1 
| | 30 | 40 | 50 | 65 | 75 | 
IHEXIL p---t------ ł------- ł------- ł------- ł------ 1 
p------------------  q7-7-7--------------- 1 | a |-1535 |-322 | 0.7 | 4.1 | 3.7 | 
| R.M.S. relative | Maximum relative | | b | 1441 | 355 | 73 | 17.5 | 10.7 | 
error/exponent |  error/exponent | | c | 1082 | 269 | 42 | 8.9 | 5.2 | 
| *10**15 | *10**15 | | a'| 1055 | 180 | 78 | 20.0 | 11.9 | 
— À——»——  "—— I------------------ A AS Sede ee A secca 
| 0.0995 | 1.73 | 
T m dan ceci da J 


Other Information: 


and Exceptional Conditions: IHEXIS: For large exponents, tor example, 


those greater than 1023, it is general- 


P: x = 0 with n s 0 ly faster and more accurate to use the 
module IHEXXS rather than IHEXIS, pass- 

I : OVERFLOW, UNDERFLOW ing the exponent as a floating-point 
Since x**(-m), where m is a positive argument. However, it should be noted 
integer, is evaluated as 1/(x**m), that IHEXXS will not accept a negative 

the OVERFLOW condition may occur when first argument, and thus it is neces- 


m is large, and the UNDERFLOW 
tion when x is very small. 


condi- sary to pass the absolute value of this 
argument, and also, in cases where the 
exponent is odd, to test the sign of 
the argument in order to be able to 
attach the correct sign to the numer- 
ical result returned. 


Implementation: 


e Module size: IHEXIS 152 bytes 


IHEXIL 152 bytes 


e Execution time: 
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Let M = number of significant bits in 
the exponent General Floating-Point Exponentiation 
N = number of 1 bits in the expo- 
nent Module Names and Entry Points: 
Module Entry 
Arqument name point 
Short float  IHEXXS IHEXXSO 
Long float IHEXXL IHEXXLO 
Function: 
TO calculate x**y, where x and y are 


floating-point numbers. 


Method: 


When 


if y > 0, and an error message if y < 


x = 0, the result x**y = 0 is given 


0. 


When x # 0 and y = 0, the result x**y = 1 


is given. Otherwise x**y is computed as. 
EXP (y*LOG(x)), using the appropriate 
mathematical function routines. 
Error and Exceptional Conditions: 
P: x = 0 with y s 0 
O : a. X < 0 with y + 0 : error caused in 
LOG routine 
b. y*LOG(x) > 174.673: error caused 
in EXP routine l 
Implementation: 
e Module size:  IHEXXS: 144 bytes 
IHEXXL: 152 bytes 


e Execution times: 


Approximate execution 


times in micro- 


seconds for the System/360 models given 


below are oktained from the table: 


ESI NACER = E MAA Terra Ween te te 
[Module | 30 | 40 | 50 | 65 | 75 
| Name | | | | | 

------- f------ $------ +------ +-----4---- 
| IHEXXS | 9809 | 2861 | 902 | 236 | 143 
-------}------ +------ +------ 4-----4---- 
|IHEXXL |30444 | 7453 | 1579 | 358 | 203 
——— Ls Li PAE AA Ba o o 


Shift-and-assign, Shift-and-load (fixed 


decimal) 


Module Name: IHEAPD 


Entry Points: 


Entry 
Operation point 
Shift and assign IHEAPDA 
Shift and load IHEAPDB 
Function: 
IHEAPDA: To convert a real fixed decimal 


number with precision (£4,q3,) to preci- 


sion (E2,d2), where Pa £ 31 and 
Fa < 15. 
IHEAPDB: To convert a real fixed decimal 


number with precision (p1,q14) to preci- 


sion (31,92), where p, < 31. 


Chapter 2. 


Method: 


The argument scale factor is subtracted 
from the target scale factor. The argu- 
ment is converted to precision 31 in a 
field with a shift equal to the magnitude 
of the difference between the scale fac- 
tors; the shift is to the left if the 
difference is positive and to the right 
if negative. 


If entry point IHEAPDB is used, the field 
is moved unchanged to the target. If 
entry point  IHEAPDA is used, the result 
is checked for  FIXEDOVERFLOW and then 
assigned to the target with the specified 
precision. The assignment may cause the 
SIZE condition to be raised. 


Error and Exceptional Conditions: 


I : FIXEDOVERFLOW or SIZE 


Implementation: 


e Module size: 360 bytes 


e Execution times: 


Let S = (Target scale factor) 
(argument scale factor) 
fa = 0 if S 2 31 
= 1 if S < 31 
fa 0 if ABS(S) 2 p 


1 if ABS(S) < p 
(where p is the precision of 
the argument) 
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Then the approximate execution times in 
microseconds for the System/360 models 
below are obtained 


given 


lowing formulas: 


IHEAPDA: 
(i) S > 
(ii) S = 
(iii) S « 
IHEAPDB: 
(i) S > 
(ii) S = 
(iii) S < 
(errs qe ==. 
| | 30 
E------ ł------ 
| ti | 1314 
ļ------ 4------ 
| ta | 451 
[------ i------ 
| ta | 1117 
ļ------ form 
| ta | 1220 
p------ +------ 
| ts | 502 
p----- ł------ 
| te | 1996 
t----~~ ł------ 
| t | 1799 
ļ------ $------ 
| te | 1902 
ottica dull 


P—+4+++-+-+-+—+—4 


te + tatf 
t7 


t *f 


ct 
Pp 


ta*f 


ct 
Uu 


ct 
E 


PR +—+— ee a 


COMPLEX. OPERATIONS 


i 


2 


a 


from the fol- 


Multiplication/Division (fixed binary) 


Module Name: 
Entry Points: 


Mathematical 
Operation 


Za*Z2 
Z1/Z2 


Function: 


To calculate z4,*zz2 Or 21/22, where z4 and 


IHEMZU 


Entry 
point 


IHEMZUM 
IHEMZUD 


Za are fixed-point binary 


ters. 


24 


complex num- 


Method: 
Let 2, = a + bl and za = c + dI. Then, 
for multiplication, an incorporated sub- 


routine is used to compute a*c - b*d and 
b*c + a*d; these are tested for FIXED- 
OVERFLOW and then stored as the real and 
imaginary parts of the result. 

used to 


For division, the subroutine is 


compute a*žc + b*d and b*c - a*d. The 
expression c**2 + d**2 is computed and 
the real and imaginary parts of the 
result are then obtained by division. 
The subroutine computes the expressions 
u*x + v*y and v*x - u*y. 
Error and Exceptional Conditions: 
I FIXEDOVERFLOW in either routine, 
ZERODIVIDE in the division routine. 
Implementation: 


e Module size: 240 bytes 
e Execution times: 


Approximate execution times in micro- 
seconds for the System/360 models given 


below are obtained from the following 
formulas and tables: 
[crc pee qoom a meme quem xsona 1 
[Entry | 30 | 40 | 50 | 65 | 75 | 
| Point] | | | | | 
~------}------ A 
| IHEMZUM| 2421 | 689 | 256 | 56.6| 37.1 | 
SARI Lisa La Loomis Less line 4 
IHEMZUD 
Let M = number of significant bits in 
Z2 *CONJG (Z2) 
N = FLOOR(M/4 - 8) 


Then the approximate execution times in 
microseconds for the System/360 models 
given below are obtained from the fol- 
lowing formula: 


a + b + c + N*d 

fr Tors qocc-—-—-— TS > Torre Ts ott 1 
| | 30 | 40 | 50 | 65 | 75 | 
ļ---4------ | Yo +------- +------ 1 
| a | 3021 |1340 | 420 { 94.1 | 64.2 | 
| b | 338 | 79 | 30 { 5.7 | 2.7 | 
| coc] 78 | 24 | 18 | 6.6 | 2.1 | 
| d | 213 | 56 | 22 | 6.2 | 3.8 | 
p2lcEiclooro bosses boone Looe sek deli 4 
Note: b= 0 if M < 31 

c= d = 0 if M< 32 


Multiplication/Division (fixed decimal) 


Module Name: IHEMZV 
Entry Points: 
Mathematical Entry 
Operation point 
Za*Z2 IHEMZVM 
Za/Z2 IHEMZVD 
Function: 


To calculate Z,*z2 Or Z4/22 where z, and 


Za are fixed-point decimal complex num- 
bers. 
Method: 

Let 2, = a + bI and Zzg = c + di. The 
products a*c, b*c, a*d and b*d are com- 
puted. Then the required result is 
oktained as follows: 
Multiplication: 

Real part a*c - b*d 

Imaginary part b*c + a*d 
Division: 

Real part (a*c + b*d)/(c*c + d*d) 


Imaginary part (b*c - a*d)/(c*c + d*d) 


Error and Exceptional Conditions: 


I : FIXEDOVERFLOW in either routine, 
ZERODIVIDE in the division routine. 
Implementation: 


e Module size: 672 bytes 


e Execution time: 


Let (p,q), (r,s) = precisions of the 
operands 
L, = FLOOR(p/2) + 1 
La = FLOOR(r/2) + 1 
T = c**2 + d**2 2 10**(15 - 2*s) 


(i.e., T = 1 or 0 depending on 
whether the relation is true 
or false) 
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Then approximate execution time 
microseconds 
given below are obtained from the 


lowing formulas: 


S 


in 


for the System/360 models 
fol- 


IHEMZVM: a + b*L4 + C*L2 + d*L4*La 
IHEMZVD: et f*L4 + g*La + h*L4*La * 
j*La**2 + T*(k + m*L5) 

CERERE pa a qr e NE Wo. ge Neo ee TU 1 
| | 30 | 40 | 50 | 65 | 75 | 
pot rn po Jr pra 1 
|a | 2605 | 956 | 457 | 128  |103 | 
| b | 246 | 105 | 49 | 17 | 13 | 
ic] 30]| 33 | 19 | 5.3] 3.5 | 
fa] 112| 15 | 0 | 4.0 | 4.7 | 
| e [14525 |4203 |1144 | 626 |508 | 
| £ | 246 | 105 | 49 | 17 | 12 | 
ig | 168 | 103 | 41 | 16 [11 | 
|h | 112| 15 | 0 | 54.0 | 4.7 | 
lj! 56] 7.5] 0 | 2.0] 2.3 | 
| k | 731 | 251 [131 | 41.2 | 30.2 | 
{mf 12] 16 | 4.3 | 0.8 | 0.3] 
E A Lowen e o drm A J 


Other Information: 


It should be noted from the timings for 
multiplication that where the operands 
differ in precision, it is faster to 


present the longer operand as the 
argument rather than the first. 


Multiplication (floating-point) 


Module Names and Entry Points: 


second 


Module Entry 
Argument name point 
Short float IHEMZW IHEMZWO 
Long float IHEMZZ IHEMZZO 
Function: 
To compute z4*za4 in floating-point, when 


Za = a + bI and za = c + dI. 


Method: 
The real and imaginary parts of 
result are computed as a*c - b*d 


b*c + a*d, respectively. 


Error and Exceptional Conditions: 


I : Exponent OVERFLOW and UNDERFLOW 
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and 
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Implementation: 


IHEMZW 64 bytes 
IHEMZZ 64 bytes 


e Module size: 


e Execution times: 


Approximate execution time in micro- 
seconds for the System/360 models given 
below are oktained from the table: 


pS EOS Tn TAS E AS TIOS 1 
[Module | 30 | 40 | 50 | 65 | 75 | 
| Name | | | | | | 
------- £------4------4------d-----1----- 
|IHEMZW | 1979 | 550 | 172 | 41.9] 23.3 | 
nn T--7--MEÉIM—MEÉ-—L-T 
|IHEMZZ | 5115 | 1307 | 251 | 62.3] 31.3 | 
Lll SE Lasse Ladera Leone ds or Looe J 
Division (floating-point) 
Module Names and Entry Points: 
Module Entry 

Argument name Point 

Short float IHEDZW IHEDZWO 

Long float IHEDZZ IHEDZ2Z0 
Function: 


To compute Z4/Za in floating-point, when 
Z4 = a + bI and za = c + dI. 
Method: 
1.  ABS(c) 2 ABS(d) 
Compute q = d/c 


then REAL (2,/z223) = 
IMAG (Z1/Za2) = 


(a + b*q)/(c + d*q) 
(b - a¥q)/(c + d¥q) 


2. ABS(c) < ABS(d) 


(a + bI)/(c + dI) = (b - al)/(d - cl), 
which reduces to the first case. 


The comparison between ABS(c) and ABS(d) 
is adequately performed in short preci- 
sion in both modules. 


Error and Exceptional Conditions: 


I : OVERFLOW, UNDERFLOW and ZERODIVIDE 
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Implementation: 


IHEDZW 
IHEDZZ 


104 bytes 
104 bytes 


e Module size: 


e Execution times: 


Approximate execution times in micro- 
seconds for the System/360 models given 
below are obtained from the table: 


ER er E mu aia ae ae 4. a ER TA VEIT CE We» vv. eed 
|Module | 30 | 40 | 50 | 65 | 75 | 
| Name | | | | | | 
------- 442-4222441 
|IHEDZW | 3546] 875 | 221 | 60.8| 35.7 | 
------- ł------ł------ł------}-----}------ 
|IHEDZZ | 11741] 2515 | 234 | 92.5| 51.1 | 
Coo Li SS: lA Le ee ees List ål EA d 


Positive Integer Exponentiation (fixed 
binary) 


Module Name: IHEXIU 


Entry Point: IHEXIUO 
Functiori: 


Tô calculate z**n, where n is a positive 
integer less than 2**31. 


Method: 


The contents of the target field are set 


to the value of z. The final result is 
obtained by repeated squaring of the 
contents of the target field or squaring 


and multiplying by z. Multiplication is 
performed by the complex multiplication 
routine IHEMZU. 


Range: 
Q0 < n < 2**31. 


The precision rules of PL/I impose a 
further restriction in that if z has a 
precision (p, q), this module may only be 
called if n*(p + 1) - 1 < 31. This 
implies that n < 32/(p + 1) < 16 for all 
such cases. 


Implementation: 
e Module size: 120 bytes 


e Execution time: 


Let M= number of significant bits in 
the exponent 
N = number of 1 bits in the expo- 


nent 


Then the approximate execution times in Then the approximate execution times in 


microseconds for the System/360 models microseconds for the System/360 models 

given below are obtained from the fol- given below are obtained from the fol- 

lowing formula: lowing formula: 

-a + b*M + c*N -a * b*M * c*N 

pog qom queens cmm qme mem quee 1 pescarese qc quum que s 1 
| | 30 | 40 | 50 | 65 | 75 | | | 30 | 40 | 50 | 65 | 75 | 
O ł------- ł------- ł------- ł------4 eo ł------- hor j------- }------ 1 
| a | 4169 | 1183 | 405 | 101 | 69.5 | | a | 9200 | 3100 | 1300 | 450 | 370 | 
| b | 2409 | 822 | 306 1 73.3 | 48.3 | | b | 5000 | 1700 | 730 | 250 | 200 | 
| c | 2553 | 738 | 276 | 62.7 | 42.1 | | c | 5500 | 1800 | 750 | 250 | 200 | 
Lal. dl es door tas J Esla. poer kassens lanii Lozza 4 


Integer Exponentiation (floating-point) 
Positive Integer Exponentiation (fixed 


decimal) 
Module Names and Entry Points: 
Module Name: IHEXIV 
Module Entry 
Entry Point: IHEXIVO Argument name point 
Short float IHEXIW IHEXIWO 
Function: Long float IHEXIZ IHEXIZO 


TO calculate z**n, where n is a positive Function: 

integer less than 2**31. 
To calculate z**n, where n is an integer 
between -2**31 and 2**31 - 1 inclusive. 


Method: 
Method: 
The contents of the target field are set 
to the value of the argument. The final If the exponent is 0 and the argument 
result is obtained by repeated squaring non-zero, the answer 1 is returned 
of the contents of the target field or immediately. If the exponent is non- 
squaring and multiplying ky the argument. zero, the contents of the target field 
Multiplication is performed by the com- are set to the argument value. The 
plex multiplication routine IHEMZV. exponent is made positive and the 


argument raised to this positive power by 
repeated squaring of the contents of the 


Range: target field or squaring and multiplying 
by the argument. Multiplication is per- 
The precision rules of PL/I impose the formed by a branch to the complex multi- 
restriction that if z has a precision plication routine. Then, if the exponent 
(p,q), this module may only be called if was negative, the reciprocal of the 
n*(p*1)-1s 15. This implies that result is taken, otherwise it is left 
ns 16/(p + 1) < 8 for all such cases unchanged. 
and, in fact, this module will operate 
only for the range 0 <n s 8. Error and Exceptional Conditions: 
P: z= 0 with n < 0 
Implementation: 
I : OVERFLOW, UNDERFLOW 
e Module size: 192 bytes Since x**(-m), where m is a positive 
integer, is evaluated as 1/(x**m), 
e Execution time: the OVERFLOW condition may occur when 
m is large and the UNDERFLOW condi- 
Let M = number of significant bits in tion when x is very small. 


the exponent 
H : OVERFLOW or UNDERFLOW in complex 
N = number of 1 bits in the expo- multiplication routine  (IHEMZW or 
nent IHEMZZ) 
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Implementation: 


IHEXIW 
IHEXIZ 


e Module size: 256 bytes 


256 bytes 
® Execution time: 


Let M = number of significant bits in 
the exponent 


N = number of 1 bits in the expo- 
nent 


Then the approximate execution times in 
microseconds for the System/360 models 
given below are obtained from the fol- 
lowing formulas: 


-a * b*M * c*N for positive exponents 


-a'+ b*M + c*N for negative exponents 


IHEXIW 
A qmm n quem Looe 1 
| | 30 | 40 | 50 | 65 | 75 | 
[---+------ +------- t------- 4------- 4------ 
| a | 1446 | 582 | 128 | 35.2 | 17.6 | 
| kt | 2125 | 565 | 174 | 45.8 | 26.9 | 
| c | 1782 | 484 | 147 | 36.1 | 21.1 | 
| | | | | | | 
| a'] -142 | -31 | -12 | 4.6 | -6.3 | 
EOE EES dolo cn dota ss. ue J 
IHEXIZ 
(Ss eo quee Treno Tem Tc 1 
| | 30 | 40 | 50 | 65 | 75 | 
}---+------ ł------- ł------- ł------- ł------ 1 
| a | 8524 | 2042 | 275 | 60.0 | 30.8 | 
| b | 5393 | 1397 | 289 | 59.9 | 34.9 | 
| c | 4918 | 1245 | 226 | 50.2 | 29.1 | 
| | | | | | | 
| a'| 95 | -374 | -53 | -5.1 | -4.8 | 
Llc enm d ciem mee Lor illa rta 4 
General Floating-Point Exponentiation 
Module Names and Entry Points: 
Module Entry 

Argument name point 

Short float IHEXXW IHEXXWO 

Long float IHEXXZ IHEXXZO 
Function: 

To calculate z,**Zz., where z, and Z2 are 


28 


complex numbers of the same precision. 


Method: 


When 2, = 0, the result 0 is returned if 
REAL(Z2) > 0 and IMAG(z,) = 0. 
Otherwise, z,**za4 is computed as 


EXP (z2*LOG(z,)), 


with the proviso that if IMAG(z,) = 0 
then LOG(ABS(z,)) is calculated by a call 
to the real LOG routine, not to the 


complex LOG routine. 


Error and Exceptional Conditions: 


P : Z4 = 0 with 
IMAG(z2) * 0 


either REAL(z2) < 0 or 


O : a. REAL(Z2*LOG(z,)) > 174.673: error 
caused in IHEEXS or IHEEXL 
b. IHEXXW: 
ABS (IMAG (z2*LOG(z,))) = 2**18*pi: 
error caused in SIN routine 
(IHESNS) 
IHEXXZ: 


ABS(IMAG(za*LOG(z4))) 2 2**50*pi: 
error caused in SIN routine 
(IHESNL) 


Implementation: 


e Module size: 


e Execution ti 


Approximate 
seconds for 
below are 


IHEXXW 280 bytes 
IHEXXZ 280 bytes 


mes: 


execut 


ion 


times in micro- 
System/360 models 
oktained from the table: 


IMAG (22) 


a = IMAG(z,) k = 

c = REAL(Z,) d = REAL(Zz2) 
| ana SS qme po. fåes 1 
| | 30 | 40 | 50 | 65 | 75 | 
RA ote AMES Eonia Iosu he coset EDER 4 
IHEXXW 
pos. ca as in en t ciii spin n e tn iu Pos seep anims mem 1 
ja = 0 | 20606| 5834| 1816 | 480| 291 | 
jc > 0 | | | | | | 
p------- +------ ł------ ł------ +----- ł------ 1 
ja = 0 | 21750| 6171| 1929 | 509| 311 | 
jd <0 | | | | | | 
[------- ł------ ł------ qe +----- +------ i 
ja #0 | 27448] 8022| 2414 | 687| 417 | 
|b= 0 | | | | | | 
poseecces pee ł------ {------ ł----- ł------ 4 
|a #0 | 28263] 8229| 2576 | 711] 424 | 
|b# 0 | | | | | | 
eee Ho ton Perde bott E RE beet J 
LHEXXZ 
[^69 == dad ai > ra iaa 1 
|a = 0 | 62440| 15325| 3206 | 745| 426 | 
Ic > 0 | | | | | | 
ļ------- ł------ ł------ ł------ +-----4------ 1 
ja = 0 | 65208| 16062| 3366 | 781| 450 | 
la < 0 | | | E | | 
}------- }------ }------ +------ +-----4------ 4 
ja #0 | 90418] 21611| 4524 | 1056| 604 | 
[b= 0 | | | | | | 
}------- +------ +------ +------ +-----4------ 4 
|a #0 | 92809| 22197| 4623 | 1077| 616 | 
|b# 0 | | | | | | 
El. go uc Lor Y MEET EN Ps Sot A ES J 
FUNCTIONS WITH REAL ARGUMENTS 
ADD (Fixed decimal) 


Module Name: i 


Entry Point: I 


Function: 


ADD (xX4,Xa2,r,q) 


HEADD 


HEADDO 


where 


x, and xg are real 
fixed-point decimal numbers, 


the required precision of the result. 


Chapter 2. 


given 


and (p,q) is 


Method: 


If both arguments are non-zero, a call to 
the module IHEAPD is used to shift the 
one with the larger scale factor to give 


it the scale factor of the other, and 
convert it to precision 31. The argu- 
ments are added together, and  IHEAPD is 


used to convert the sum to the specified 
precision and to assign it to the target 
field. 


If one of the arguments is zero, the 
other is treated as the sum above. 
Error and Exceptional Conditions: 
H : FIXEDOVERFLOW or SIZE may occur in 
IHEAPD. 
Implementation: 


e Module size: 216 bytes 


e Execution time: 


Let S4 time for IHEAPDA with argument 
equal to the sum of the two 
arguments to IHEADD, precision 
31 and scale factor equal to 
the minimum of the scale fac- 


tors of the two arguments 


time for IHEAPDB with argument 


equal to the argument to 
IHEADD with the larger scale 
factor 


Then the approximate execution time in 
microseconds for the System/360 models 
given below are obtained from the fol- 
lowing formulas: 


(i) Both arguments non-zero: 


tı + S4 + S3 


(ii) At least one argument zero: 
ta + Sa 

[romeo porum ESA Toena T pete 1 
| 30 | 40 | 50 | 65 | 75 | 
}------+------ 4-----4------- +------ +------ 4 
ta | 1895 | 648 | 320 | 85.0 | 65.6 | 
ļ------ ł------ ł-----ł------- ł------ +------ 4 
ta | 1619 | 534 | 260 | 71.5 | 53.0 | 
Leone sarao doce D d rc LE J 
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MAX, MIN 


Module Names and 


Argument 


Fixed binary 


Fixed decimal 


Short float 


Long float 


Function: 


Entry Points: 


PL/I Module 
function name 
MAX IHEMXB 
MIN 
MAX IHEMXD 
MIN 
MAX IHEMXS 
MIN 
MAX IHEMXL 
MIN 


Entry 
point 


IHEMXBX 
IHEMXBN 


IHEMXDX 
IHEMXDN 


IHEMXSX 
IHEMXSN 


IHEMXLX 
IHEMXLN 


To find the maximum or the minimum of a 
group of arithmetic values. 


All arguments 


Method: 


must have the same base, 
Scale and precision. 


IHEMXB,IHEMXS,IHEMXL: 
current maximum or minimum is set to 
the value of the first argument; it is 
then compared algekraically with the 


next argument and replaced by 


appropriate. 


The 


until a test on 

indicates that all 
been processed, when the current value 
is stored as the result. 


The value 


of the 


it if 


process is repeated 
the argument list 
source items have 


IHEMXD: The address of the current maxi- 
mum or minimum is set to the address of 
the first argument; 
then compared algebraically with the 


next argument, and the address 
latter replaces 


appropriate. 


keen processed, 


this argument is 


The process is 
until a test on 
indicates that all 


of the 


that of the former if 
repeated 


the argument list 
source items have 


when the 


moved into the target field. 


Implementation: 


30 


e Module sizes: 


IHEMXB 96 bytes 
IHEMXD 120 bytes 
IHEMXS 96 bytes 
IHEMXL 96 bytes 


result is 


Execution time: 


IHEMXB 
Let N = the number of source arguments 


Then the average execution times in 
microseconds for the System/360 models 
given below are obtained from the fol- 
lowing formulas: 


IHEMXBX 
N 
a + b*N + ct DD (17i) 
i=2 
IHEMXBN 
N 
a' + b*N + c* 2, (1/i) 
i-2 
Thee Ter Se querere quee fræse 7 
| 30 | 40 | 50 | 65 | 75 | 
f------ t------- ¢-------}------- $-----+4 
| 391 | 135.7 | 51.8 | 12.1 | 7.4 | 
| 186 | 70.7 | 30.0 | 11.0 | 8.06| 
| 54 | 21.3 | 8.0 | 2.5 | 0.9 | 
| | | | | | 
| 356 | 123.8 | 47.3 | 10.7 | 6.3 | 
Ellos Las E E AAA desta J 
IHEMXD 
Let N = the number of source arguments 


(p,q) = the precision 


L = FLOOR(p/2) + 1 


Then the average execution times in 
microseconds for the System 360 models 
given below are obtained from the fol- 
lowing formulas: 


IHEMXDX 


N 
a + (b + d*L*N + c* DD (1/1) 


i=2 
IHEMXDN 
N 
at + (b + d*L)*N + c* © (1/i) 
i=2 
i BER ES era ow Wa gars a DEO E ER quei qoo PTT. 1 
| 30 | 40 | 50 | 65 | 75 | 
quee ee qe pesata === 1 
| 670 | 217 | 84.1 | 20.9 | 12.3 | 
| 232 | 92.5 | 46.4 | 16.4 | 13.5 | 
| 44 | 16.9 | 65| 1.9 | 1.5 | 
| 5| 2.5] 2.1 [| 0.4 | 0.4 | 
| | | | | | 
| 635 | 205 | 79.6 | 19.6 | 11.3 | 
dee, dee p doe dimen d 


IHEMXS, IHEMXL 


Let N = the number of source arguments 


Then the average execution times in 
microseconds for the System/360 models 
given below are obtained from the fol- 
lowing formulas: 
IHEMXSX, IHEMXLX 
N 
a + b*N + c* 2, (1/1) 
i=2 
IHEMXSN, IHEMXIN 
N 
at + E*N + ct i’ (17i) 
i=2 
IHEMXS 
ENS TR RR qose quum qon eee TARO 1 
| | 30 | 40 | 50 | 65 | 75 | 
AS $------- j------- 4------- ł------ 1 
| a 367 | 131.2 | 44.6 | 11.5 | 7.3 | 
| b {| 219 | 73.3 | 29.1 | 10.5 | 7.7 | 
| c | 55 | 21.3 | 7.3 | 2.5 | 1.9 | 
| | | | | | | 
| a"f 332 | 118.6 | 40.1 | 10.2 | 6.2 | 
oo tee cece Bacca delusa dicci Lara J 
IHEMXL 
ESAS” A RARAS Tess ad puce 1 
| | 30 | 40 | 50 | 65 | 75 | 
ļł---ł------ IS e a 4 
| a 367 | 133.2 | 46.4 | 11.7 | 7.3 | 
| b | 251 | 81.2 | 31.4 | 10.6 | 7.71| 
| cl 71 | 26.3 | 9.3 | 2.7 | 1.9 | 
| | | | | | | 
|a'| 332 | 120.6 | 41.9 | 10.4 | 6.2 | 
A AA Lino tass La Lope deluse J 


FUNCTIONS WITH COMPLEX ARGUMENTS 


ADD (Fixed decimal) 


Module Name: IHEADV 

Entry Point: IHEADVO 

Function: 
ADD (Za,ZaePsq) where Z4 and zg are com- 
plex fixed-point decimal numbers, and 
(p,q) is the required precision of the 
result. 

Method: 


The real parts of each argument are added 
and the sum is assigned to the target 
field by using the real fixed decimal ADD 
module (IHEADD). The imaginary parts are 
treated similarly. i 


Chapter 2. 


Error and Exceptional Conditions: . 


H : FIXEDOVERFLOW or 
IHEAPD. 


SIZE may occur in 


Implementation: 
e Module size: 96 bytes 
e Execution time: 


Let T4, = execution time for IHEADD with 
the real parts of the argu- 
ments for IHEADV as arguments 


To = execution time for IHEADD with 
the imaginary parts of the 
arguments for IHEADV as  argu- 
ments 


Then the approximate execution times in 
microseconds for the System/360 models 
given below are obtained from the fol- 
lowing formula: 


-——— —À I ee er 


MULTIPLY (fixed binary) 


Module Name: IHEMPU 
Entry Point: IHEMPUO 
Function: 
MULTIPLY (z4,2Z2,P,q) where z, and zz are 


complex fixed-point binary numbers, and 
(p,q) is the required precision of the 
result. 


Method: 


Let the arguments be z, = a + bI and 
Za 7 C + dI. 


Then REAL(Z24*za4) = a*c - b*d 

IMAG (24 *Z2) = b*c + a*d 
The real and imaginary parts of the 
product are computed. These numbers are 
then shifted to give them the required 
scale factor(q). 


The results of the shifts are tested for 
FIXEDOVERFLOW and truncated by left 
shifts. 
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Error and Exceptional Conditions: 
I : FIXEDOVERFLOW 
Implementation: 
e Module size: 240 bytes 


e Execution time: 


Implementation: 


e Module size: 288 bytes 
e Execution time: 


Let (pi,da)(pa,qd2) = the precisions of 
the arguments: 


Li = FLOOR(p,/2) 

T = Time to shift-and-assign the 
result using  IHEAPDA, with 
argument of precision 
(pa + Pa + 1,91 + 92) and 


result (p,q) 


Let gi = scale factor of the first 
argument 
G2 = scale factor of the second 
argument 
Q = scale factor of the target 
Then the approximate execution times in 


microseconds for the System/360 models 


given below are obtained from the 
table: 
k = qı + q2 
[ces TS per pesare qr queo TT DOT 1 
| | 30 | 40 | 50 | 65 | 75 | 
E------- ho $------ O i------ 4 
| Q= k | 2816 | 847 | 329 | 75.7 | 49.3 | 
l------- ł------ ł------ ł----- ł------ ł------ 1 
| Q>k | 3098 | 910 | 354 | 81.5 | 54.9 | 
[------- ł------ ł------ ł----- Ho 4------ 1 
| Q < k | 3394 | 1052 | 404 | 96.0 | 66.8 | 
RESTI A Tissi E dee luca 4 
MULTIPLY (fixed decimal) 
Module Name: IHEMPV 
Entry Point: IHEMPVO 
Function: 
MULTIPLY(zi,Z2,p,q) where z, and za are 
complex fixed-point decimal numbers, and 
(p,q) is the required precision of the 
result. 
Method: 
Let z = a + LI and Za = c + dI, then: 
REAL(24*23) = a*c - b*d. 
IMAG(z4*za3) = b¥c + a*d, 
The real and imaginary parts are calcu- 


lated and then each is assigned to the 
target with precision (p,q) by separate 
calls to the entry point IHEAPDA of the 
decimal shift and assign module IHEAPD. 


Error and Exceptional Conditions: 
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H : FIXEDOVERFLOW or SIZE in IHEAPD. 


Then the appro 
microseconds f 
given below a 
lowing formula 


2*T + a + 


ESISTERE TREE 
|| 30 | 40 
}-+------ 4-------- 
|al 3164 | 1170 
}-+------ j-------- 
|b| 358 | 89.7 
t-+------ ł-------- 
[c| 142 | 18.4 
hat +-------- 
{aj 112 | 15.0 
Llc-:cc-- llutrtonzan 


ximate execution times in 
or the System/360 models 
re obtained from the fol- 


b*L, + C#*L, + d*L,*La 


DIVIDE (fixed binary) 


Module Name: IHE 


Entry Point: IHE 


Function: 


DIVIDE(Z,,Z2,Pr4 
complex fixed-p 
(p,q) is the req 
result. 


Method: 
Let Z, = a + bI, 


REAL (Z1/ Z2) = 
IMAG(z4/22) = 


The expressions 

C**2 + d**2 are 

of 63. The de 
shifted to preci 
or left shift. 


quem qon qp 1 
| 50 | 65 | 75 | 
$------- $------- 4------- 

| 540 | 168.1 | 131.7 | 
+------- +------- 4------- 

| 46.4 | 25.2 | 17.3 | 
}------- ł------- poo 

| 18.3 | 8.8 | 7.9 | 
j------- $------- ł------- 

| 0 | 4.0 | 4.7 | 
Lares rs Lats 4 
DVU 

DVUO 

) where z4 and Za are 


oint binary numbers, and 
uired precision of the 


and Za = c + dI, then: 


(atc + b*d)/(c**2 + d**2) 
(b*c - a*d)/(c**2 + d**2) 


a*c * b*d, b*c - a*d, and 
computed with a precision 
nominator, c**2 + d**2 is 
sion 31 by either a right 


Two calls are then made to an incorporat- 
ed subroutine which accepts a numerator 
and shifts it so that it has two insigni- 
ficant leading digits. It then divides 
Ly c**2 + d**2 and shifts the quotient to 
the required scale factor (q). 


Error and Exceptional Conditions: 


I : FIXEDOVERFLOW or ZERODIVIDE 


Implementation: 
e Module size: 408 bytes 


e Execution time: 


Let N4 number of significant bits in 


the expression a*c * b*d 


N2 = number of significant bits in 
the expression b*c - a*d 


Na = number of significant bits in 
the expression c**2 + d**2 
F4 = FLOOR ((61 - N,)/4) if N, < 61 
= 0 if N, > 61 
F2 = FLOOR ((61 - N2)/4) if No £ 61 
= 0 if N. > 61 
Fa = FLOOR ((Na - 32)/4) 
Sa = (61 - N4) = F,*4 if Ni < 61 
= 0 if N, > 61 
S2. = (61 = N) E” F2*4 if Na < 61 
= 0 if Na > 61 
Chapter 2. 


Then the approximate execution times in 
microseconds for the System/360 models 
given below are obtained from the fol- 
lowing formulas: 


(i) na > 32: ty + tat (Fy + Fa) + 
ta*(S, + So) + ty*Fs 


(ii) na = 32: ts + to*(F, + Fo) + 
ta*(S4 + S2) 


(iii) na « 32: te * ta*(F, + F2) + 
ta* (S4 + Sg) + t7*Na 


[77y2-2--7 T SSP ERER > ÓN Tosse Trenen 1 
| | 30 | 40 | 50 | 65 | 75 | 
porton fo ł------- ł------- ł------- 1 
|t4| 5924| 1835 | 609 | 144 | 104 | 
I--4----- ł-------- ł------- ł------- ł------- 1 
Ital 261| 73.5 | 26.0 | 7.6 | 5.1 | 
[--1----- ł-------- ł------- ł------- ł------- 1 
It3| 251] 73.5 | 25.5 | 7.8 | 5.0 | 
[--4----- ł-------- ł------- ł------- ł------- 1 
[tel 213] 55.6 | 21.5 | 6.0 | 3.8 | 
I--1----- ł-------- ł------- ł------- }------- 1 
|ts| 5714| 1780 | 582 | 138 | 100 | 
I--1----- ron ł------- ł------- ł------- 1 
|te! 3712| 1288 | 473 | 115 | 86.5 | 
I--1----- ł-------- ł------- ł------- ł------- 1 
jt7| -111| -38.8 | -13.3 | -3.4 | -2.6 | 
DIIcplibll d Ico T NND SETE SEE qc cles pcc cei. J 
DIVIDE (fixed decimal) 
Module Name: IHEDVV 
Entry Point: IHEDVVO 
Function: 

DIVIDE(Z1,Z2,P,q) where z, and Za are 


complex fixed-point decimal numbers, and 
(p,q) is the required precision of the 
result. 


Method: 
Let z, = a + bI, and za = c + dI, then 


REAL(z4/z2) = (a*c + b*d)/(c**2 + d**2) 
IMAG (24/22) = (b*c - a*d)/(c**2 + d**2) 


The expressions a*c + b*d, b*c - afd, and 
C**2 + d**2 are computed. Leading zeros 
are removed from the denominator 
(c**2 + d**2) by truncation on the left 
and a left shift if necessary. If the 
denominator is still more than 15 digits 
long it is truncated on the right to 15 
digits. 


Two calls are then made to an incorporat- 


ed subroutine which accepts a numerator 
and shifts it to precision 31 with 2 
leading zeros by calling IHEAPD 


(via 


Arithmetic: Functions with Complex Arguments 33 


entry point IHEAPDB). It then divides by 
C**2 + d**2 and calls IHEAPD (via entry 
point IHEAPDA) to assign the quotient to 
the target field with the required preci- 
sion (p,q). 
Error and Exceptional Conditions: 
I : ZERODIVIDE 
H : FIXEDOVERFLOW or SIZE in IHEAPD 
Implementation: 
e Module size: 576 bytes 
* Execution time: 


Let pa4,E2 = the precisions of the argu- 


ments 
L4, = FLOOR(p4/2) 
L = FLOOR(p2/2) 


Then the approximate execution times in 
microseconds for the System/360 models 
given below are oktained from the fol- 
lowing formula: 


a + b*L, + c*L, + d*L¿*L2 + e*L,**2 


uni see Wee TC qoi Y qd 2E quao cero 1 
| | 30 | 40 | 50 | 65 | 75 | 
A ł------- ł------ ł------- ł------ { 
| a | 28016 |7646 |3495 |1146 |975 | 
}---+------- +------- +------ +------- ł------ 1 
| t | 368 | 94.7 | 58.7| 25.8 | 17.9 | 
}---4------- ł------- ł------ ł------- ł------ 1 
| c | 382 | 67.5 | 60.4] 23.2 | 19.9 | 
}---4------- $------- ł------ ł------- ł------ 1 
[aj 112 | 15.0 | O | 4.0 | 4.7 | 
}---4------- ł------- ł------ ł------- ł------ 
| e | 56] 7.5 | 0 | 2.0 | 2.3 | 
Licda Lc es es ass d--cc- Loe dll J 


ABS (fixed binary) 


Module Name: IHEABU 


Entry Point: IHEABUO 
Function: 


TO calculate ABS(z) = SQRT(x**2 + y**2), 
where z = x + yl. 


Method: 
if X= y, result is x*SQRT(2). Other- 
wise, 
let X1 = MAX(ABS(x),ABS(y)) 
Y1 = 


MIN(ABS(x),ABS(y)). 
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Then ABS(z) is computed as 

X1*SORT(1 + (Y1/X1)**2), 
where the fixed binary calculation of 
SQRT(g) for 1 < g < 2 is included within 
the module. 


The first approximation to the square 
root is taken as 


g/ (ltg) + (1+g)/4, 

with maximum relative error 1.8*2**-10. 
One Newton-Raphson iteration gives maxi- 
mum relative error 1.6*2**-20, and suffi- 
ces if X1 « 2**(15-q) where q is the 
Scale factor of z. 
Otherwise a second iteration is used, 
with theoretical maximum relative error 
of 1.3*2**-40. 

Error and Exceptional Conditions: 
I : FIXEDOVERFLOW 

Implementation: 


e Module size: 184 bytes 


e Execution times: 


Approximate execution times in micro- 

seconds for System/360 models given 

below are obtained from the table: 

a = 2**(15-q) 

persona qucm o ES rene quere quum 1 
| xi | 30 | 40 | 50 | 65 | 75 | 
------- ł------ł}------ł------ł}-----4------4 
| <a | 3809 | 1218 | 320 | 79.4| 52.3 | 
ļ------- ł------ ł------ ł------ ł-----}------ 1 
| za | 4601 | 1473 | 372 | 93.1] 59.8 | 
tias uso KASSE HS een La | EP NE J 


ABS (fixed decimal) 


Module Name: IHEABV 


Entry Point: IHEABVO 
Function: 


TO calculate ABS(z) = SQRT (x**2 + y**2) 
where z = x + yl. 


Method: 
x and y are converted to binary, with 


appropriate scaling if either exceeds 9 
Significant decimal digits. 


Let X1 be the maximum, and Y1 the  mini- Then the approximate execution times in 

mum, of the absolute values of the two microseconds for the System/360 models 

binary numbers thus obtained. shown below are obtained from the fol- 
lowing formulas: 

Then if X1 = Y1 = 0, result 0 is 


returned. Otherwise, an approximation to L<5and D2< 7 : a 
ABS(z) is computed as 
L < 5, 7 < Di < 10 
X1*SORT(1 + (Y1/X1)**2), and Da > 7 : b+f*L+g*L**2 
where the fixed binary calculation of 5 < L s 8 and Da <7 : c 
SQRT (g) for 1s g < 2 is included within 
the module. 5 <Ls 8, 7 < Di 10 


< 
and Da > 7 : d*f*L*g*L**2 
The first approximation to the square 
root is taken in the form 5 < LS 8 
and 10 < Di 3515: e+f*L+tg*L**2 
A + B*(1 +g) - A/(1 + g) 


[reve poc mme TESS y mi pst ds qoc ccm Toe 1 
with maximum relative error 2.17*10**-4, | | 30 | 40 | 50 | 65 | 75 | 
and one Newton-Raphson iteration then I---1------ +------- +------- I------- Homo 4 
gives a value with maximum relative error | a | 6220 |1971 | 656 | 169 {116 | 
2.35*10**-8. | b [13001 |3785 {1101 | 460 [352 | 
| c | 6666 |2200 | 737 | 190 {132 | 
Muitiplication by X1 produces a value for | d {13447 (4014 [1182 | 481 1368 | 
ABS(z) which is rounded and converted to | e |13918 [4194 |1279 | 509 |391 | 
decimal, and this suffices if it has not | £ | 82 | 61.6 | 40.1 | 7.9 | 5.2 | 
more than 7 significant decimal digits. | gl 56 | 7.5 | 0.0 | 2.0 | 2.3 | 
Otherwise, this approximation is scaled Lib 1------- 1------- 1—------- l-.----- 4 
if necessary and used in a final Newton- 
Raphson iteration for SQRT(x**2 + y**2) 
in decimal, with theoretical maximum 
relative error 2.76*10**-16. 
Error and Exceptional Conditions: ABS (floating-point) 
I : FIXEDOVERFLOW 
Module Names and Entry Points: 
Implementation: Module Entry 
Argument name point 
e Module size: 544 bytes 
Short float IHEABW IHEABWO 
e Execution times: Long float IHEABZ IHEABZO 
Let (p,q) » the precision of the argu- 
rent . Function: 
L = CEIL((p+1)/2), i.e., the To calculate ABS(z) = SQRT (x**2 + y**2), 
length in bytes of each of the where z = X + yI. 
real and imaginary parts of 
the argument 
Method: 
Da = maximum number of significant 
digits in real and imaginary Let z = X + yI. If x = y= 0, answer is 
parts of the argument 0. 
D2 = number of significant digits Otherwise let Z1 = MAX (ABS(x),ABS(y) ) 
in result and 22 = MIN(ABS(x),ABS(y)). 


Then the answer is computed as 


Z1*SORT(1 + (22/21)**2). 
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Accuracy: 
IHEABW 
poe ree > er SEE 1 
| Arguments | Relative Error | 
| *10**6 | 
p A vr 
| Range {Distribution| RMS | Maximum | 
}~---------}------------ +------- ł--------- 1 
|Full range|Exponential | 0.833 | 2.02 | 
| radially, | | i 
| | uniform | | | 
| |round origin| | | 
Las Li IL Le sor arios J 
IHEABZ 

CEROTTI ES puce ee 1 
| Arguments | Relative Error | 
| *10**15 | 
Se ¡RUE Pana occu. presen 
| Range |Distribution| RMS | Maximum | 
}----------}------------ }------- de 4 
|Full range|Exponential | 0.828 | 3.38 | 
| | radially, | | | 
| | uniform | | | 
| |round origin| | | 
ese ae ee a ee Ls to Luo Lon J 
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Error and Exceptional Conditions: 


I : OVERFLOW 


Implementation: 


e Module size:  IHEABW 128 bytes 


IHEABZ 128 bytes 
e Execution times: 
Approximate execution times in micro- 


seconds for the System/360 models given 
below are obtained from the table: 


FEE STE po TAGS D esae p OG ERIS o v E NONSE 1 
[Module | 30 | 40 | 50 | 65 | 75 | 
| Name | | | | | | 
------- E 
|IHEABW | 5595| 1493 | 447 | 129 | 79.9| 
------- }------}------}------}-----}------ 
| IHEABZ | 14318] 3191 | 695 | 174 | 104 | 
to di RS Lo T ee ee doo... do J 


The Library supports all float arithmet- 
ic generic functions and has separate 
modules for short and long precision real 
arguments and also for short and long 
precision complex arguments where these are 
admissible. 


Since the calling sequence generated in 
compiled code is the same as that required 
for passing the same arguments to a PL/I 
procedure, it is permissible to pass the 
names of any of the float arithmetic gener- 
ic functions as arguments between proce- 
dures, according to the normal rules for 
entry names. 


Any restrictions on the admissibility of 


arguments are noted under the headings 
"Range" and ‘Error and Exceptional Condi- 
tions.' 

Range: This states any ranges of arguments 


which a module assumes to have been exclud- 
ed prior to its being called. 


Error and Exceptional Conditions: These 
cover conditions which may result from the 


use of a routine; they are listed in four 
categories: 


conditions in the module 
Programmed tests are 


P -- Prograrmed 
concerned. 


made where this is not too costly 
and, if an invalid argument is 
found, a branch is taken to the 


entry point IHEERRC of the execution 
error package (EXEP). This results 
in the printing of an appropriate 
message and in the ERROR condition 
being raised. 


I -- Interrupt conditions in the module 
concerned. For those routines where 
SIZE and FIXEDOVERFLOW are detected 
by programmed tests or where hard- 
ware interruptions may occur, the 
OVERFLOW and UNDERFLOW conditions 
pass to the ON handler (IHEERR) and 
are treated in the normal way. The 
machine is assumed to be enabled for 
all interruptions except signifi- 
cance, which is masked off. 
O -- Programmed conditions in modules 
called by the module concerned. 
These occur when invalid arguments 
are detected in the module called. 


H -- As I, but the interrupt conditions 
occur in the modules called by the 
module concerned. 
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Speed 


The average execution times given are 
based on the IBM System/360 Instruction 
Timing Information, Form  A22-6825. These 
times include times for the modules called. 


Accuracy 


In order to appreciate properly the 
meaning of the statistics for accuracy 
given with each module, some consideration 
of the limits and implications of these 
Statistics is required. Because the size 
of a machine word is limited, small errors 
may ke generated by mathematical routines. 
In an elaborate computation, slight inaccu- 


racies can accumulate and become large 
errors. Thus, in interpreting final 
results, errors introduced during the var- 


ious intermediate stages must be taken into 
account. 


The accuracy of an answer produced by a 
routine is influenced by two factors: (1) 
the accuracy of the argument and (2) the 
performance of the routine. 


Most arguments contain errors. ..An error 
in a given argument may have accumulated 
over several steps prior to the use of the 
routine. Even data fresh from input con- 
version may contain slight errors. The 
effect of an argument error on the accuracy 
of an answer depends solely on the nature 
of the mathematical function involved and 
not on the particular coding by which that 


function is computed within a routine. In 
order to assist users in assessing the 
accumulation of errors, a guide on the 


propagational effect of argument errors is 
provided for each function. Wherever pos- 
Sible, this is expressed as a simple formu- 
la. 


The performance statistics supplied in 
this document are based upon the assumption 
that the arguments are perfect (i.e., with- 
out errors, and therefore having no argu- 
ment error propagation effect upon 
answers). Thus the only errors in answers 
are those introduced by the routines them- 
selves. 


For each routine, accuracy figures are 
given for the valid argument range or for 
representative segments of this. In each 
case the particular statistics given are 
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those most meaningful to the function and 


range under consideration. 


For example, the maximum relative error 
and the root-mean-square of the relative 
error of a set of answers are generally 
useful and revealing statistics, but are 
useless for the range of a function where 
its value becomes 0, since the slightest 
error of the argument value can cause an 
unbounded fluctuation in the relative mag- 
nitude of the answer. Such is the case 
with SIN(x) for values of x close to pi; in 
this range it is more appropriate to dis- 
cuss absolute errors. 


derived from random 
distributions of 5000 arguments per  seg- 
ment, generated to be either uniform or 
exponential, as appropriate. It must be 
emphasized that each value quoted for the 
maximum error refers to a particular test 
using the method described akove, and 
should be treated only as a guide to the 
true maximum error. 


The results were 


This explains, for example, why it is 
possikle that the maximum error quoted for 
a segment may be greater than that found 


from a distribution of different arguments 
over a larger range which includes the 
former. 


Hexadecimal Truncation Errors 


While the use of hexadecimal numbers in 
System/360 has led to increased efficiency 
and flexibility, the effect of the variable 
number of significant digits carried by the 
floating-point registers must be noted in 
making aliowance for truncation errors. In 
the production of the PL/I Library, special 
care was taken to minimize such errors, 
whenever this could be accomplished at 
minor cost. As a result, the relative 
errors produced by some of the Library 
routines may ke considerakly smaller than 
the relative error produced in some instan- 
ces by a single operation such as multi- 
plication. 


Representations of finite length  entail 
truncation errors in any number system. 
With binary normalization, the effect of 
truncation is roughly uniform. With hexa- 
decimal normalization, however, the effect 
varies by a factor of 16 depending on the 
Size of the mantissa; in a chain of  compu- 
tations, the worst error committed in the 
chain usually prevails at the end. 


In short-precision 
number has between 21 
binary digits. 


representation, a 
and 24 significant 
Therefore, the truncation 
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errors range from 2**-2u to 2**-20 
(5.96*10**-8 to 9.5*10**-7). Assuming 
exact operands, a product or quotient is 
correct to the 24th binary digit of the 
Mantissa. Hence truncation errors contrib- 
uted by multiplication or division are no 
more than 2**-20. The same is true for the 
sum of two operands of the same sign. 
Subtraction, on the other hand, is the 


commonest cause of loss of significant 
digits in any number system. For short- 
precision operations, therefore, a guard 


digit is provided which helps to reduce 
such loss. 
In long-precision representation, a 


number has between 53 and 56 significant 
binary digits. Therefore truncation errors 
range from 2**-56 to 2**-52 (1.39*10**-17 
to 2.22*10**-16). Assuming exact operands, 
a quotient is correct to the 56th binary 
digit of the mantissa. Therefore,  trunca- 
tion errors resulting from division are no 
more than 2**-52. The accuracy of a pro- 
duct, on the other hand, depends on the 
necessity for post-normalization. If the 
mantissas of both operands are close to 1, 
the truncation error of a product is about 


2**-56. If the product of the mantissas is 
about 1/16, the truncation error is about 
2**-52. On the other hand, if the mantis- 


sas of both operands are close to 1/16, the 
intermediate product has 7 leading zeros, 
and post-normalization introduces 4  trail- 
ing zeros. In this case, the truncation 
error can be close to 2**-48 
(3.55*10**-15). In particular, multi- 
plication by 1 in the long-precision form 
has the effect of erasing the last hexa- 
decimal digit of the multiplicand. 


Normal care in numerical analysis should 
be exercised for addition and subtraction. 
In particular, when two algorithms are 
theoretically equivalent, it usually pays 
to choose the one which avoids subtraction 
between operands of similar size. There is 
no guard digit for long-precision additions 
and subtractions. 


Hexadecimal Constants 


Many of the modules described below 
discriminate between algorithms or test for 
errors by comparisons involving hexadecimal 
constants; it must be realized that where 
decimal fractions are used in the descrip- 
tions the fractions are only quoted as 
convenient approximations to the hexadeci- 
mal values actually employed. 


Terminology 


Maximum and root-mean-square values for 
the relative and (where necessary) the 
absolute errors are given for each module. 
These are defined thus: 


Let f(x) = the 
tion 


correct value for a func- 


g(x) = the result obtained from the 
module in question 


Then the absolute error of the result is 
ABS(f(x) - g(x)), 

and the relative error of the result is 
ABS ((£(x) - g(x))/£(x)). 


oktained 
of the 


results 
root-mean-square 


Let the number of sample 
be N; then the 
absolute error is 


SQRT (QU i (ABS (£ (Xi) - g(x4)) **2)/N), 
of the 


and the root-mean-square relative 


error is 


SORT (22; (ABS ((£(x4) - g (xq) )/£ (34) ) +42) /N). 


The Library mathematical modules are 
summarized in Figures 4 and 5. 


| 
a a a a a es we +----------- — Sa yess 
| SORT | IHESQS | IHESOL | 
EXP |  IHEEXS |  IHEEXI | 
| LOG, LOG2, LOG10 | IHELNS | IHELNL | 
SIN, COS,SIND,COSD| IHESNS | IHESNI | 
| TAN, TAND | IHETNS |  IHETNL | 
| ATAN, ATAND | IHEATS | IHEATI | 
|SINH, COSH |  IHESHS |  IHESHL | 
| TANH |  IHETHS |  IHETHI | 
| ATANH | IHEHTS |  IHEHTI | 
| ERF, ERFC |  IHEEFS | IHEEFL | 
EC. aw Ii aUa E m re eel Me MS ur Uu Lec ad 


Mathematical Functions with Real 
Argurents 


Figure 4. 
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4 
Short Float | Long Float | 


Function | 
laica nce rca pee 
| SORT |  IHESQW |  IHESQZ | 
| EXP |  IHEEXW |  IHEEXZ | 
| LOG |  IHELNW |  IHELNZ | 
|SIN,COS,SINH,COSH | IHESNW | IHESNZ | 
[TAN, TANH |  IHETNW | IHETNZ | 
| ATAN, ATANH | | IHEATW |  IHEATZ | 
sesta la cc c A, em metà 


Mathematical Functions with Com- 
plex Arguments 


Figure 5. 


FUNCTIONS WITH REAL ARGUMENTS 


SQRT (short floating-point real) 


Module Name: IHESQS 


Entry Point: IHESQSO 
Function: 

To calculate the square root of x. 
Method: 

If x = 0, SQRT(x) = 0. Otherwise, let 
X = 16**(2*p + q)*f, 


where p is an integer, 
1/16 < f < 1. Then 


q=0 or 1, and 


SORT (x) 
where z 
Z 


16** (p + q)*z, 
SQRT (f) if q 
SQRT (f)/4 if q 


1 


1. 


An initial approximation, Yo, is taken in 
the hyperbolic form a + b/(c + f) with 
different sets of constants for the two 
cases: 

1. q= 0 1.80713 
-1.57727 


a 
b 
c 0.954182 


The maximum relative error in this range 
is then less than 2**(-5.44), with an 
exact fit at f= 1 to guard as far as 
possible against loss of the last hexa- 
decimal digit when f is nearly 1. 
2. q=1 0.428795 
-0.3430368 
0.877552 


a 
b 
c 
The maximum relative error in this range 


is less than 2**(-6)*f**(-1/8). 
Then Ya = 16**(p + Dry. 
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Two Newton-Raphson iterations then yield: 


Ya = (ya + x/ya)72 
SQRT(xX) = ya + (x/Ya ~ Y2)/2 
For case q = 0, the final relative error 
from this algorithm is less than 


2**(-20.7), and, for case 


than 2**(-29). 


q = 1, less 


Effect of Argument Error: 
The relative error caused in the result 
is approximately half the relative error 
in the argurent. 


Accuracy: 


| Arguments Relative Error | 
*10**6 | 

---------- A OS 
| Range | | Maximum | 


Error and Exceptional Conditions: 
P: x < 0 
Implementation: 
* Module size: 168 bytes 
e Execution time: 
Approximate execution times in micro- 


seconds for the System/360 models given 
below are obtained from the table: 


ICONE ANE E Mia SA A ea ee UE XESG EN 7 5.2 — 
[Module | 30 | 40 | 50 | 65 | 75 | 
| Name | | | | | | 
------- dtt 
|IHESQS | 3140 | 793 | 227 | 68.4] 40.7 | 
Lowes Lis illc Llloleze de Lom J 


SQRT (long floating-roint real) 


Module Name: IHESQL 


Entry Point: IHESQLO 
Function: 


TO calculate the square root of x. 


Method: 
If x= 0, SQRT(x) = O. Otherwise, let 
x = 16**(2*p - q)*f, where p is an inte- 
ger, q = 0 or 1, and 1/16 f < 1. Then 


SORT(x) = 16**p*2**(-2*q)*SORT(£f). 


40 


An initial approximation, yo, is taken by 
using (2/9 + 8/9*f£) for SQRT(f). Multi- 
plication by 2**(-2) when q = 1 is accom- 
plished by using the  HALVE instruction. 
twice. The maximum relative error of 
this approximation is 1/9. 


Four Newton-Raphson iterations of the 
form Yn+1 = (yn + X/Yn)/2 are then 
applied, two in short precision and two 
in long precision,. the last being comput- 
ed as 


SQRT(x) = ya + (X/ys - Y3)/2 


to minimize the truncation error. 


The maximum relative error in the result 
from this algorithm is 2**(-65.7). 


Effect of an Argument Error: 
The relative error caused in the result 


is approximately half of the relative 
error in the argument. 


Accuracy: 

a ee at a MM p c EE CLE qucm 
| Arguments | Relative Error | 
| | *10**15 | 
eo Lr | EOS TSE 
| Range |Distribution| RMS | Maximum | 
p----------4----------— +-------4--------- 
|x > 10** |Exponential | 0.0276| 0.124 | 
| (-52) | | | | 
A SEE ESKE AAA l ass] 


Error and Exceptional Conditions: 
P:x< 0 
Implementation: 
e Module size: 160 bytes 
e Execution times: 
Approximate execution times in micro- 


seconds for the System/360 models given 
below are obtained from the table: 


GIA MES te el drei i ee ge EN tae ina | 
[Module | 30 | 40 | 50 | 65 | 75 | 
Pci RSS NUR um t DIE 
|IHESQL | 8282 | 1733 | 376 | 97.7| 57.2 | 
leccata dul. Le Lira Lu <= = dica J 


EXP (short floating-point real) 


Module Name: IHEEXS 


Entry Point: IHEEXSO 


F 


M 


E 


A 


r 


F 


E 


unction: To calculate e to the power x. 
ethod: 


If x < -180.2, a zero result is returned 
immediately. 


Otherwise, EXP(x) is calculated as: 
2** (x*LOG2 (e) ) 


The calculation is performed as follows: 


x*LOG2(e) = f + N, 
where N = 4h + g, h is an integer such 
that g = 0, 1, 2 or 3, and O < f < 1.0. 


Then, by suktracting 0.5, this is reduced 
to the range -0.5 < £ < 0.5. Next, 2**f 
is calculated as: 


la + E/(c + x*x)+x)/(a + b/(c + x*x)-x). 


This is multiplied by 2**0.5 and then 
shifted in the appropriate direction to 
give the effect of multiplication by 
2**g. Finally, the exponent of the 
result is oktained from h. 


ffect of Argument Error: 


The relative error caused in the result 
is approximately equal to the aksolute 
error in the argument, i.e., to the 


argument relative error multiplied by x. 
Thus for large values of x, even the 
round-off error of the argument causes a 
substantial relative error in the answer. 


ccuracy: 
SoS Se e TET messen ae 
Arguments | Relative Error | 
| *10**6 | 
ARRE nay ras isla e I 
Range |Distribution| RMS | Maximum | 
A A ptc 
-1 < x < 1| Uniform | 0.132 | 0.490 | 
wtih cur Gu qi S FS UE QE Sa SES + ee eee +--------- 
Full Range| Uniform | 1.29 | 2.61 | 
IL PPP E ME BS ee ieee AN A LS xj 


rror and Exceptional Conditions: 


I : OVERFLOW if x > 174.673 
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Module Name: 
Entry Point: 


Function: To 


Implementation: 


* Module size: 232 bytes 


e Execution times: 
Approximate execution times in micro- 


seconds for the System/360 models given 
below are obtained from the table: 


ERES RP EP RSS Wm SATIRO tee ee | 
[Module | 30 | 40 | 50 | 65 | 75 | 

Name | | | | | | 
------- ł------ł------}------ł}-----ł}------4 
|IHEEXS | 3847 | 1172 | 356 | 90.0| 58.0 | 
tose be d'auccsta Ll == ==. do Lose J 


EXP (long floating-point real) 


IHEEXL 
IHEEXLO 


calculate e to the power x. 


Method: 


If x < -180.2183, return zero result. 
x/LOG (2) 


Otherwise let y = 
= 4*a - b - c/16 - d 


where a, b and c are integers, 0 < b < 3, 
0 < c < 15 and 0 < d< 1/16. 
Then EXP(x) = 2**y 
= 16**a*2**(-b)*2**(-c/16) 
*2**(-D. 
Compute 2**(-d) by using the Chebyshev 


interpolation polynomial of degree 6 over 
the range 0 < d < 1/16, with maximum 
relative error 2**(-57). 


If c > 0, multiply 2**(-d) by 2**(-c/16). 
The constants 2**(-c/16), 1s c < 15, are 
included in the subroutine. 

If b > 0, halve the result b times. 


Finally, multiply by 16**a by adding a to 
the characteristic of the result. 


Effect of an Argument Error: 


The relative error caused in the result 
is approximately equal to the absolute 
error in the argument, i.e., to the 


argument relative error multiplied by x. 
Thus for large values of x, even the 
round-off error of the argument causes a 
substantial relative error in the answer. 
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Accuracy: Then m a*(1+y)/(1-y) and ABS(y) < 1/3. 


p7-7--------------------- T----------------- 1 Now x = 2**(ü*p-b)*(1*y)/(1-y). 
| Arguments | Relative Error | 
| | *10**15 | Therefore 
}----------7------------ }------- T--------- 1 
| Range Iristribution| RMS | Maximum | LOG(x) = (ü*p-b)*LOG(2) + 
[---------- I------------ +------- I--------- 1 LOG((1+y)/(1-y)). 
[71 <x < 1| Uniform | 0.0674| 0.216 | 
=--------- 4$- 4-------4---------] LOG((1+y)/(1-y)) is computed by using the 
|Full fange| Uniform | 0.867 | 2.30 | Chebyshev interpolation polynomial of 
-—-----2-2---A---2--------- de nico doin degree 4 in y**2 for the range 
O < y**2 < 1/9, with maximum relative 
Error and Exceptional Conditions: error 2**-27.8. LOG2(x) or LOG10(x) are 
calculated by multiplying LOG(x) by 
I : OVERFLOW if x > 174.673 LOG2(e) or LOG10(e) respectively. 


Implementation: 
Effect of Argument Error: 
e Module size: 448 bytes 
The aksolute error caused in the result 


e Execution times: is approximately equal to the relative 
error in the argument. Thus if the 
Approximate execution times in micro- argument is close to 1, even the round- 
seconds for the System/360 models given off error of the argument causes a 
below are obtained from the table: substantial relative error in the answer, 
since the function value there is very 
Sa A SR LEI LI E small. 
[Module | 30 | 40 | 50 | 65 | 75 | 
| Name | | | | | | 
p------- j------ 4------ i------ i----- i------ j Accuracy: 
| IHEEXL | 12131| 2901 | 616 | 343 | 194 | 
m isc A E ES E SSO PA o a ----------------- 
| Arguments | Relative Error | 
| *10**6 
}----------7------------ ł------- T--------- 1 
| Range [Distribution | RMS | Maximum | 
LOG, LOG2, ICG10 (short floating-point LL-2--------l------------ i-2------— dino J 
real) ; IHELNSE 
MER Ne n aM ea a Traino ad SS 
| Excluding | | | | 
Module Name: IHELNS J0.5 < x |Exponential | 0.0320| 0.577 | 
p 2.0 | | | | 
Entry Points: .  —» —— LYVi---------- di --------- i------- 1-------~-- 4 
tiene) 
Mathematical PL/I Entry |  p---------- JOR T------- T--------- 4 
function name point |Excludin g | | | | 
10.5 < x |Exponential | 0.342 | 0.754 | 
Log x to the base e LOG (x) IHELNSE |< 2. 0 | | | | 
Log x to the kase 2 LOG2 (x) IHELNS2 t---~------- d-_---------- doom. L-_____--- 4 
Log x to the kase 10 LOG10 (x) IHELNSD ana 
mm queen he ee oe 
Function: To calculate log x. | Excluding | | | | 
|0.5 < x |Exponential | 0.170 | 0.963 | 
Method: |< 2.0 | | | | 
—— Liz ia l-24 bonos 4 


Let x = m*16**p where 1/16 < m < 1 and p 
is an integer. 


Two constants, a (= base point) and E (= 
-LOG2(a)), are defined as follows: 


1/16 m « 1/8, a = 1/16, b= 4; 
1/8 <m < 1/2, a = 1/4, b= 2; 
172 sm <i, a= 1, b= 0. 


Let y = (m-a)/(m*a) 
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e e a —À—— — a e 3 


PASAS AS AAA T 
Arguments | Absolute Error | 
| *10%**6 | 
os ==> NAE ATE 
Range |Distribution| RMS | Maximum | 
¡AAA A eeu mA Loco esos itouch J 
IHELNSE 
(SATA AA THAT q === 
10.5 < x | Uniform | 0.0960] 0.394 | 
|< 2.0 | | | | 
"PI Pesos we AE edes moni 
IHELNS2 
putem m ree eee quee Trees 
[0.5 < x | Uniform | 0.177 | 0.842 | 
p< 2.0 | | | | 
AAA ba cole ol Ss dul cw È UN SOIN eee EIOS 4 
IHELNSD 
A Ce eh rr ia 
10.5 € x | Uniform | 0.0526] 0.164 | 
|< 2.0 | | | | 
Los sois. Ace tus logos LR SE J 


Implerentation: 
e Mcdule size 
e Execution t 


Approximate 


: 256 bytes 


imes: 


execution times in micro- 
seconds for the System/360 models given 
below are obtained from the appropriate 
entry point in the table: 


m ti ae qo UU NEL eA ds ver. UN ru S T A 
[Entry | 30 | 40 | 50 | 65 | 75 | 
| Point| | | | | | 
------- ł------ł------ł}------ł4-----}------4 
| IHELNSE| 4669 | 1238 | 385 | 173 | 95.7} 
-------}------ ł------ł}------ł}-----}------ 
|IHELNS2| 5041 | 1342 | 417 | 180 | 101.3] 
------- +------4------1------4-----4------ 
|IBELNSD| 5054 | 1366 | 417 | 180 | 101.3] 
A Do o o o e o le e o o SE o e ls o o o o e e o o e o o do o o o e o a el 


LOG, LOG2, IOG10_ (long floating-point real) 


Module Name; IHELNL 
Entry Points: 
Mathematical PL/I Entry 
Function name point 
Log x to the base e LOG(x) IHELNLE 
Log x to the base 2 LOG2 (x) IHELNL2 
Log x to the kase 10 LOG10 (x) IHELNLD 


Function: To calculate log x. 
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Method: 


Let x = 16**p*2**(-q)*m where p is the 
exponent, q is an integer such that 
0 < q < 3, and 1/2 <m < 1. 


Two constants, a (= base point) and b (= 
-LOG2(a)), are defined as follows: 


1/2 < m < 1/SQRT(2: a = 172, b= 1 
1/SQRT(2) <m < 1: a = 1, b= 0 
Let y = (m - a)/(m + a). 

Then m = a*(1 + y)/(1 - y) and ABS (y) 
< 0.1716. 
Now x = 2**(U*p - q- b)*(1 + y)/(1 - y) 
Therefore 

LOG (x) = (4*p - q - b)*LOG(2) + 


LOG((1 + y/(1 - y)). 


LOG((1 + y)/(1 - y)) is computed by using 
the Chebyshev interpolation polynomial of 
degree 7 in y**2 for the range 
O < y**2 < 0.02944, with maximum relative 
error 2**(-59.6). 


LOG2(x) or LOG10(x) is calculated by 
multiplying the result by LOG2(e) or 
LOG10(e) respectively. 


Effect of an Argument Error: 


The absolute error caused in the result 
is approximately equal to the relative 
error in the argument. Thus if the 
argument is close to 1, even the round- 
off error of the argument causes a 
substantial relative error in the answer, 
since the function value there is very 
small. 
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Accuracy: 


PEPE SE A SE Tet SEERE AS 1 
| Argurents | Relative Error | 
| | *10**15 | 
}~---------7------------ +------- t---------] 
| Range [Distribution] RMS | Maximum | 
llli Lions Loss LS J 
IHEINIE 
pium rp Tope e mE he ee ee ed 
| Excluding | | | | 
10.5 < x |Exponentiai | 0.0530| 0.329 | 
|< 2.0 | | | | 
AAA dis ÉT eee dida J 
IHELNI2 
Se A eg ts ree eee od 
|Excludin g | | | | 
|0.5 < x |Exponential | 0.443 | 2.60 | 
|< 2.0 | | | | 
S Mites A A umen Sid 
CHEERED 
pore n qune quee foo ET 1 
| Excluding | | | | 
(0.5 « x |Exponential | 0.155 | 0.402 | 
|< 2.0 | | | | 
bull cese | src ESS decli be ce es 3 
fee memet us Teor eo ee 1 
| Arguments | Absolute Error | 
| *10**15 | 
Fr" Sa ee ds presse 4 
| Range [Distribution] RMS | Maximum | 
Lali FT AA ds oceania J 
IHELNIE 
bcd qoe he a er] 
10.5 < x | Uniform | 0.192 | 0.507 | 
|< 2.0 | | | | 
Up Beep ee er ee aec. douce È IRPI E EEE J 
IHELNI2 
possono ec ane tec a on [= ong Se 1 
[0.5 « x | Uniform | 0.245 | 0.466 | 
|< 2.0 | | | | 
Loc2ca2encse Lo ascos Los 4 
IHELNLD 
masses poor en non nn nnn pn 
[0.5 < x | Uniform | 0.0318| 0. 0625] 
|< 2.0 | | | | 
(CLS Liso aos das Loca 3 
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Implementation: 


e Module size: 360 bytes 


e Execution times: 


Approximate execution times in micro- 
seconds for the System/360 models given 
below are obtained from the appropriate 
entry point in the table: 


B Tessa a Wears ce ere, FE AERE AU quoc TOGO 1 
[Entry | 30 | 40 | 50 |65 | 75 | 

Point | | | | | | 
~------}------ }------}------}-----}------ 
|IHELNLE| 16216| 3926 | 788 | 178 | 98.5| 
~------4------ }------4-~----}-----}------4 
|IHELNL2| 17315| 4196 | 834 | 190 | 107 | 
p------- +------ $------ ł------ $-----}------ { 
|IHELNED| 17284| 4192 | 828 | 188 | 105 | 
| a ee Oe nN ee eee SES: Lor mms Al-- das J 


SIN, SIND, COS, COSD (short floating-point 
real) 


Module Name: IHESNS 
Entry Points: 
Mathematical PL/I Entry 
function name point 
Sin(x radians) SIN(x) IHESNSS 
Sin(x degrees) SIND (x) IHESNSZ 
Cos(x radians) COS (x) IHESNSC 
Cos(x degrees) COSD (x) IHESNSK 
Function: To calculate sin x or cos x. 
Method: 
Let k = pi/4 
Evaluate p = ABS(x)*(1/k) if x is in 
radians 
Or p = ABS(x)*(1/45) if x is in 
degrees, 


using long-precision multiplication to 
safeguard accuracy. 


Separate p into integer part q and frac- 


tional part r, i.e, p=q + r where 
Q rex. 
Define gi = q if SIN or SIND is required 
and x 2 0; 
qa = q + 2 if COS or COSD is 
required; 
qa = q +4 if SIN or SIND is 
required and x < 0. 
Then for all values of x each case has 
been reduced to the computation of 


SIN(k* (q4+r)) = SIN(t) say, where t 2 0. 


Let ga = MOD(q,,8). 

If ga = 0, SIN(t) = SIN(k*r) 

If q2 =1, SIN(t) = COS(k*(1-r)) 
If q2 = 2, SIN(t) = COS(k*r) 

If qa = 3, SIN(t) = SIN(k*(1-r)) 
If q2 = 4, SIN(t) = -SIN(k*r) 

If gz = 5, SIN(t) = -COS(k*(1-r)) 
If 2 = 6, SIN(t) = -COS(k*r) 

If q2 = 7, SIN(t) = -SIN(k*(1-r)). 


Thus it is necessary to compute only 
SIN(k*r,) or COS(k*r,) where r, = or 
1 - r and 0 < r, < 1. 


This is performed by using the Chebyshev 
interpolation polynomials of degree 3 in 
r4**2, with maximum relative error of 
2**-28.1 in the sine polynomial and 
2**-20.6 in the cosine polynomial. 


Effect of an Argument Error: 


The absolute error of the answer is 
approximately equal to the absolute error 
in the argument. Hence, the larger the 
argument, the larger its absolute error 
and the larger the aksolute error of the 
result. Since the function diminishes 
periodically for koth sine and cosine, no 
consistent  ccntrol of the relative error 
can be maintained outside the range -fi/2 
to pi/2 radians (or -90 to +90 degrees). 


Accuracy: 

N ae IET qoom er e 1 
| Arguments | Absolute Error | 
| | — *10**6 | 
p----------1---7-7-------—- jeeceeceqooeec ucc 1 
| Range [Distribution] RMS | Maximum | 
AAA A eme d AS Loans J 
IHESNSS 

AARÓN ION: toral jar md 
{ABS(x) < | Uniform | 0.0557] 0.126 | 
|pi/ 2 | | | | 
---------- +----------- }------- 4--------- 1 
[pi/2< | | | | 
| ABS (x) | Uniform | 0.0553] 0.148 | 
|< 10 | | | | 
---------- ł------------ł-------ł---------4 
|10 < | | | | 
|ABS(x) | Uniform | 0.0560| 0.143 | 
|< 100 | | | | 
Listas Lili. isis i ate J 
IHESNSC 

DELA Tore ala cm Uc tcm) 
¡O < x < pil Uniform | 0.0553| 0.149 | 
~--------- 444 
l-10 sx f | | | 
|< 0, | Uniform | 0.0571| 0.154 | 
lei < x | | | ME 
|< 10 | | | | 
---------- ł------------ł-------4---------4 
|10 < en | | | 
| ABS (x) | Uniform | 0.0553| 0.142 | 
|< 100 | l | | 
rc ica location AA J 


"UE e eee E PROS SE mmm 
| Arguments | Relative Error | 
| *10**6 | 
}----------7------------ +------- 1--------- 1 
| Range [Distribution | RMS | Maximum | 
——— — — AAA RES] 
IHESNSS 
mo merci queccmeec ce qe COEUR 
[ABS GO < | Uniform | 0.198 | 1.40 | 
[Pid | | | | 
IUE tossico ia 


Error and Exceptional Conditions: 


P : IHESNSS, IHESNSC: 
ABS(x) 2 2**18*pi 


IHESNSZ, IHESNSK: 
ABS(x) 2 2**18*180 
Implementation: 
e Module size: 320 bytes 
e Execution times: 
Approximate execution times in micro- 
seconds for the System/360 models given 


below are obtained from the appropriate 
entry point in the table: 


RARE E DERE WT SA ee Wy ety ye gee E Gee te em  d 
[ABS (x) | 30 | 40 | 50 | 65 | 75 | 
——— — A E A A AA A | 
sinti 

D TORT Sete EN TR RE T 
|< pi/4 | 4091 | 1120 | 333 | 85.0] 50.6 | 
------- — —M—M—— M 
|2 pi/4 | 4386 | 1190 | 362 | 92.5| 53.8 | 
——— cs dle———————-————-cl-ezclccccclocclcJ 
IHESNSC 

CIARA A A! che ee FOÀ A 
|< pi/4 | 4078 | 1115 | 329 | 83.6] 49.9 | 
------- ł------ł}------ł}------}-----}------ 
|2 pi/4 | 4373 | 1184 | 357 | 91.0| 53.1 | 
Lo Lom mic dicci dom eee me dicci Lui == J 
IHESNSZ 

AAA a E ad ee E ee 
|< 45 | 4026 | 1132 | 338 | 86.3| 51.6 | 
E------- ł------ ł------ ł------ 4-----4------ 
|2 45 | 4421 | 1202 | 366 | 93.7| 54.8 | 
| A. lc LE a Las. Lens ili inno inn J 
IHESNSK 

NI DC IM > ee A eI, fa pe ee di 
|< 45 | 3693 | 1127 | 334 | 84.8| 51.0 | 
p------- j------ j------ 4------ 4-----4------ 4 
|2 45 | 4408 | 1196 | 362 | 92.3| 54.1 | 
bs decano Lo 1 sòsà du Los J 
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SIN, SIND, COS, COSD (long floating-point 


real) 


Module Name: IHESNL 

Entry Points: 
Mathematical PL/I Entry 

function name point 

Sin(x radians) SIN(x) IHESNLS 
Sin(x degrees) SIND(x) IHESNLZ 
Cos(x radians) COS (x) IHESNIC 
Cos (x degrees) COSD (x) IHESNLK 


Function: To calculate sin x or cos x. 
Method: 
Let y = ABS(x)/(pi/4) for x in radians, 
or y = ABS(x)/45 for x in degrees, 
and y = q + r, q integral, 0 < r < 1. 
Take qa = q for SIN or SIND with positive 
or zero argument, 
dqa = q + 2 for COS or COSD, 
qı =q + 4 for SIN or SIND with 
negative argument, 
and qz = MOD(q,,8). 


SIN(ABS(x) + pi/2) 


Since COS(x) = 
) = SIN(ABS(x) + pi), 


and SIN(-x 


it is only necessary to find 


SIN(pi/4* (q2 + r)), for 0 < qa £ 7. 
Therefore compute: 

SIN(pi/4*r), if q2 = 0 or 4, 

COS(pi/4*(1 - r)), if qz = 1 or 5, 

COS(pi/4*r), . if q2 = 2 or 6, 

SIN(pi/4*(1 - r)) if qa = 3 or 7. 


SIN(pi/U*r4)/r;, where ra is r or 
(1 - r), is computed by using the Cheb- 
yshev interpolation polynomial of degree 
6 in ri**2, in the range 0 < r,**2 s 1, 
with maximum relative error 2**(-58). 


COS (pi/4*r,) is 
Chebyshev interpolation 


computed by using the 
polynomial of 


degree 7 in r4,**2, in the range 
0 < ra,**2 < 1, with maximum relative 
error 2**(-64.3). 

Finally, if qa 2 4 a negative sign is 


given to the result. 
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Effect of an Argument Error: 


The absolute error of the answer is 
approximately equal to the absolute error 
in the argument. Hence, the larger the 
argument, the larger its absolute error 
and the larger the absolute error of the 
result. Since the function diminishes 
periodically for both sine and cosine, no 
consistent control of the relative error 
can be maintained outside the range -pi/2 
to pi/2 radians (or -90 to +90 degrees). 


Accuracy: 

IHESNLS 
[PODIO Se ee ASS tC oe SEE 1 
| Arguments | Relative Error | 

| *10**15 | 
CHANGE scies seta bMS E LOHN MON 
| Range [Distribution| RMS | Maximum | 
da ls dq {------------ ---~---+4--------- 
|-pi/2 < x | Uniform | 0.0542]| 0.381 | 
[< pi/2 | | | 
paraa A A er — esd 

IHESNLC 
RS EE QUEEN IE ae ee 1 
| Arguments | Absolute Error | 

| *10**15 | 
AER He ep eR TE 
| Range |Distribution| RMS | Maximum | 
E $------------4---------------- 
|-pi/2 < x | Uniform | 0.0604] 0.168 | 
|< pi/2 | | | | 
mc T Mins eo a A BE RE EDER] 


Error and Exceptional Conditions: 


IHESNLC: 
2 2**50*pi 


P : IHESNLS, 


ABS (x) 


IHESNLZ, IHESNLK: 
ABS(x) 2 2**50*180 


Implementation: 
e Module size: 416 bytes 

e Execution times: 
Approximate execution times in micro- 
seconds for the System/360 models given 


below are obtained from the appropriate 
entry point in the table: 


| ARE E arenes oe ae Tq ara ‘tot oe qr ERR MERE 1 
|Entry | 30 | 40 | 50 | 65 | 75 | 
| Point| | | | | | 
-------}------ dente 42-41 
|IHESNLS| 13654| 3290 | 661 | 155 | 85.3 | 
-------4------ $------ j------ j----------- 
|IHESNLC| 13641] 3288 | 654 | 153 | 84.3 | 
~------}------}------ A 
|IHESNLZ| 13689| 3302 | 665 | 157 | 86.3 | 
~------}------}------ ł------ł-----4------ 
| IHESNLK| 13676| 3300 | 659 | 155 | 85.3 | 
base eda an an a en an ae a nam om om om ow A---—-—- dc din im me J 


TAN, TAND (short floating-point real) Of pi/2, an argument error will produce a 
large relative error in the result. 


Module Name: IHETNS 


Accuracy: 
Entry Points: 
SSS ens See Em me Teo oo mem eT 
Mathematical PL/I Entry | Arguments | Relative Error | 
function nare point | | *10**6 

MM ļ---------- q------------ ł------- T--------- 4 
Tan(x radians) TAN (x) IHETNSR | Range |Distribution| RMS | Maximum | 
Tan(x degrees) TAND (x) IHETNSD Lion Sad nn 1------- dorm J 

IHETNSR 
ANA quoc mre quee qus menm 1 
Function: To calculate tan x. |ABS(x) < | Uniform | 0.319 | 1.92 | 
|pi/u | | | | 
Method: p---------- +------------ +------- I--------- 1 
lpi/4 < | | | | 
Evaluate p = (4/pi)*ABS(x) if x is in | ABS (x) | Uniform | 0.465 | 1.24 | 
radians, |< 1.5 | | | | 
or p = (1/45)*ABS(x) if x is in ----------- I------------ +------- I--------- 1 
degrees, |pi/4 < | | | | 
[ABS (x) | Uniform | 3.14 | 170* | 
using  long-precision multiplication to |< pi/2 | | | | 
safeguard accuracy. p---------- Ho e nono +------- I-]-------- d 
|pi/2 < | | | | 
Let q andr ke respectively the integral | ABS (x) | Uniform | 1.25 | 70.6* | 
and fractional rarts of p. |s 10 | | | | 
---------- ł------------ł-------ł}---------] 
If q is even, put s = 1; 110 < | | | | 
if q is odd, put s = i-r. | ABS (x) | Uniform | 3.57 | 205* | 
|= 100 | | | | 
Let qı = MOD(q,4). Then dr l1------------— d___-_-- L________- 4 
If qa = 0, TAN(ABS(x)) = TAN(pi*s/H) *These maximum errors are those encoun- 
If Ga = 1, TAN(ABS(x)) = COT(pi*s/4) tered in a sample of 5000 points; each 
If Ga = 2, TAN(ABS(x)) = -COT(pi*s/4) figure depends very much on the particu- 
If qı = 3, TAN(ABS(x)) = -TAN(pi*s/4) lar points encountered near the singular- 


ities of the function. 
Compute TAN(pi*s/4) and COT(pi*s/ü) as 
the ratio of two polynomials: 
Error and Exceptional Conditions: 
TAN (pi*s/4) 
COT (pi*s/4) 


S*P(s**2)/Q(s**2) 
Q(s**2)Y (s*P(s**2)) P: IHETNSR: ABS(x) 2 2**18*pi 
IHETNSD: ABS(x) 2 2**18*180 


where P(s**2) 


212.58037 - 12.559912 

*S**2 I: IHETNSR: OVERFLOW 
270.665736 - 71.645273 IHETNSD: OVERFLOW 
*g**2 + s*xu, 


Q(s**2) 


Finally, if x < 0, put 
TAN(x) = -TAN(ABS(x)). 
Effect of an Argument Error: 


The -absolute error of the answer is 
approximately equal to the absolute error 
of the argument multiplied by 
(1 + TAN(x)**2). Hence if x is near an 
odd multiple of pi/2, an argument error 
will ¡produce a large absolute error in 
the answer. 


The relative error in the result is 
approximately equal to twice the absolute 
error in the argument divided by 
SIN(2*x). Hence, if x is near a multiple 
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Implementation: 


e Mcdule size: 


280 bytes 


* Execution times: 


Approximate execution times 


in micro- 


seconds for the System/360 models given 
below are obtained from the appropriate 
entry point in the takle: 


PH 


P—_+4 


IE—-4-—4 


TAN, TAND (long floating-roint real) 


Module Name: 


Entry Points: 


Function: 


Meth 


I 


Mathematical 


fu 


nction 


HETNL 


Tan(x radians) 
Tan(x degrees) 


od: 


Evaluate 


or 


Le 


E 
E 


t q 


PL/I 
name 


TAN( 
TAND 


TO calculate tan x. 


T 
| 65 | 75 
FETI CA 4 
qo equ ee 
| 85.8] 51.0 | 
4-----1------ 
| 95.1] 55.0 | 
ins Lati 4 
quoe qoem 1 
| 87.0| 52.1 | 
1----- Fon 1 
| 96.3| 56.0 | 
y nro N st 4 
Entry 
point 
x) IHETNIR 
(x) IHETNID 


(4/pi)*ABS(x) if x is in radians 
(1/45) *ABS(x) if x is in degrees. 


and r ke respectively the integral 
and fractional parts of p. 


If q is even, 


put s = r; 


If q is odd, put s= 1 - r. 
Let qa = MOD(q,4). Then 
If ga = 0, TAN(ABS(x)) = 
If qa = 1, TAN(ABS(x)) = 
If qa = 2, TAN(ABS(x)) = - 
If qı = 3, TAN(ABS(x)) = - 
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TAN (pi*¥s/4) 
COT (pits/4) 
COT (pi*¥s/4) 
TAN (pi*s/4) 


Compute TAN(pi*s/4) and COT(pi*s/4) 
the ratio of two polynomials: 


as 


TAN(pi*s/4) = s*P(s**2)/Q(s**2)) 
COT (pits/4) = Q(s**2)/(s*P(s**2)) 
where P(s**2) is of degree 3 and Q(s**2) 


is of degree 4 in s**2, and maximum 
relative error is 3.4*10**-19, 


Finally, if x < 0, TAN(x) -TAN(ABS(x)). 


Effect of an Argument Error: 


The absolute error in the result is 
approximately equal to the absolute error 
in the argument multiplied by 
(1*TAN(x)**2). Hence, if x is near an 
odd multiple of pi/2, an argument error 
will produce a large absolute error in 
the result. 


The relative error in the result is 
approximately equal to twice the absolute 
error in the argument divided by 
SIN(2*x). Hence, if x is near a multiple 
of pi/2, an argument error will produce a 
large relative error in the result. 


Accuracy: 
IHETNLR 
[ESSE Gc EE Qa CO qu TAN 1 
| Arguments | Relative Error | 
| | *10**15 | 
}---------- Til +e 4 
| Range |Distribution| RMS | Maximum | 
---------- }------------}-------}--------- 
|ABS(x) < | Uniform | 0.091| 0.530| 
|pi/4 | | | | 
Sam 4----------- ł------- ł--------- { 
|pi/4 < | | | | 
|ABS(x) | Uniform | 0.437| 2.31 | 
|< 1.5 | | | | 
---------- ł-------- ----}------- ł--------- 1 
|eis < do | | | 
|ABS(x) | Uniform | 7.75 | 416* | 
|< pi/2 | | | | 
---------- ł------------ ł------- ł--------- 4 
|pi/2 < | | | | 
|ABS(x) | Uniform | 18.3 | 1140* | 
|= 10 | | | | 
~--------- det} 
[10 < | | | | 
| ABS (x) | Uniform | 271 [13400* | 
|= 100 | | | | 
==> cua AAA P ER I zs 32s J 
*These maximum errors are those encoun- 
tered in a sample of 5000 points; each 
figure depends very much on the particu- 


lar points encountered near the singular- 
ities of the function. 


Error and Exceptional Conditions: Method: 


P: IHETNLR: ABS(x) 2 2**50*pi 1. ATAN(y,x) 
IHETNID: ABS(x) 2 2**50*180 
If x = 0 or ABS(y/x) 2 2**24, the answer 


I : IHETNLR: OVERFLOW SIGN(y)*pi/2 is returned except for the 
IHETNID: OVERFLOW error case x = y = 0. Otherwise 
Implementation: ATAN(y,x) = ATAN(y/x) if x > 0 
or ATAN(y,x) = ATAN(y/x) + SIGN(y)*pi 
e Mcdule size: 352 if x < 0. 
e Execution times: Hence the computation is now reduced to 


the single argument case. 
Approximate execution times in micro- 


seconds for the System/360 models given 2. ATAN(x) 
below are obtained from the appropriate 
entry point in the table: The general case may be reduced to the 
range 0 < x < 1 since 
(= que perm quem Tuer qoc 1 
JABS(x) | 30 | 40 | 50 | 65 | 75 | ATAN(-x) - -ATAN(x), and 
L.------— dim d--_--- Lomo. doo d------ 4 ATAN(1/ABS(x)) = pi/2 - ATAN(ABS(x)). 
IHETNIR 
------- T-----*1------1-----<1-----1------1 A further reduction to the range ABS(x) < 
[€ pi/4 | 15440| 3622 | 687 | 154 | 87.2 | TAN (pi/12) is made by using 
------- ł------4------ł------ł-----4------ 
|2 pi/4 | 16130| 3817 | 747 | 169 | 93.9 | ATAN(x) = pi/6 + ATAN((SQRT(3)*x - 1)/ 
LL-2-----L4------ 4------ dono dro l1------ J (x + SQRT(3))). 
IHETNID 
------- T------1------1------1-----1------1 Care is taken to avoid the loss of 
|< 45 | 15475] 3634 | 691 | 155 | 88.2 | Significant digits in computing 
p------- ł------ ł------ 4------ f----- ł------ 1 
|> 45 | 16165| 3829 | 751 | 170 | 94.9 | SORT(3)*x - 1. 
User Eccc Ecl ccu: TERES Elo lecce J 
For the basic range ABS(x) s TAN(pi/12), 
use an approximation formula of the form 
ATAN(X), ATAND(X), ATAN (Y,X), ATAND (Y,X) ATAN(x)/x 7 a * b*x**2 * c/(d * x**2) 


(short floating-point real) 
with relative error less than 2**-27.1. 


Module Name: IHEATS 3. ATAND(x) and ATAND(y,x) 
Entry Points: The treatment is as above with the  addi- 
tion of a final conversion of the result 
Mathematical PL/I Entry to degrees. 
function name point 
Arctan x (radians) ATAN(x) IHEATS1 Effect of an Argument Error: 
Arctan(y/x) (radians) ATAN(y,x) IHEATS2 
Arctan x (degrees) ATAND (x) IHEATS3 Let t = x or y/x; then the absolute error 
Arctan(y/x) (degrees) ATAND(y,x) IHEATS4 of the answer approximates to the abso- 
lute error in t divided by (1 + t**2). 
Function: Hence, for small values of t, the two 
errors are approximately the same; howev- 
To calculate arctan x or arctan(y/x). er, as t becomes larger the effect of the 
The result range is: argument error on the answer error dimin- 
ishes. 
Arctan x (radians) t pi/2 
Arctan(y/x) (radians)+ pi 
Arctan x (degrees) + 900 


Arctan(y/x) (degrees)t 1809 
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Accuracy: 


POR MM 
| Arguments 
| 
p----------7------------ 
| Range [Distribution] 
IHEATS1 


presser sm ditte an fan um am DE A O a 


{Full range |Tangents of 
{numbers 
|uniformly 


[Full range|y 


Urata eae ss d erre dadi 


[distributed 


| Eetween 
|-pi/2 and 


|x 


sin and 


cos of 


{numbers 
| uniformly 
[distributed 
| ketween 
|-pi/2 and 


|ri^ 


2 


I ——————-—4 


I ————————4 


NEUE TET 1 
Relative Error | 
*10**6 | 
------- T--------- 
RMS | Maximum | 
Ill cz J 

ED sus aoe ane sunk aam) Cms pem mmm 
0.443 | 0.958 | 
| | 

| | 

| | 

| | 

| | 

| | 
ro losas 
— TESTA 
0.449 | 1.42 | 
| | 

| | 

| | 

| | 

| | 

| | 

| | 

—— do id 


Error and Exceptional Conditions: 


P 4 


IHEATS2, 


Implerentation: 


IHEATS4: x = y = 0 


e Mcdule size: 


408 bytes 


e Execution times: 


Apprcximate 


execution times in micro- 
seconds for the System/360 models given 
below are obtained from the appropriate 
entry point in the takles: 


TAND(15) 


50 


F—-4 


I —-4-———--—-—-—--—4 


SP La ee | 
50 | 65 | 75 | 

ass cas eas ee d.a ades an in on me d 
Vase at toe = tid 3 

79 | 78. ul 48.4 | 
------ ł-----4------]4 
| | | 

47 | 97.8] 58.0 | 

| | | 

NOE NAME ME 
81 |108 | 64.1 | 

| | | 

------ $-----}------4 
13 | 88.7] 54.5 | 

te an a mm dimo e oe om oe ow nw ean ol 


IHEATS2 


ee eS tte: Stee | 


PR — — + — — — e 41 


— am a as — 


PH — 51 


E-—--———--—-——-4--—4 


F—-4-———--—-—-—--—4 


400 P—+_——+-——+-4 


I—-4-—-——-r——-—-4-— 


—4 


PFP +———+———+ 


SF re ae ger 
I 
] 
i 

I —-4- — — +M 


"ig Sacer Pn EL BER 
I 
I 
P—+-——+-——+—4 


gere rre rara greg 
i 
I 
| 

PR + — — — + — — — 4+ — 4 


ATAN(X), ATAND (X), ATAN (Y,X), ATAND (Y,X) 
(long floating-point real) 


Module Name: 


Entry Points: 


I 


HEATL 


Mathematical 


function 


Arctan x (radians) 


Arctan(y/x) 


(radians) 


Arctan x (degrees) 


Arctan(y/x) 


(degrees) 


PL/ 


nam 


ATAN 
ATAN 
ATAN 
ATAN 


I 
e 


(x) 
(y,x) 
D(x) 
D(Y,X) 


Entry 
point 


IHEATL1 
IHEATL2 
IHEATL3 
IHEATL4 


Function: 
To calculate arctan x or arctan(y/x). 
The result range is: 
Arctan x (radians) t pi/2 
Arctan(y/x) (radians)t pi 
Arctan x (degrees) i 909 


Arctan(y/x) (degrees) +180° 


Method: 


1. ATAN(y,x) 
If x = 0 or ABS(y/x) 2 2**56, the answer 
SIGN(y)*pi/2 is returned except for the 


error case x y = 0. Otherwise 


ATAN(y,x) 
or ATAN(y,x) 


ATAN(y/x) if x > 0 
ATAN(y/x) + SIGN(y)*pi 
if x < 0. 


Hence the computation is now reduced to 


the single argument case. 
2. ATAN(x) 


The general case may be reduced to the 
range 0 s x < 1 since 


ATAN(-x) 
ATAN(1/ABS(x)) 


- ATAN(x), and 
pi/2 - ATAN(ABS(x)). 


A further reduction to the range ABS(x) 
< TAN(pi/12) is made by using 


ATAN(x) = pi/6 + ATAN((SQRT(3)*x - 1)/ 
(x + SORT(3))) 
Care is taken to avoid the loss of 


Significant digits in computing 
SQRT(3)*x - 1 


For the basic range ABS(x) < TAN(pi/12), 
use a continued fraction of the form 


ATAN(x)/x = 1 + a4,*x*x/(b, + x*x + az2/ 
(ba + xxx + a3/ (ba + X*X + 


ay/(Dy + x*x)))) 
with relative error less than 2**(-57.9). 


3. ATAND(x) and ATAND(y,x) 


The treatment is as above with the addi- 
tion of a final conversion of the 
to degrees. l 


result 


Charter 3. 


Effect of an Argument Error: 


Let t = x or y/x; then the absolute error 
of the answer approximates to the abso- 
lute error in t divided by (1 + t**2). 
Hence, for small values of t, the two 
errors are approximately the same; howev- 
er, aS t becomes larger the effect of the 
argument error on the answer error dimin- 
ishes. 


Accuracy: 


Relative Error | 
*10**15 | 


Error and Exceptional Conditions: 


P: IHEATL2, IHEATL4: x = y = 0 


Implementation: 


e Module size: 544 bytes 


e Execution times: 


Approximate execution times in micro- 
seconds for the System/360 models given 
below are obtained from the appropriate 
entry point in the table: 


mme qToT—eeeeq-e—————qee————-————T[45--7--73 
[Entry | 30 | 40 | 50 | 65 | 75 | 
| Point| | | | | | 
p-------i------ +------4------4-----4------ 1 
|IHEATL1| 20472| 4389 | 826 | 181 | 100 | 
porron poo ł------ ł----- ł------ 
|IHEATL2| 23523| 5042 | 967 | 217 | 123 | 
Tec ł------4-----4------ 
|IHEATL3| 21575| 4667 | 871 | 190 | 105 | 
—------4------4------ 4------4-----+----- 
¡AREA 200021 5305 | 1008 | 226 | 128 | 
ee a bene ci duc eho ee eho oes 
SINH, COSH (short floating-point real) 
Module Name: IHESHS 

Entry Points: 

Mathematical PL/I Entry 
Hyperbolic sin x SINH (x) IHESHSS 
Hyperbolic cos x COSH (x) IHESHSC 
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Fu 


Me 


Ef 


nction: 


To calculate 
bolic cos x. 


hyperbolic sin x or hyper- 


thod: 
For IHESHSS, if ABS(x) < 1, use a poly- 
nomial approximation of the seventh 
degree. 
Otherwise, 
SINH(x) = EXP(x)/2 - 0.5/EXP(x), 
COSH(x) = EXP(x)/2 + 0.5/EXP(x). 
These two versions of EXP (x) /2 + 


0.5/EXP(x) are preferable to the equival- 
ent versions of (EXP(x) - 1/EXP(x))/2 
because, in floating-point, 0.5 has three 
more significant bits than 1.0. 


fect of Argument Error: 


The relative error caused in the result 
is approximately as follows: 


SINH: The absolute error in 
divided try TANH(x), 
order of the absolute 
argument for large 
relative error in the 
small x. 


the argument 

i.e., Of the 
error in the 
x, or of the 
argument for 


The absolute error in the argument 
multiplied by TANH(x), i.e., of the 
order of the absolute error in the 
argument. 


COSH: 


Thus, for large values of x, even the 
round-off error of the argument causes a 
substantial relative error in the answer. 


Accuracy: 
FOSS mmm id 1 
| Arguments | Relative Error | 
| | *10**6 | 
I[---------- q------------ 4------- T--------- 1 
| Range |Distribution| RMS | Maximum | 
LESS A A Ea le ar J 
IHESHSS 
fru Sor TSS To > TERRASSA que—m——m 1 
|o < | i | | | 
|ABS(x) < 1| Uniform | 0.200 | 0.932 | 
---------- ł------------ł-------4---------] 
|1 < | i | | | 
|ABS(x) < 2| Uniform | 0.221 | 0.950 | 
Licciana Lo ciclo decicion | sm J 
IHESHSC 
f= === queen TIT == ToT 1 
JABS(x) < 1| Uniform | 0.367 | 0.908 | 
}---------- }------------ ł------- ł--------- 1 
|1 < | | | | 
|ABS(x) « 2|] Uniform | 0.192 | 0.700 | 
LS ee A A A S oe di J 


Error and Exceptional Conditions: 


H : OVERFLOW 


(IHEEXS) . 


Implementation: 


e Module size: 


in real EXP routine 


216 bytes 


e Execution times: 


Approximate execution times 


in micro- 


seconds for the System/360 models given 
below are obtained from the appropriate 
entry point in the table: 


| IHESHSS | 
| ABS (x)<1| 2544 
| ABS (x)>1/5647 


cino cl — — — 


P—d4-——+-——4d 


SINH, COSH (long floating-point real) 


Module Name: I 
Entry Points: 
Mathematical 


function 


Hyperbolic sin x 
Hyperbolic cos x 


Function: 
To 
bolic cos x. 

Method: 


For IHESHSS, 


calculate 


i 


sasso FOSS da nane 

40 | 50 | 65 | 75 | 

| | | | 

a (eo DNE RE 

692 | 228 | 56.3] 33.0 | 

1693 | 526 |144 | 91.8 | 

------ $------}-----}------4 

1648 | 509 |139 | 88.4 | 

PN y EEE AAN 
HESHL 

PL/I Entry 

name point 

SINH (x) IHESHLS 

COSH (x) IHESHLC 


hyperbolic sin x or hyper- 


f ABS(x) < 0.3465736, com- 


pute SINH(x)/x using a polynomial approx- 


imation 


of degree 5 in x**2, with rela- 


tive error less than 2**-61.9 


Otherwise, 


COSH (x) 
SINH(x) 


compute s = 


EXP(ABS(x)); then 


(s + 1/s)/2 
SIGN(x) *(s - 1/s)/2 


Effect of an Argument Error: 


The relative error caused in the result 
is approximately as follows: 


SINH: The absolute error in the argument 
divided by TANH(x), i.e., of the 
order of the aksolute error in the 
argument for large x, or of the 
relative error in the argument for 
small x. 


COSH: The aksolute error in the argument 
argument multiplied by TANH (x), 
i.e., Of the order of the absolute 
error in the argument. 


Thus, for large values of x, even the 
round-off error of the argument causes a 
substantial relative error in the answer. 


Accuracy: 
poete SE SES E Ce Tone em 1 
| Arguments | Relative Error | 
| | *10**15 | 
ES ł------- q---------4 
| Range ]Distribution| RMS | Maximum | 
espe MD NV A Lutti ny a eae ee a J 
IHESHLS 
d RE RUE eee emper eec 
| ABS (x) < | Uniform | 0.0530| 0.217 | 
[0.34657 | | | | 
E---------- Y 1------- 1--------- 1 
[0.34567 « | | | 
xod oid < 5| Uniform | 0.0870] 0.359 | 
ccc ccm A O IS SS 22g 
etus 
ENT IRE PNE q----7------Q(-------q--------- 
[ABS (x) € 5| Uniform | 9.123 | 0.429 | 
> d.eec-nseeeseedueeeuseedeuercuesuuj 


Error and Exceptional Conditions: 


H : OVERFLOW in real EXP routine 
(IBEEXI). 


Implementation: 
e Module size: 264 bytes 


e Execution times: 


Approximate execution times in micro- 


seconds for  System/360 model 


S 


given 


below are obtained from the appropriate 


entry point in the table: 


T T T 

[ABS (x) | 30 | 40 | 50 | 65 | 75 | 
i nia mune mE lhe O AA en eo ecm PA ee oe ae ed 
Mec 

T E sew e We MERLO ONG le OA ee een dS 
|< 0.347] 9024] 2279 | 450 | 101 | 59.0] 
ļ------- ł------ ł------ ł------ $-----$------ { 
[0.347 <| | | | | | 
|ABS(x) | 18634] 4338 | 938 | 215 | 125 | 
|< 174. et | | | | | 
== 23 ==> Benen m UE m Al1-----L1-.-----1-----l1-2.-----—J 
sane 

DIM Db DISEASES Ae ars) SOW muU IUNGERE SE 
is 174. I 18493] 4300 | 924 | 211 | 123 | 
es e o o o o co lo o o o e o Loca PO OS 


TANH_(short floating-point real) 


Module Name: IHETHS 


Entry Point: IHETHSO 


Function: To calculate hyperbolic tan x. 


Method : 
1. ABS(x) s 2**-12 
Return x as result. 


2. 2**-12 < ABS(x) < 0.54931 


Use a transformed continued fraction of 


the form: 


TANH(x)/x = 1-((x**2 + a)/(x**2 + b + 


C/x**2)) 


with relative error less than 2**-27. 


3. 0.54931 < x < 9.011 


Use TANH(x) = 1 - 2/(EXP(2*X) + 1). 


4. xz 9.011 
Return result 1. 
5. x < -0.54931 


TANH (x) = -TANH(-x). 
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Effect of an Argument Error: 


The relative error caused in the result 
is approximately twice the absolute error 
in the argument divided by SINH(2*x). 
Thus for small values of x it is of the 
order of the relative error in the argu- 
ment, and as x increases the effect of 
the argument error is diminished. 


Accuracy: 
precem mo eee E mem 
| Arguments | Relative Error | 
| | = *10**6 | 
pesesasserae REESE ERR ec puse a 
| Range |Distribution| RMS | Maximum | 
-----~---- fł------------4}-------4--------- 
|-0.5 < x | Uniform | 0.174 | 0.867 | 
|< 0.5 | | | | 
---------- d------------4---------------- 
|-9 «x < 9| Uniform | 0.0720| 0.782 | 
Pecos sees D-———— d---—ee- La J 
Implerentation: 


e Module size: 200 bytes 
e Execution times: 
Approximate execution times in micro- 


Seconds for the System/360 models given 
below are obtained from the table: 


prom Te quoe fo M queue 1 
| | 30 | 40 | 50 | 65 | 75 | 
So E fo=--- ł----- punto +------ 1 
|ABS(x) < | 791 | 263 [102 | 28.7 | 21.7 | 
|2**-12 | | | | | | 
[eem cue test t----- peste ssa quee 1 
|2**-12 < ]3033 | 785 |231 | 64.1 | 43.9 | 
{ABS (x) <0.5| | | | | | 
om Fo pp p<] 
[0.5 <= x 15934 |1805 |562 |152 1117 | 
|< 9 | | | | | | 
---------- ¢-----}-----}--~-+------4------4 
|x 29 |1095 | 363 [139 | 40.5 | 35.2 | 
ARI lecca dirai dese cim mim d ee dme 4 
TANH (long floating-point real) 
Module Name: IBETHL 
Entry Point: IHETHLO 
Function: To calculate hyperbolic tan x. 
Method: 

1. ABS(x) < 0.54931 

Compute TANH (x) /x using a rational 

approximation, with relative error less 


than 2**-64.5 
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2. 0.54931 < X < 20.101 


1 - 2/(EXP(2*x) + 1). 


TANH (x) 
3. x 2 20.101 
Return result 1. 
4. x < -0.54931 

TANH(x) = -TANH(-x) 

Effect of an Argument Error: 

The relative error caused in the result 
is approximately twice the absolute error 
in the argument divided by SINH(2*x). 
Thus for small values of x it is of the 
order of the relative error in the argu- 


pent, and as x increases the effect of 
the argument error is diminished. 


Accuracy: 
fae n mum To cc mecs 1 
| Arguments | Relative Error | 
| | *10**15 | 
pe LER VEA SUN AQ pene ee 
| Range |Distribution| RMS | Maximum | 
---------- ł------------}-------}--------- 
|ABS(x) < | Uniform | 0.0440] 0.211 | 
[0.54931 | | | | 
---------- fo----------- +} 
[0.54931 < | | | 
{ABS(x) < 5| Uniform | 0.0250] 0.199 | 
Lacio E AA Lon diLzssssssà J 
Implementation: 
e Module size: 280 bytes 


e Execution times: 


Approximate execution times in micro- 
seconds for the System/360 models given 
below are obtained from the table: 


FESTER alias Mee Wet m Wu Secor 2 FISSE, ea ay 4 
JABS(x) | 30 | 40 | 50 | 65 | 75 | 
------- ł------}------ł------ł-----}------ 
|< 0.549] 12745] 3030 | 564 |123 | 67.9| 
------- j------4------d------4-----4------4 
[0.549 s| | | | | | 
JABS(x) | 16400| 3918 | 878 {205 | 119 | 
|< 20.1 | | | | | | 
------- dont 224244 
|2 20.1 | 1239| 372 | 135 | 39.3] 25.5] 
bores Lia do a mm mmm m do. dici d 


ATANH (short floating-point real) 


Module Name: IHEHTS 


Entry Point: IHEHTSO 


Function: To calculate hyperbolic arctan x. 
Method: 
1. ABS(x) s 0.2 
Use a rational approximation of the form: 
ATANH(x) = x + x ** 3/7 (a + b*x**2) 
2. 0.2 < ABS(X) < 1 
ATANH(x) = -SIGN(x)*0.5*LOG((0.5 - 
ABS(x/2))/(0.5 * ABS(x/2))) 
Effect of an Argurent Error: 


The absolute error caused in the result 


is approximately equal to the absolute 
error in the argument divided by 
(1 - x**2). Thus as x approaches 


*1 or -1, relative error increases rapid- 
ly. Near x = 0, the relative error in 
the result is of the order of that in the 
argument. : 


Accuracy: 
pue ccu cu Mieres qoem oa ETE: 1 
| Arguments | Relative Error | 
| | (*#10**6) | 
=== 42 t MERVUENPA CHINE NES D RRR | Ea 
| Range  |Distribution| RMS | Maximum | 
---------- NS IO 
|-0.8 <x | Uniform | 0.440 | 1.32 | 
|= 0.8 | | | | 
—-------- A | 
{-0.9 <x | Uniform | 0.389 | 1.14 | 
|< 0.9 | | | | 
Lunes cms ada sesos Locas J 
Error and Exceptional Conditions: 
P: ABS(x) 2 1 
Implementation: 
e Module size: 192 bytes 
e Execution times: 
Approximate execution times in micro- 


seconds: for the System/360 models given 
below are oktained from the table: 


p Em pepe Nem qum qoom COSTARE 1 
|ABS(x) | 30 | 40 | 50 | 65 | 75 | 
------- from fra pp p< 
[s 0.2 | 2520 | 667 | 208 | 52.0| 31.3 | 
aan) SER 1------ 1------ t----- 1------ j 
10.2 < | | | | | | 
|ABS(x) | 7091 | 1829| 606 |163 | 94.8 | 
|< 1 | | | | | | 
[EFE SE KE dia Eee: doce pc ila J 
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ATANH (long floating-point real) 


Module Name: IHEHTL 


Entry Point: IHEHTLO 


Function: To calculate hyperbolic arctan x. 


Method: 
1. ABS(x) < 0.25 


Use a Chebyshev polynomial of degree 8 in 
x**2 to compute ATANH(x)/x. 


2. 0.25 « ABS(x) « 1 
ATANH(x) = -SIGN(x)*0.5*LOG((0.5 - 
ABS (x/2))/(0.5 + ABS(x/2))) 
Effect of an Argument Error: 


The absolute error caused in the result 
is approximately equal to the absolute 


error in the argument divided by 
(1 - x**2). Thus as x approaches +1 or 
-1, relative error increases rapidly. 
Near x = 0, the relative error in the 
result is of the order of that in the 
argument. 
Accuracy: 
A ive enm eer ee me qoum mmm mri 
| Arguments | Relative Error | 
| | *10**15 | 
[rece Torone ERR q ea 
| Range |Distribution| RMS | Maximum | 
---------- SS fr 
|ABS(x) < | Uniform | 0.0650| 0.223 | 
[0.25 | | | | 
---------- ł------------ 4-24] 
|ABS(x) < | Uniform | 0.133 | 0.397 | 
[0.95 | | | | 
lec eed Lc ss rr SS J 


Error and Exceptional Conditions: 


P : ABS(x) 2 1 
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Implementation: 


e Module size: 272 bytes 


e Execution times: 
Approximate execution times in micro- 


seconds for the System/360 models given 
below are obtained from the table: 


fT Sat TH TTF = ross T====== ===" ===" 1 
JABS(x) | 30 | 40 | 50 | 65 | 75 | 
------- ł------ł------ł}------4-----4------] 
[s 0.25 | 12252| 3037 | 562 | 121 | 68 | 
------- łt------ł4------4------ł-----4------4 
10.25 < | | | | | | 
|ABS(x) | 20448| 4900 | 1040 | 242 | 137 | 
|< 1 | | | | | | 
bese dcus Et d auras T PES Lacuna 4 
ERF, ERFC (short floating-point real) 
Module Name: IHEEFS 
Entry Points: 
Mathematical PL/I Entry 
function name point 
Errcr function (x) ERF(x) IHEEFSF 
Complement cf error ERFC (x) IHEEFSC 


function (x) 


Function: 


To calculate the error function of x or 
the complement of this function. 


Method: 
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1. OS x s 1.317 


Compute ERF(x)/x ky using a Chebyshev 
interpolation polynomial of degree 6 in 
x**2, with relative error less than 
2%*-24. 

ERFC(x) = 1 - ERF(x) 

(ERFC(x) > 1/16 in this range). 
2. 1.317 < x S ky where k = 2.04000092 
Compute ERFC(x) by using a Chebyshev 


interpolation polynomial of degree 7 in 
(x-k), with absolute error less than 
1.3 + 2**-30. 


ERF(x) = 1 - ERFC(x). 
(ERFC(x) > 1/256 in this range). 


3. k < x < 13.306 


ERFC (x) *x*EXP(x**2) is computed by using 
a Chebyshev interpolation polynomial of 
degree 6 in x**-2, with relative error 
less than 1.2 * 2**-23. 


If x « 3.9192, ERF(x) = 1 - ERFC(x) 
If x2 3.9192, ERF(x) - 1 
4. x 2 13.306 
Results 1 and 0 are returned for ERF(x) 
and ERFC(x) respectively. 
5. x< 0 
ERF (x) - ERF(-x) 


and ERFC(x) 2 - ERFC(-x). 


Effect of an Argument Error: 


result 
equal to the absolute 
argument multiplied by 


The aksolute error caused in the 
is approximately 
error in the 
EXP(-x**2). 


ERF(x): As the magnitude of the argument 
increases from 1, the effect of an argu- 
ment error diminishes rapidly. For small 
x, the relative error of the result is of 


the order of the relative error of the 
argument. 
ERFC(x): For x > 1, ERFC(x) is approxi- 


mately EXP(-x**2)/(2*x). Thus the rela- 
tive error in the result is approximately 
equal to the relative error in the argu- 
ment multiplied by 2*x**2. For negative, 
or small positive, values of x, the 
relative error in the result is approxi- 
mately equal to the absolute error in the 
argument multiplied by EXP(-x**2). 


Accuracy: 
PA e EO D TES SES EQ E 
Arguments | Relative Error | 
; | *10**6 | 
pee RE exe tee qeu emm 1 
Range |Distribution| RMS | Maximum | 
i--2l-cocces y ———— delia Lasa Sea J 

IHEEFSF 

FEER OO quee TISSER SEE qoc rd 
| ABS (x) | Uniform | 0.139 | 0.934 | 
|s 1.3 | | | | 
---=------ dn 
[1.3 < | i | | | 
| ABS (x) | Uniform | 0.0372| 0.263 | 
Is 2 | | | | 
~--------- ł-----+------ł-------}---------4 
2 < ES | | | 
|ABS (x) | Uniform | 0.0347] 0.0605 | 
[= 3.9 | | | | 
lossas M ESE dcin E J 


IHEEFSC 


ere Ey NUTUS SIS Me LR SIMA UNT 
[-3.8 <x | Uniform | 0.297 | 0.930 | 
MN AE Se 
po <x | Uniform | 0.505 | 3.80 | 
ii ipod 
[1.3 <x | Uniform | 0.314 | 4.08 | 
ss — | 
|2 <x | Uniform | 0.367 | 1.45 | 
CEL en DES 
13.9 <x | Uniform | 9.09 | 17.1 | 
|£ 13.3 | | | i | 
| ES 0 o a Li A has. Lor ms e e co e e e o e dJ 
Implementation: 


e Module size: 376 bytes 
e Execution times: 


Approximate execution times. in  micro- 
seconds for the System/360 models given 
below are obtained from the appropriate 
entry point in the table: 


Coe en ne Wes eee o COLLE Lm 1 
|ABS(x) | 30 | 40 | 50 | 65 | 75 | 
ao E Lens do Lia as de o cn o o Allcl-- d 
IHEEFSF 
TERA a a MS 
|< 1.32 | 4354] 1191 | 392 |100 | 57.2 | 
------- ł------4------4------}-----}------4 
{1.32 < | | | | | | 
[ABS(x) | 4613| 1266 | 418 |110 | 62.0 | 
|< 2.04 | | | | | | 
~------ a ee a os BEEN, 
{2.04 < | | | | | | 
|ABS(x) | 10013| 2843 | 868 |228 |140 | 
|< 3.92 | | | | | | 
------- dint 4- 
[2 3.92 | 1530| 473 | 183 | 50.7| 32.0 | 
Los mar dicon A weeee be Lon dlc J 
IHEEFSC 
DRETAR WESS ST MY ee EA 
|< 1.317| 4412| 1204 | 385 {103 | 58.8 | 
------- deco ee +} 
11.32 < | | | | | 
|ABS(x) | 4582| 1262 | 405 |110 | 63.4 | 
Is 2.04 | | | | | | 
------- ł------ł------ł}------ł}-----4------4 
{2.04 < | ~ | | | | | 
|ABS(x) | 9982| 2839 | 854 |228 |142 | 
|< 3.92 | | | | | | 
------- 4444] 
[3.92 < | | | | | | 
|ABS(x) | 10168| 2897 | 879 |236 |147 | 
|< 13.31] | | | | | 
d 4------ +------ +------ ł-----ł------ 1 
[2 13.31| 1598| 499 | 182 | 54.5] 35.5 | 
lasco dini Al---—— Lor de o hee J 


ERF, ERFC (long floating-point real) 


Module Name: IHEEFL 


Entry Points: 


Mathematical PL/I Entry 
function name point 
Error function (x) ERF (x) IHEEFLF 
Complement of error ^ ERFC(x) IHEEFLC 
function(x) 
Function: 


TO calculate the error function of x or 
the complement of this function. 


Method: 
1. OSXS 1 


Compute ERF(x)/x by using a Chebyshev 
interpolation polynomial of -degree 11 in 
x**2, with relative error less than 
1.07*2**-57. 


ERFC(x) = 1-ERF(x). 
(ERFC(x) > 1/16 in this range). 


2. 1 <x s 2.04000092 


Compute ERFC(x) by using a Chebyshev 
interpolation polynomial of degree 18 in 
(x - 1.999999), with absolute error less 
than 1.5*2**-61. 


ERF(x) = 1-ERFC(x). 
(ERFC(x) » 1/256 in this range). 


3. 2.04000092 < x < 13.306 


ERFC(x) is computed by using a Chebyshev 
interpolation polynomial of degree 20 in 
x**-2 for ERFC (x) *x*EXP(x**2), with rela- 
tive error ranging from 2**-53 at 
2.04000092 to 2**-51 at 13.306. 


If x « 6.092, ERF(x) 
If x 2 6.092, ERF(x) 


1-ERFC(x). 
1 


4. x 2 13.306 


Results 1 and 0 are returned for ERF(x) 
and ERFC(x) respectively. 


5. x <0 


ERF (x) 
and ERFC (x) 


-ERF(-x) 
2-ERFC(-x). 
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Effect of an Argument Error: 


The absolute error caused in the result 
is approximately equal to the absolute 
error in the argument multiplied by 
EXP (-x**2). 


ERF (x): As the magnitude of the argument 
increases from 1, the effect of an argu- 
ment error diminishes rapidly. For small 
x, the relative error of the result is of 
the order of the relative error of the 
argument. 


ERFC(x): For x > 1, ERFC(x) is approxi- 
mately EXP(-x**2)/(2*x). Thus the rela- 
tive error in the result is approximately 
equal to the relative error in the argu- 
ment multiplied by 2*x**2. For negative, 
or small positive, values of x, the 
relative error in the result is approxi- 
Mately equal to the absolute error in the 
argument multiplied by EXP(-x**2). 


Accuracy: 
Fre ee ory oe oe ee ERE 1 
| Arguments | Relative Error | 
| | *10+*15 | 
[--7-—------1------------ i-—---7--1------—- 
| Range [Distribution] RMS |Maximum | 
bee eo La pence use dicci LESS 4 
IHEEFLF 
ea RON EK DENIS NNNM MG 
| ABS GO < | Uniform | 0.0280 | 0.202 | 
11.317 | | | | 
-—------- ł------------}--------4--------]4 
|1.317 < | | | | 
| ABS (x) | Uniform | 0.0107 | 0.0291| 
|< 2.04 | | | | 
---------- dee 
[2.04 < | | | | 
| ABS (x) | Uniform | 0.00803| 0.0170| 
I< 6.092 | | | | 
—À— D c, € TERN, 
ii 
ASSE | cele on MAE 
|-6 < x | Uniform | 0.0684 | 0.188 | 
|< 0 | | | | 
---------- ł------------ł--------4--------4 
¡O < x | Uniform | 0.0762 | 0.352 | 
|< 1.317 | | | | 
--—-------- ł------------ł--------}--------4 
|< 2.04 | | | | 
~--------- E 4224 
}2.04 <x | Uniform | 1.24 | 4.02 | 
|< 4 | | | | 
T-—--7----- O a 
14 < x | Uniform | 1.40 | 5.02 | 
|< 13.3 | | | | 
bisse K EEE T — l-— — J 
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implementation: 


e Modul 


e size: 


744 bytes 


e Execution times: 


Approximate execution times 


in 


micro- 


seconds for the System/360 models given 
below are obtained from the appropriate 


entry 


point in the table: 
-T 
30 | 40 | 50 
ic aus ima ii ce io Alcl--zedllo2-- 
ATRIA ES 
16567| 4154| 805 
ui us qi a Eni ai NET i ma 
24413| 6095] 1175 
| | 
sanno ni azioni 
45574| 1105] 2142 
| | 
a ec I------1- ai Marcio ame QE 
2493] 707| 222 
sadaa Li di 
pa pa Se ee 
16905| 4240| 837 
a DOT ae 
24263] 6065] 1166 
| | 
eee oe as 
45424 | 10974| 2133 
| | 
eee ¡o oia 
45624] 1t040| 2160 
| | 
—— I------14------ 
2408] 693| 219 
"E PER, dius um an m om 


FUNCTIONS WITH COMPLEX ARGUMENTS 


SQRT (short floating-point complex) 


Module Name: i 
Entry Point: I 
Function: 

To 


square root of z, 


HESQW 


HESQWO 


of result < pi/2. 


i.e., 


75 | 

ons a es i ee E J 
PATATA 
103 | 
------ { 
| 

147 | 
| 
pe | 
269 | 
| 

36.3| 
pm J 
HRS DES 1 
106 | 
------ 1 
| 

147 | 
| 
------ { 
. | 
270 | 
| 
------ { 
| 

274 | 
| 

37.5] 


calculate the principal value of the 
-pi/2 < argument 


Method: Implementation: 


Let z = x + yI, and e Module size: 152 bytes 
SCRT(z) = u + vI. 
e Execution times: 
1. x=y=0 


Approximate execution times in micro- 
Then u = v = 0. seconds for the System/360 models given 
below are obtained from the table: 


2. x20 
pre oes T€ pasa moss a i 1 
Then u = SQRT((ABS(x) + ABS(x + y1))/2) {Module | .30 | 40 | 50 | 65 | 75 | 
and v = y/(2*u). | Name | | | | | | 
I------- ł------ ł------ ł------ ł-----ł------ 4 
3. x< 0 |IHESQW | 11130| 3023 | 1006 | 265 | 164 | 
lc dica Ls d wee deem Loa J 
Then u = y/(2*v) 
and v = S(y)*SQRT((ABS(x) + 


ABS(x + yI))/2) 


where S(y) 1 if y 20 


-1 if y < 0 


SQRT (long floating-point complex) 


Effect of an Argument Error: 
Module Name: IHESQZ 
Let Z 
SQRT (z) 


r*EXP(hI), and 
S*EXP(kI). Entry Point: IHESQZO 


Then the relative error in s is approxi- 

mately half the relative error in r, and Function: 

the relative error in k is approximately 

equal to the relative error in h. To calculate the principal value of the 
square root of Z, i.e., -pi/2 < argument 
of result < pi/2. 


Accuracy: 
Method: 
poeseos eem que mcer Exe 1 
| Arguments | Relative Error | Let z = x + yI, and 
| | *10**6 | SQRT(z) = u + vi. 
pm qecemcosecrce t------- Trosososse 
| Range |Distribution| RMS | Maximum | 1. x=y= 0 
AA AAA 4-------4--------- 
|Full range|Exponential | 0.513 | 1.51 | Then u= v= 0. 
| | radially, | | | 
| | uniform | | | 2. x20 
| |round origin| | | 
LL--------- d_________-- d-____-- d_______- J Then u = SQRT((ABS(x) + ABS(x + yI))/2) 
and v = y/(2*u). 
Error and Exceptional Conditions: 3. x < 0 
I : OVERFLOW Then u = y/(2*v) 
and v = S(y)*SQORT((ABS(x) + 
H : OVERFLOW in complex ABS routine ABS (x + y1))/2) 
(IHEARW) 
where S(y) = 1if y2 0 
= -1 if y< 0 


Effect of an Argument Error: 


Let Z 
SORT (z) 


r*EXP(hI), and 
S*EXP (kI). 


Then the relative error in s is approxi- 
mately half the relative error in r, and 
the relative error in k is approximately 
equal to the relative error in h. 
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Accuracy: 


O eee qu Dmm s AO 
| Arguments | Relative Error | 
| | *10*415 | 
panas A Rio AO EAS pe ==ec==i 
| Range |Distribution| RMS | Maximum | 
}----------}------------ ł------- ł--------- 1 
|Full range|Exponential | 0.263 | 1.54 | 
| | xadially, | | | 
| | uniform | | | 
| [round origin| | | 
— AAA A Issa 4 
Errcr and Exceptional Conditions: 

I : OVERFLOW 

H : OVERFLOW in complex ABS routine 

(IHEABZ) 

Implementation 


e Mcdule size: 144 bytes 


e Execution times: 


Approximate execution times in micro- 
seconds for System/360 models given 
below are oktained from the takle: 


TARA esr on ICE idrato son eo | 
[Module | 30 | 40 | 50 | 65 | 75 | 
| Nare | | | | | | 
------- ł------ł------ł------ł-----4------4 
|IHESCZ | 26996| 5957 | 1352 | 341 | 203 | 
Er as Li Al---——- bs donc diaz 4 


EXP (short floating-point complex) 


Module Name: IHEEXW 


Entry 


Point: IHEEXWO 


Function: To calculate e to the power z. 


Method: 


Let 


Then REAL(EXP(z)) 


and 


z = x + yI. 


EXP (x) *COS (y) 
EXP (x) *SIN(y). 


IMAG (EXP (z) ) 


Effect of an Argurent Error: 


Let EXP(x + yI) = s*EXP(kI). 


Then k = y, and the relative error in s 


is 


approximately equal to the absolute 


error in x. 
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Accuracy: 


BEST mecs qo emm 
| Arguments | Relative Error | 
| | *10**6 | 
arras iO TS Tass ===] 
| Range |Distribution| RMS | Maximum | 
~----=---- ł------------ f-------}---------4 
| ABS (x) | Uniform | 1.32 | 3.14 | 
|= 170 | | | | 
[ABS(y) — | | | | 
|£ pi/2 | | | | 
---------- }-~---------- A 
| ABS (x) | Uniform | 1.31 | 3.34 | 
|< 170 | | | | 
|riz2« f| | | | 
|ABS (y) | | | | 
|s 20 | | | | 
[TNR NUR RS dí lr clau Plz Loria J 


Error and Exceptional Conditions: 


O : ABS(y) 2 2**18*pi : error caused in 
real SIN routine (IHESNS) 


H : OVERFLOW in real EXP routine (IHEEXS) 


Implementation: 
e Module size: 136 bytes 
e Execution times: 
Approximate execution times in micro- 


seconds for the System/360 models given 
below are obtained from the table: 


CIARA a e S We m e O 
[Module | 30 | 40 | 50 | 65 | 75 | 
| Name | | | | | | 
------- $------}------}+-----}-----}------4 
|IHEEXW | 14240] 4058 | 1244 | 326 | 200 | 
[ ma d ia m be lic d a d eme J 


EXP (long floating-point complex) 


Module Name: IHEEXZ 


Entry Point: IHEEXZO 
Function: To calculate e to the power z. 


Method: 
Let z= x + yI. 


Then REAL(EXP(z)) 
and IMAG(EXP(z)) 


EXP (x) *COS (y) 
EXP (x) *SIN(y). 


Effect of an Argument Error: 
Let EXP(x + yI) = s*EXP(kI). 
Then k = y, and the relative error in s 


is approximately equal to the aksolute 
error in x. 


Accuracy: 
(Sate Soe er n Tm qoe ala ia 
| Arguments | Relative Error | 
| *10**15 | 
Fs q AA toe 1 
| Range |Distribution| RMS | Maximum | 
---------- O 
|ABS(x) < 1| Uniform | 0.136 | 0.478 | 
| ABS (y) | | | | 
|< pi/2 | | | | 
ļ---------- 4------------ ł------- ł--------- 1 
| ABS (x) | Uniform | 1.28 | 2.29 | 
|« 20 | | | | 
| ABS (y) | | | | 
|< 20 | | | | 
A doncs cocos Lella PL aw 4 


Error and Exceptional Conditions: 


O : ABS(y) 2 2**50*pi : error 
real SIN routine (IHESNI) 


caused in 


H : OVERFLCW in real EXP routine (IHEEXL) 


Implementation: 
e Module size: 136 bytes 
e Execution times: 
Approximate execution times in micro- 


seconds for System/360 models given 
below are obtained from the table: 


IEEE AERA Ve wc xcd ae ede ds ee | 
|Module | 30 | 40 | 50 | 65 | 75 | 
| Nare | | | | | | 
------- ł------ł------ł------ł-----}------4 
|IHEEXZ | 42838] 10560| 2174 | 505 | 287 | 
bose base boss de Al--—- All2l-- J 


LOG (short floating-roint complex) 


Module Name: IHELNW 


Entry Point: IHELNWO 
Function: 
TO calculate the principal value of the 


natural log of Z, i.e., -pi < imaginary 
part of result < pi. 


Chapter 3. 


Method: 


Let LOG(x + yl) = u + VI. 


Then u = LOG(ABS(x + yI)) 
= LOG(SQRT(x**2 + y**2)) 
= LOG(x**2 + y#*2)/2 
and v = ATAN(y,x). 


In computing u, the exponents of x and y 
are modified if necessary to avoid OVER- 
FLOW or UNDERFLOW, with the appropriate 
correction being applied after the logar- 
ithm has been taken. 


Effect of an Argument Error: 


Let z = r¥EXP( hI) and LOG(z) = u + vI. 


Then the absolute error in u is approxi- 
mately equal to the relative error in r. 
For the absolute error in 
v(= h = ATAN(y,x)), see corresponding 


paragraph for module IHEATS. 


Accuracy: 
(Seer te Sea ee TET Te ee SSR 1 
| Arguments | Relative Error | 
l *10**6 | 
eses re == A peo] 
| Range |Distribution| RMS | Maximum | 
[----------}------------ ł------- }--------- 1 
|Full range|Exponential | 0.150 | 2.27 | 
| except | radially, | | | 
[within | uniform | | | 
|10**-6 of |round origin| | | 
|1*0I | | | | 
———À É——— Lin. AS J 
Error and Exceptional Conditions: 
O: x = y = 0, error in real LOG routine 


(IHELNS) 
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Implementation: Accuracy: 


® Module size: 272 bytes p7---------------------- T----------------- 1 

| Arguments | Relative Error | 

2 Execution times: | *10**15 | 

p---------- T------------ +------- T--------- 1 

Approximate execution times in micro- | Range |Distribution| R.M.S.| Maximum | 

seconds for the System/360 models given [---------- I------------ I-]------ I-]-------- 1 

below are obtained from the table: |Full range|Exponential | 0.0558| 1.46 | 

| except | radially, | | | 

r------- q------ q------ T77-7-7-- T----- T------ 1 [within | uniform | | | 

| | 30 | 40 | 50 | 65 | 75 | |10**-6 of |round origin | | | 

p------- ł------ ł------ł}------ ł-----4------ 4 |1+01 | | | | 

| Gi) | 215112] 3414 | 1078 | 308 | 183 | diri d-______-_--- d__-_--- d_______-- J 
------- f------}------}------+-----4------4 
| (ii) | 11688] 3489 | 1104 | 318 | 190 | 

------- T-2-----4------4------4-----4------4 Error and Exceptional Conditions: 

|(iii) | 11800] 3520 | 1117 | 321 | 193 | 

L-L-.-----L-.-----— 1------ do l-2----— dono 4 O: x = y= 0, error caused in log rou- 


tine (IHELNL) 
(i) ABS(x) and ABS(y) < SQRT(8)*16**31 
and either ABS(x) or ABS(y) 2 16**-30 


Implementation: 
(ii) Either ABS(x) or ABS(y) 
2 SQRT(8)*16**31 e Module size: 288 bytes 
(iii) ABS(x) and ABS(y) « 16**-30 e Execution times: 


Approximate execution times in micro- 
Seconds for System/360 models given 
below are obtained from the table: 


LOG (long floating-point complex) 


(eres TO yc TASA Pe ro” 1 
| | 30 | 40 | 50 | 65 | 75 | 
Module Name: IHELNZ I------- +------ +------ +------ tooo f------ 1 
| Ci) | 44101] 10166| 2086 | 480 | 274 | 
Entry Point: IHELNZO I------- +------ I------ I------ $----- I------ d 
| (ii) | 44398| 10316| 2155 | 507 | 290 | 
Function: I------- I------4------14------ +----~ I------ 1 
[(iii) | 44438] 10325] 2156 | 507 | 290 | 
TO calculate the principal value of natu- LL-.----- domo do don i-----— i------ 4 
ral log of z, i.e., -pi < imaginary part 
of result s pi. (i) ABS(x) and ABS(y) < SQRT(8)*16**31 
and either ABS(x) or ABS(y) 2 16**-26 
Method: 
(ii) Either ABS(x) or ABS(y) 
Let LOG(x + yl) = u + VI. 2 SORT(8)*16**31 
Then u LOG(ABS(x + yI)) (iii) ABS(x) and ABS(y) < 16**-26 


LOG(SQRT(x**2 * y**2)) 
LOG(x**2 + y**2)/2 
ATAN(y,x). 


and v 


In computing u, the exponents of x and y SIN, SINH, COS, COSH (short floating-point 
are modified if necessary to avoid  OVER- complex) 


FLOW or UNDERFLOW, with the appropriate 
correction being applied after the logar- 


ithm has keen taken. Module Name: IHESNW 
Effect of an Argument Error: Entry Points: 

Let z = r*EXP(hI) and LOG(z) = u + vi. Mathematical PL/I Entry 
function name point 

Then the absolute error in u is approxi- 

mately equal to the relative error in r. Sin z SIN(z) IHESNWS 

For the absolute error in Hyperbolic sin z SINH(z) IHESNWZ 

v(= h = ATAN(y,x)) see the corresponding Cos z cos (z) IHESNWC 

paragraph for module IHEATL. Hyperbolic cos z COSH (z) IHESNWK 
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Function: Accuracy: 


TO calculate sin z or hyperbolic sin 2, p[T---------------2-2------ 


q7---------------- 
or cos z or hyperkolic cos z. | Arguments | Relative Error | 
| *10**6 | 
——————————p------------ Po pc e o e 
Method: | Range [Distribution] RMS | Maximum | 
lenses AA Losses Lino. J 
Let z= x + yl. IHESNWS 
l toc x te ae pmo quem mpm > TARA 1 
Then REAL(SIN(z)) = SIN(x)*COSH (y) [ABS (x) | Uniform | 0.721 | 2.07 | 
and IMAG(SIN(z)) = COS(x)*SINH(y); [s 10, | | | | 
TABE s 1| | | | 
REAL(COS(z)) = COS(x)*COSH(y)  t---------- 4—--~-------— L1l2-----—— d_______ J 
and IMAG(COS(z)) = -SIN(x)*SINH (y); Tease 
SS ae Po A O 
REAL(SINH(z)) = COS (y)*SINH(xX) [ABS (x) | Uniform | 0.561 | 1.86 | 
and IMAG(SINH(z)) = SIN(y)*COSH(x); [s 10, | | | | 
| ABS Cy? $ 1| | | | 
REAL(COSH(z)) = COS(Cy)*COSH(X) | |  L---------- dono dro mo —— Lai mme J 
and IMAG(COSH(z)) = SIN(y)*SINH(x). sante 
ti en i E Tease sae ini SEES arnt 1 
To avoid making calls to evaluate SINH [ABS (x) | Uniform | 0.546 | 2.00 | 
and COSH separately, and thus frequently js 10, | | | | 
having to evaluate EXP twice for the same | ABS (y) | | | | 
argument, SINH(u) is computed as follows: iS 20 | | | | 
a A AA AAN 
1. u > 0.3465736 SEEK 
PUR uM Te qux cix fee SES dix cuenta A 
SINH(u) = (EXP(u) - 1/EXP(u))/2. [ABS GO | Uniform | 0.558 | 2:35. 4] 
Is 10, | | | | 
2. 0 < u < 0.3465736 |ABS (y) | | | | 
[s 20 | | | | 
SINH(u)/u is approximated by a polynomial t---------- dr non mk da d-_______ 4 
of the form a + a4*u**2 + ag*u**ü (which 
has a relative error of less than 
2**-26.4). Error and Exceptional Conditions: 
3. uso O : IHESNWS, IHESNWC: 
ABS(x) 2 2**18*pi : error caused in 
SINH(u) - -SINH(-u). Then real SIN routine (IHESNS) 
COSH(u) = SINH(ABS(u)) + 1/EXP(ABS(u)). 


IHESNWZ, IHESNWK: 
ABS(y) 2 2**18*pi : error caused in 
real SIN routine (IHESNS) 
Effect of an Argument Error: 
H : OVERFLOW in real EXP routine (IHEEXS) 

Combine the effects on SIN, COS, SINH and 
COSH according to the method of evalua- 
tion described in the above paragraph. 
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Implementation: 


e Module size: 320 bytes 


® Execution times: 


Approximate execution times in micro 
seconds for the System/360 models give 
below are obtained from the arpropriat 
entry point in the tables: 


SIN, COS : ABS(y) > 0.3465736 
SINH,COSH: ABS(x) > 0.3465736 
orrore I Too meses For 
|Entry | 30 | 40 | 50 | 65 | 75 
| Point| | | | | 
~--~-~--}------}------ +------4-----+------ 
|IHESNWS| 15826| 1508 | 648 | 363 | 223 
------- O O 
|IHESNWC| 15898| 1518 | 653 | 366 | 225 
~------}------}------ $------}-----}------ 
|IHESNWZ| 15930| 1520 | 653 | 366 | 226 
-------}------}------ $------+-----++----- 
| IHESNWK| 15900] 1519 | 655 | 367 | 227 
A NS A ES AE PS 
SIN, COS : ABS(y) < 0.3465736 
SINH,COSH: ABS(x) < 0.3465736 
pao === qe quee T^--T—-- a A TASAS 
| IHESNWS| 16896| 1585 | 674 | 381 | 232 
~------}------}------ f------}-----4------ 
|IHESNWC| 16968| 1595 | 679 | 384 | 234 
~~-----}------}------ +------4-----4+------ 
| IHESNWZ| 17000| 1596 | 679 | 384 | 234 
A }------}-----}------ 
|IHESNWK| 16970| 1595 | 681 | 384 | 235 
LEA ET BÆR as NA d uvm Lia Limia Liar 
SIN, SINH, COS, COSH (long floating-point 
complex) 
Module Name: IHESNZ 
Entry Points: 
Mathematical PL/I Entry 
function name point 
Sin z SIN(z) IHESNZ 
Hyperbclic sin z SINH(z) IHESNZ 
Cos Z COS (2) IHESNZ 
Hyrerkolic cos z COSH (z) IHESNZ 


Function: 


To calculate sin z or hyperkclic sin z 
or cos z or hyperkolic cos z. 
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Method: 


Let z = x + yl. 


Then REAL(SIN(z)) SIN (x) *COSH (y) 


and  IMAG(SIN(z)) COS (x) *SINH (y) ; 
REAL(COS(z)) =  COS(x)*COSH(y) 
and IMAG(COS(z)) = -SIN(x)*SINH(y); 
REAL(SINH(z)) = COS (y) *SINH(x) 
and IMAG(SINH(z)) = SIN(y)*COSH(x); 
REAL(COSH(z)) = COS (y) *COSH(x) 
and IMAG(COSH(z)) = SIN(y)*SINH(x). 


To avoid making calls to evaluate SINH 
and COSH separately, and thus frequently 
having to evaluate EXP twice for the same 
argument, SINH(u) is computed as follows: 


1. u > 0.3465736 
SINH(u) = (EXP(u) - 1/EXP(u))/2 

2. 0 < u < 0.3465736 

SINH(u)/u is approximated by a polynomial 


of the fifth degree in u**2 which has a 
relative error of less than 2**-61.8 


3. u<O 
SINH(u) = -SINH(-u). Then 
COSH(u) = SINH(ABS(u)) + 1/EXP(ABS(u)). 


Effect of an Argument Error: 


Combine the effects on SIN,COS,SINH and 
COSH according to the method of evalua- 
tion described in the above paragraph. 


Accuracy: 


| ABS (x) 
|< 10, 
LARS (9 


| ABS (x) 
|< 10, 
| ABS (y) 


| ABS (x) 
|s 10, 
| ABS Cy) 
la 20 


PL ——« 


H——— 


Pd k= 


k—— 


Implementation: 


e Module size: 368 bytes 


Relative Error | 


Error and Exceptional Conditions: 


O : 


IHESNZZ, 
ABS(y) 2 2**50*pi: 


IHESNZS, 
ABS(x) 2 2**50*pi: 


IHESNZC: 


*10**15 | e Execution times: 
pres q---------] 
RMS | Maximum | Approximate execution times in micro- 
ll----—---- J seconds for the System/360 models given 
below are obtained from the appropriate 
wees T---------1 entry point in the tables: 
2.11 | 82.4 | 
| | SIN, COS : ABS(y) > 0.3465736 
| | SINH,COSH: ABS(x) > 0.3465736 
undue e issssssssså ra ida is ta 
jEntry | 30 | 40 | 50 | 65 | 75 | 
------- T-—------ | Point| | | | | | 
0.180 | 2.31 | |[-------1------1------ +------ +----- +------ 4 
| | |IHESNZS| 46584| 11267| 2363 | 552 | 313 | 
| | F------- ł------ }------}------ +-----}------ : 
------- dro md |IHESNZC| 46656| 11294| 2373 | 555 | 315 | 
Sp pon t-----+------ 4 
------- T--------- |IHESNZZ| 46726| 11317| 2378 | 557 | 316 | 
0.389 | 6.33 |  [L-------4------1t4------1------ I-----1------ 4 
| | | IBESNZK | iio 113041 2374 | 556 | 316 | 
| | bebe ed y EA tas J 
ele licia 
------- JOR SIN, COS : ABS(Y) < 0.3465736 
1.11 | 19.4 | SINH, coon: ABS(X) < 0.3465736 
| LS O ES yes pss 1 
| | | IHESNZS| 54173] 13141] 2656 | 612 | 345 | 
| | de ł------ +-----4------ 1 
------- d-______-J | IHESNZC| 54245| 13168| 2666 | 615 | 347 | 
}------- ł------ sper EAE oa as 4 
|IHESNZZ| 54325| 13191| 2671 | 617 | 347 | 
~------}------}------}------ $-----}-+----4 
|EHESRAK| STENT SEI 2667 | 618 | 348 | 
oii AS ee oi cn Seto Sod 


error 


real SIN routine (IHESNL) 


IHESNZK: 


error caused 


real SIN routine (IHESNL) 


caused in 


in 


TAN, TANH (short floating-point complex) 


H : OVERFLOW in real EXP routine (IHEEXI) 
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Module Name: IHETNW 
Entry Points: 
Mathematical PL/I Entry 
function name point 
Tan z TAN(z) IHETNWN 
TANH (z) IHETNWH 


Hyperbolic tan z 


Function: 


TO calculate tan z or hyperbolic tan z. 
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Method: Implementation: 


Let z = x + yl. e Module size: 184 bytes 
Then REAL(TAN(z)) - e Execution times: 
TAN(x)*(1 - TANH(y)**2)/ 
(1 + (TAN(x)*TANH(y))**2), Approximate execution times in micro- 
seconds for the System/360 models given 
IMAG(TAN(z)) = below are obtained from the formulas: 
TANH(y)*(1 + TAN(x)**2)/ 
(1 + (TAN(x)*TANH(y))**2). IHETNWN: a + time for execution of 


IHETHS with argument y 
TANH(z) = - (TAN(ZI))I. 
IHETNWH: b + time for execution of 
IHETHS with argument x 
Effect of an Argument Error: 


IT DEA Tq em T ECT | eras E enero 1 
The aksolute error caused in the result | | 30 | 40 | 50 | 65 | 75 | 
is approximately equal to the absolute p---+------ +------- +------- I------- lI------ 4 
error in the argument divided by | a | 9094 | 2310 | 696 | 186 | 111 | 
ABS (COS (z)**2) for IHETNWN, or divided Ly | t | 9197 | 2454 | 716 | 191 | 115 | 
ABS (COSH(z)**2) for IHETNWH. The rela- London doin o 1------- dl l1.-----— 4 


tive error caused in the result is 
approximately twice the absolute error in 
the argument divided by ABS(SIN(2*2)) for 
IHETNWN, or divided by ABS(SINH(2*z)) for 


IHETNWH. 
TAN, TANH (long floating-point complex) 
Accuracy: 
Se e ooo T-----------------1 Module Name: IHETNZ 
| Arguments | Relative Error | 
| | *10**6 | Entry Points: 
Lsecnecen Toe A =>] : 
| Range |Distribution| RMS | Maximum | Mathematical PL/I Entry 
L-L--------- d_----------- it------- 1--------- 4 function name point 
IHETNWN 
p-7--------- T------------ T------- T--------- 1 Tan z TAN (z) IHETNZN 
|ABS(x) « 1| Uniform | 9.430 | 1.67 | Hyperbolic tan z TANH (z) IHETNZH 
JABS(y) < 9| | | | 
—— rt — Ica ea 4 
IHETNWH Function: 
(FS a : aed {= qeu 1 . 
|ABS(x) « 9| Uniform | 0.430 | 1.45 | To calculate tan z or hyperbolic tan z. 
¡ABS (y) < 1| | | | 
Lasses Ll Les Lo J 
Method: 
Error and Exceptional Conditions: Let z = x + yl. 
I : OVERFLOW Then REAL(TAN(z)) = 
TAN(x)*(1 - TANH(y)**2)/ 
O : ABS(u) 2 2**18*pi, where (1 + (TANGO*TANH(y))**2), 
u = x for IHETNWN, 
u = y for IHETNWH. IMAG(TAN(z)) = 
TANH(y)*(1 + TAN(x)**2)/ 
H : OVERFLOW in real TAN routine (IHETNS) (1 + (TAN(x)*TANH(y))**2). 
TANH(z) = - (TAN(ZI))I. 


66 


Accuracy: 


FE FOSS = ue e mE qoe e cT 1 
| Arguments | Relative Error | 
| | *103*15 | 
}----------7------------ }------- q 1 
| Range [Distribution] RMS | Maximum | 
raus MeL lotti dl cosceuicecectoc 
IHETNZN 

oe re SOT Se dad ia 
| ABS (x) < 1| Uniform | 0.139 | 1.11 | 
[ABS (y) « 9| | | | 
SOLE NN ES oe ee Oe EE AA 
ani 

sE IO: A To aa 
| ABS (x) «9| Uniform | 0.137 | 0.980 | 
[ABS Ty] < 1| | | | 
€ A E mee ad A AS 


Error and EXCEPTIONAL Conditions: 


I OVERFLOW 


O 


ABS(u) 2 2**50*pi, where 
u - x for IHETNZN, 
u = y for IHETNZH. 


ee 


OVERFLOW in real TAN routine (IHETNL) 


Implementation: 
e Module size: 184 bytes 
e Execution times: 
Approximate execution times in micro- 
seconds for the System/360 models given 


below are obtained from the appropriate 
entry point in the takle: 


=esezesmra===e Ue er 
[ABS (y) | 30 | 40 | 50 | 65 | 75 | 
ila ee E RS A RA 
EINEN 

A quen mmm qum mie m Ape m erp o em mg 
|< 0.549] 40843| 9625 | 1844 | 411 | 233 | 
------- ł------ł------4------4-----4------4 
[0.549 <| | | | | | 
|ABS(y) | 44498[10513 |10365 | 493 | 284 | 
|< 20.1 | | | | | | 
]------- ł------ 4------4------ fr fm] 
|2 20.1 | 29337| 6967 | 1415 | 327 | 190 | 
boss se TORNA boos boos sione E OPNS pare J 
[7 —-—-——— quee esere 112237 T^---- TS 1 
|ABS(x) | 30 | 40 | 50 | 65 | 75 | 
towels Ioia fee epee N MEI Lu Ls 4 
IHETNZH 

A ST TEA NAAA 
[< 0. 549 | 411221 9709 | 1871 | 419 | 236 | 
~------}------ $------}------4-----4------4 
[0.549 <| | | | | | 
|ABS(x) | 44777|10597 | 2185 | 501 | 287 | 
|< 20.1 | | | | | | 
------- O E OI E 
12 20.1 | 29616] 7051 | 1442 | 334 | 193 | 
DS Bo Lio sd ah ei 
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ATAN, ATANH (short floating-point complex) 


Module Name: IHEATW 
Entry Points: 
Mathematical PL/I Entry 
function name point 
Arctan z ATAN (Zz) IHEATWN 
Hyperbolic arctan z ATANH(z) IHEATWH 


Function: 
TO calculate arctan z or hyperbolic arc- 
tan Z. 

Method: 
Let z= x + yl. 


Then REAL(ATANH(z)) (ATANH(2*x/ 


(1+x*x+y*y)))/2 


IMAG (ATANH (z) ) (ATAN(2*y, 


(1-x*x-y*y)))/2 


and ATAN(z) = -(ATANH(ZI))I. 


Effect of an Argument Error: 


The absolute error in the result is 
approximately equal to the absolute error 
in the argument divided by (1 + z**2) in 
the case of IHEATWN, or by (1 - z**2) in 
the case of IHEATWH. Thus the effect may 
be considerable in the vicinity of 
z = t 11 (IHEATWN) or + 1 (IHEATWH). 


Accuracy: 
(CSS err a ee quoc es TA 1 
| Arguments | Relative Error | 
| *10**6 | 
}----------7------------ 4------- T--------- 1 
| Range [Distribution | RMS | Maximum | 
DUM MUSS dll aus è REMO Leccosisse J 
IHEATWN 
HEU. MAG NE E 
(Full Fange mxponentiat | 0.216 | 2.88 | 
ps eo DUE EE ee od 
PARANA 
Ss FN au 1 
|Full range| Exponential | 0.208 | 1.18 | 
| NNI NEA MANU c E M A J 


Error and Exceptional Conditions: 


P: IHEATWN: 
IHEATWH: 


1I 


+ 
= +1 


> I+ 


Z 
Zz 
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Implerentation: 


e Module size: 304 bytes 


e Execution times: 


Approximate execution times in micro- 
seconds for the System/360 models given 
kelow are oktained from the appropriate 
entry point in the takle, where 


= ABS (2*u/(1+x*x+y*y) ) 
u = y for IHEATWN 
= x for IHEATWH 


Corn ee ee e qe armel 1 
| | 30 | 40 | 50 | 65 | 75 | 
tasca Lasses Die aera en uic dccus J 
IHEATWN 
poe to ok ips qe eS a RAN 1 
ja < 0.2] 12235] 3306 | 1033 | 279 173 | 

------- ff ff rr] 
(0.2 < al 16056| 4454 | 1408 | 462 | 276 | 
p< 1 | | | | | | 
lola I dirai d aca aia ai y T Tert J 
IHEATWH 
pee qe E que mm que pecore 1 
ja s 0.2] 12100| 3267 | 1017 | 275 | 171 | 
I------- o t------ 1----- ł------ 1 
[0.2 < a| 15921] 4415 | 1392 | 458 | 274 | 
p< 1 l | | | | | 
AA E res doc lxi. Los La J 
ATAN, ATANH (long floating-point complex) 
Module Name: IHEATZ 
Entry Points: 

Mathematical PL/I Entry 
function name point 
Arctan z ATAN(Z) IHEATZN 
Hyperkolic arctan z ATANH (z) IHEATZH 
Function: 
To calculate arctan z or hyperbolic arc- 


tan Z. 
Method: 
Let z= x + yl. 


Then REAL (ATANH(Z) ) (ATANH (2*x/ 


(1*x*x*ty*y)))/2 


IMAG (ATANH (z) ) 


(ATAN(2*y, 
(1-x*x-y*y)))/2 


and ATAN(z) -(ATANH(zI))I. 
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Effect of an Argument Error: 


The aksolute error in the result is 
approximately equal to the absolute error 
in the argument divided by (1 + z**2) in 
the case of IHEATZN, or by (i - z**2) in 
the case of IHEATZH. Thus the effect may 
be considerakle in the vicinity of 
z = t I (IHEATZN) or + 1 (IHEATZH). 


Accuracy: 
PESSESS ec oe qe renere 
| Arguments | Relative Error | 
| *10**15 | 
ge oe S ee ee ee NUT tessere 
| Range [Distribution | RMS | Maximum | 
CSS coa Llantas AAA 4 
IHEATZN 
essre, testo Tc qune dcc M 
|Full range| Exponential | 0.141 | 7.93 | 
POR EEA VIS RI RT t Mm Es 
rimini 
uuum EA $us ci FEE SS Lr eme 
|Full range | Exponential | 0.0826| 1.20 | 
E Bie ees SR ate TUS bea ede ed 


Error and Exceptional Conditions: 


P: IHEATZN: Z 
IHEATZH: z 


+11 
11 


Implementation: 
e Module size: 296 bytes 
e Execution times: 
Approximate execution times in micro- 
seconds for the System/360 models given 


below are obtained from the appropriate 
entry point in the table, where 


a = ABS(2*u/(1*x*x*y*y)) 
u = y for IHEATZN 
= x for IHEATZH 

sogna Teena === CA 1 
| | 30 | 40 | 50 | 65 | 75 | 
— sacara E E docce Macao Lcd J 
IHEATZN 

So. Tot SS ESPE a a 
[aso. 25 | 43477| 9977 | 2006 | 455 | 260 | 
[------- ł------ ł------ ł------ ł-----4------ 4 
[0.25«a | 51673|11840 | 2406 | 576 | 329 | 
p< 1 | | E | | | 
CSS D PI | PN os Lita Lanas 4 
IHEATZH 

ERROR UE quete uc ee ae oe ee 
[aso. 25 | 43293] 9923 | 1987 | 450 | 258 | 
po===--- ł------ ł------ i------ j----- i------ 1 
[0.25«a | 51489/11786 | 2466 | 571 | 327 | 
ES 1 | | | | | | 
da Ll dr ood 


The Library supports the array built-in 
functions SUM, PROD, POLY, ALL and ANY, and 
also provides indexing routines for han- 
dling simple (i.e., consecutively stored) 
and interleaved arrays. 


input Data 


The array function modules are distingu- 
ished from the other Library modules in 
that they all accept array arguments and 
perform their own indexing, whereas the 


other modules require that indexing should 
be handled ky compiled code. Calls to 
conversion routines are included in the 


SUM, PROD and POLY modules with fixed-point 
arguments, .so that these arguments are 
converted to floating-point as they are 
accessed (it should be noted that it is a 


requirement of the language that the 
results from these modules ke in 
floating-point). On the other hand, the 
conversions necessary for the ALL and ANY 


modules (the arguments must be converted to 
bit string arrays) are not part of the 
modules and must ke carried out before the 
modules are invoked. 


Any restrictions on the admissibility of 


arguments are noted under the headings 
'Range' and 'Error and Excepticnal 
Conditions'. 

Range: This states any ranges of arguments 


which a module assumes to have been exclud- 
ed prior to its being called. 


Error and Exceptional  Conditions:. These 
cover conditions which may result from the 


use of a routine; they are listed in four 
categcries: 


P -- Programmed conditions in the module 
concerned. Programmed tests are 
made where this is not too costly 
and, if an invalid argument is 
found, a branch is taken to the 
entry point IHEERRC of the execution 
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error  package(EXEP). This results 
in the printing of an appropriate 
message and in the ERROR condition 
keing raised. 


I -- Interrupt conditions in the module 
concerned. For those routines where 
SIZE and FIXEDOVERFLOW are detected 
by programmed tests or where hard- 
ware interruptions may occur, the 
OVERFLOW, UNDERFLOW, and (when the 
conversion package is called) SIZE 
conditions pass to the ON handler 
(IHEERR) and are treated in the 
normal way. The machine is assumed 
to be enabled for all interruptions 
except significance, which is masked 


off. 

O -- Programmed conditions in modules 
called by the module concerned. 
These occur when invalid arguments 


are detected in the module called. 


H -- As I, but the interrupt conditions 
occur in the modules called by the 
module concerned. 


Effect of Hexadecimal Truncation 


See the corresponding section in the 
introduction to Chapter 3 for guidance to 
the accuracy of SUM, PROD, and POLY. If 
fixed-point arguments are passed to these 
functions, further errors may be introduced 


ky conversions. 


Speed 


The average execution times given are 


based on IBM System/360 Instruction Timing 
Information, Form A22-6825. 


A summary of the Library array modules 
is given in Figures 6 and 7. 
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Simple arrays, and 
interleaved arrays of 
variable length strings 


a eum M ee 1 
Interleaved string | 


arrays with fixed- | 
length elements | 


ns 1 
| Indexers IHEJXS IHEJXI | 
| ALL, ANY IHENL1 IHENL2 I 
PAN A EA AE Lee IS J 
Note:  IHEJXI is used for indexing through interleaved 

arithmetic arrays 

Figure 6. Bit String Array Functions and Array Indexers 
ATE A A | Cer TS ee Mm md 1 
| PL/I | Fixed-point | Floating-point arguments | 
| function | arguments pon a npn a ne ee 4 
| | | Short precision | Long precision | 
| }-------- T----------- ł-------- q----------- ł--------r- a 
| | Simple |Interleaved| Simple |Interleaved| Simple |Interleaved| 
p-------------- ł-------- rn Jo +----------- f-------- +---------—- 1 
| SUM real | IHESSF | IHESMF | IHESSG | IHESMG | IHESSH | IHESMH | 
| complex | IHESSX | IHESMX | IHESSG | IHESMG | IHESSH | IHESMH | 
| | | | | | | | 
| PROD real | IHEPSF | IHEPDF | IHEPSS | IHEPDS | IHEPSL | IHEPDL | 
| complex | IHEPSX | IHEPDX | IHEPSW | IHEPDW | IHEPSZ | IHEPDZ | 
| [-------- Lo porn i----------- T-------- dn 4 
| POLY real | IHEYGF | IHEYGS | IHEYGL | 
| complex | IHEYGX | IHEYGW I IHEYGZ | 
DUREE EDS D ——————— ean ra eee ce a een nee idi Lo DES LEDE J 
Figure 7. Arithmetic Array Functions 
ARRAY INDEXERS Method: 


Indexer for Simple Arrays 


Module Name: 


IHEJXS 


Entry Points: 


Element Entry 
address point 
Bit addresses IHEJXSI 
Byte addresses IHEJXSY 


Function: 
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To find the first and last elements of an 


array. Their addresses are returned, in 
general registers 0 and 1 respectively, 
as bit addresses  (IHEJXSI) or byte 


addresses (IHEJXSY). 


The address of the virtual origin B of 
the array (i.e., the address that would 
correspond to the element A(0,..0)) is 
obtained as a byte address for IHEJXSY, 
or a bit address for IHEJXSI, by refer- 
ring to the first word of the array dope 
vector (ADV). 


Address of first element B + 


Address of last 


B + 


$ 
M¿Lj 
ivi 
T 
n 
element > MU; 
i=1 


where Mj is the multiplier for the ith 
dimension 
Lj is the lower bound for the ith 
‘dimension 
U; is the upper bound for the ith 
dimension, and 


n is the number of dimensions. 


Range: 


0 < number of dimensions < 2**22 


Implementation: 


® Module size: 


e Execution time: 


Let n = 


104 bytes 


the number of dimensions 


Then the approximate execution times in 
onds for the System/360 models 
given below are obtained from the fol- 
lowing formulas: 


wicrosec 


IHEJXSI: 

IHEJXSY: 
¡ cst ee T 
| | 30 | 
E + 
| a | 720 | 
| b | 555 | 
| e | 377 | 
Lod. 4 


c*n 
c*n 


E———-4-—4 


E—-—-—--—4 


Indexer for Interleaved Arrays 


Module Name: 


Entry Points 


IHEJ 


Operation 


XI 


Initialization for bit 


addresse 


S 


Initialization for byte 


addresse 


S 


Elements after the first 


Function: 


PL — — + — 4 


Entry 
point 


IHEJXII 
IHEJXIY 


IHEJXIA 


TO find the next element of an array and 


to return 


its kit or byte 


general register 1. 


address in 


Entry point IHEJXII is used to initialize 
the routine for 
provide the address of the first 
in the array; IHEJXIY does the same for 


byte addre 
used ther 


address i 
the routin 


sses. 
eafter 


e. 


bit addresses 


Entry point 


and to 
element 


IHEJXIA is 


to obtain the addresses 
of subsequent elements of the array; one 
S returned for each entry into 


Method: 


Arrays are stored in row major order. 
Let Li be the lower bound and Uj the 
upper bound of the ith dimension, and n 
the number of dimensions. Starting with 
the element A(I4,L2,......Lpg), the rou- 
tine varies the subscripts through their 
ranges to A(U4,U2,......Up), changing the 
nth subscript most rapidly; in this way 
the elements are referenced in the order 
in which they are stored. 


The routine does not deal with actual 
subscript values but calculates the 
extent E ,(= Uj - Lj + 1) of each dimen- 
sion and uses this as a count that varies 
from Ej to 1 for subscript values Lj to 
Ui. A 'base address' for each dimension 
is maintained and, for the ith dimension, 
is defined as the address of the element 
with ith subscript equal to its lowest 
bound Lj and with all other subscripts at 
their current values. 


Thus initially the base addresses are all 
equal to the address of 
A(L4,La2,..-.«« Ln). Each subsequent ele- 
ment address is generated from the  pre- 
vious one by adding the multiplier Mn 
from the array dope vector (ADV) and 
reducing the subscript count by 1. When 
the count for the ith dimension has been 
reduced from Ej to 1 it is reset to Ei, 
Mi.4, is added to the (i-1)th dimension's 
base address and the count for this 
dimension is decreased by one. 


This new base is the starting point for 
further increments by Mn. When a new 
base address is calculated, the base 
addresses for all higher dimensions 
((i+1), (i+2),......n) is set equal to 
the ith base address. 


Range: 


0 < number of dimensions < 2**22 
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Implerentation: ARRAY FUNCTIONS 
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e Module size: 328 bytes 


ALL (X), ANY (X) 
e Execution time: 


Let n = the number of dimensions of Module Names: 
the array 
Module 
E,; = Uj- L; + 1 Arguments name 
n-1 Simple arrays and interleaved IHENL1 
T = T Ei arrays with variable-length 
i=1 elements 
n-1 n-i Interleaved arrays with fixed- IHENL2 
s=} T Ej length elements 
i=1 j=1 
Entry Points: 
V4 = time to get a VDA (IHESADF) 
Entry 
V2 = tire to free a VDA (IHESAFD) PL/I function point 
Then the approximate execution times in ALL(X), ANY(X), byte-aligned IHENL1A 
ricroseconds for the System/360 models IHENL2A 
given below are obtained from the  fol- 
lowing formulas: ALL(X), any alignment IHENL1L 
IHENL2L 
IHEJXII a + V4 + c*n 
ANY (X), any alignment IHENLI1IN 
IHEJXIY bk + V, + c*n IHELN2N 
IHEJXIA: Function: 
n= 1: d + Va + e*E, The argument X is a bit string array (any 
necessary conversion having been per- 
n> 1: k + V2 + h*S + formed prior to the invocation of these 
T*(f + e*E, + g*n*(n-1)) modules). The result is a scalar bit 
string of length equal to the greatest of 
The total time required to index the current lengths of the elements of X. 
through the complete array is the sum l 
of: ALL(X): the ith bit of the result is 1 if 
the ith bits of all the elements of X 
i) the time for either IHEJXII or exist and are 1; otherwise it is 0. 


IHEJXIY, and 
ANY(X): the ith bit of the result is 1 if 


ii) the time from the appropriate the ith bit of any element of X exists 
IHEJXIA formula and is 1; otherwise it is 0. 
—T-——-—-- queer T--———--— T-———-——- T-————-— 1 Method: 
| 30 | 40 | 50 | 65 | 75 | 
—-+------ +------- +------- +------- I------ 4 For byte-aligned string arrays, AND 
a | 1620 | 526 | 203 | 55.9 | 34.4 | (IHEBSA) and OR (IHEBSO) are used for ALL 
b | 1398 | 472 | 182 | 50.8 | 31.8 | and ANY respectively; for string arrays 
c | 508 | 192 | 75.8 | 20.6 | 14.3 | with any alignment BOOL (IHEBSF) is used. 
d | 122] 51.5 | 8.8 | 2.5 | 1.6 | with appropriate parameter bits. 
e | 786 | 267 | 107 | 28.8 | 19.7 | 
£ | 204 | 66.7 | 25,8 | 7.7 | 5.3 | The elements of the array are passed to 
g | 42 | 14.4 | 4.8 | 1.4 | 1.1 | IHEBSA, IHEBSO or IHEBSF one at a time, 
h | 260 | 95.3 | 36.8 | 11.1 | 8.2 | and the result is developed in the target 
k | 220 | 91.9 | 24.0 | 7.6 | 5.3 | field. For the first call to any of 
do d______- i-2------ dn dr J these logical modules the first element 
of the array serves as both first and 
Note: the fastest overall execution time second source arguments. For subsequent 
will occur when the extent of the calls, the result already developed in 
nth dimension is the largest of the the target field is the first argument 
subscript extents of the array and the next element of the array is the 
dimensions. second argument. 


Range: 


Bit 


strings 


are limited to a maximum of 


32,767 bits. 


Implementation: 


e Module size: 


e Execution time: 


Let R 


Ta 


Ta 
Ta 


Ty 


IHENL1: 280 kytes 

IHENL2: 192 bytes 
number of elements in the 
array 
time required to execute 
appropriate kit string routine 
via © IHEBSAO, IHEBSOO or 
IHEBSFO 


time to index via IHEJXSI 
time to index via IHEJXSY 


sum of times required to index 
via IHEJXII and IHEJXIA 


Then the approximate execution times in 


(X) 


SUM 


Module Names and Entry Points 


Simple Arrays 


Module 
Arguments name 
Fixed, real IHESSF 
Fixed, complex IHES SX 
Short float 
real IHES SG 
complex  IHESSG 
Long float 
real IHESSH 
complex IHESSH 


Entry 
point 


IHESSFO 
IHESSXO 


IHESSGR 
IHESSGC 


IHESSHR 
IHESSHC 


Interleaved Arrays 


Module 
Arguments name 
Fixed, real IHESMF 
Fixed, complex IHESMX 
Short float 
real IHESMG 
complex  IHESMG 
Long float 
real IHESMH 
complex  IHESMH 
Function: 


To produce a scalar with a value which is 
S of the array 


the sum of all the element 
argument. 


Method: 


Entry 
point 


IHESMFO 
IHESMXO 


IHESMGR 
IHESMGC 


IHESMHR 
IHESMHC 


microseconds for the System/360 models 
given below are obtained from the fol- 
lowing formulas: 
IHENL1 
Fixed-length elements: 

a + R*(b + Ty + T2) 
Varying-length elements: 

a + R*(c + T, + T3) 
pee quem re Totem qoem 1 
| 30 | 40 | 50 | 65 | 75 | 
fe forro +------- f------- ł------ 
| 1318 | 470 | 189 | 51.6 | 36.9 | 
| 837 | 258 | 95 | 27.3 | 18.0 | 
| 375 | 134 | 49.8 | 15.2 | 10.7 | 
LESSERE deis icona SEES Lat 4 
IHENL2 

a + Ty + R* (Eb + T4) 
qucm qoem es quere qum m 1 
| 30 | 40 | 50 | 65 | 75 | 
$------}------- 1------- 4------- i------ 
| 1341 | 479 | 188 | 53.5 | 36.6 | 
| 428 | 133 | 54.0 | 15.2 | 11.5 | 
Feeds loci Lo La ds Loidas J 


Error 


The elements of the array are added to 
the current sum in row major order. 


For fixed-point arguments each element is 
converted to floating-point by using the 
PL/I Library conversion package. 


For a complex argument, the summation of 
the real parts is performed before the 
summation of the imaginary parts is begun 
in modules IHESSG and IHESSH, while the 
two sums are developed concurrently in 
other modules. 


and Exceptional Conditions: 
I : OVERFLOW, UNDERFLOW 
H : IHESSF, IHESSX, IHESMF, IHESMX: 
ABS (element of the array) > 
7.2*10**75: SIZE condition caused in 
conversion package 
Chapter 4. Arrays: Array Functions 73 


Impl 


erentation: 
Module Sizes: 75 
Mcdule Bytes 


IHESSF 168 
IHESSX 216 
IHESSG 104 
IHESSH 104 
IHESMF 136 
IHESMX 224 
IHESMG 128 
IHESMH 128 


Execution times: 

Approximate execution times in micro- 
seconds for the System/360 models given 
kelow are obtained from the appropriate 
formula. 

Constants used in these formulas are: 

R = number of elements in the array 


T, = sur of times required to execute 
IHEJXI using IHEJXIY and IHEJXIA 


Ta = time to execute IHEJXS Ly means of 
IHEJXSY 


Tz = time for the appropriate conver- 
sion using IHEDMA 


The binary and decimal source data is 
always fixed-point; target data is 
short or long floating-roint. 


IHESSF 

Scurce Target 

binary short a * Ta * R*(e * T4) 

decimal short b + T2 + R*(e + Ta) 

binary long C + Ta + R*(f + T3) 

decimal long d+ Ta + R*(f + T3) 
A possente Tren for cimus quee 1 
| | 30 | 40 | 50 | 65 | 75 | 
ļ---ł------ fon t------~ ł------- i------ 1 
| a | 1281 | 446 | 177 | 46.6 | 31.6 | 
| E | 1367 | 469 | 186 | 58.7 | 32.2 | 
| c | 1308 | 454 | 179 | 46.6 | 31.4 | 
| a | 1394 | 477 | 188 | 48.6 | 32.0 | 
| e | 119 | 35.5 | 13.4 | 12.4 | 8.3 | 
| £ | 137 | 35.2 | 12.2 | 11.2 | 7.2 | 
A NOA y ee S NS pee dec RECITA I 4 
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I 


S 


HESSX 


ource 


binary 


decimal 


binary 


d 


Pe —————— + — e] 


P———+-5 


ecimal 


-———— + —44 -————— — + —4 


H= 


a + T2 + R* (e*2*T4) 


b + Ta + R*(e*2*T4) 


C + T2 + R*(f*2*Tg) 


d + T2 + R*(£*2*T3) 


Target 
short 
short 
long 
long 
PENE q~- 
40 | 50 
ERR I------- 
492 | 208 
515 | 216 
518 | 315 
540 | 224 
259 | 94.5 
267 | 96.1 
— —— Lo 
a + Ta + R*b 
C + Ta + R*d 
sssaaa T------- 
40 | 50 
—— € 4------- 
301 | 118 
34.9 | 12.4 
347 | 134 
67.8 | 24.8 
sana Lan ca 
a + Ta + R*b 
c + Ta + R*d 
nus T------- 
40 | 50 
PETRUS I------- 
314 | 121 
43.9 | 15.2 
372 | 142 
87.8 | 30.4 
SES d asas 


-———— + —4 


E————--—4 


$ 
| 

+ 
| 
| 
| 
| 
| 
| 

d 


-————+—4 


b--——— + — 4 


OSA CORE 1 
75 | 
------ 1 
40.7 | 
41.3 | 
41.5 | 
42.0 | 
| 

| 

y 


IHESMF 
Target 
short 


long 


E————--—34 


IHESMX 


Source 


binary 
decimal 


binary 


——— tl 


decimal 


E—————-—--—4 
-——— —— — + — 44 


IHESMG 


Real 


Complex 


E————--—4 


PL ———+—«] 


IHESMH 
Real a + T4, + R*b 
a + T4 + R*(e + Ta) Complex C + T4, + R*d 
C + Ti + R*(f + Ta) r---r------ quem pan 22 ===> perso» 1 
| | 30 | 40 | 50 | 65 | 75 | 
------- e  þ---ł------4-------ł}-------}-------}------4 
40 | 50 | 65 | 75 | la | 887 | 298 | 111 | 30.9 | 20.4 | 
------- +-------4-------|------ | b | 366 | 116 | 42.4 | 12.3 | 7.9 | 
363 | 141 | 37.3 | 24.6 | | c | 993 | 337 | 125 | 34.4 | 23.5 | 
381 | 147 | 38.5 | 25.6 | {i a | 514 | 157 | 55.9 | 16.0 | 9.8 | 
139 | 54.4 | 16.2 | 10.9 | ld d-____-- d-____-- don d___--- J 
139 | 53.2 | 15.1 | 9.8 | 
P | mene A er 
PROD (X) 
Target 
short a + Ty + R*(e+2*T3) Module Names and Entry Points: 
short b + T4, + R*(e*2*T4) Simple Arrays 
long C + T4, + R*(£*2*T4) Module Entry 
Arguments name point 
long d + T4, + R*(£+2*T3) 
Fixed, real IHEPSF IHEPSFO 
------- TT7-7-----T-------q------4 Fixed complex IHEPSX IHEPSXO 
40 | 50 | 65 | 75 | Short float 
------- $-------4-----~~-4---~---4 real IHEPSS IHEPSSO 
511 | 229 | 63.8 | 47.8 | complex IHEPSW IHEPSWO 
533 | 238 | 65.9 | 48.4 | Long float 
515 | 227 | 63.7 | 47.4 | real IHEP SL IHEPSLO 
537 | 235 | 65.8 | 48.0 | complex  IHEPSZ IHEPSZO 
232 | 87 | 25.1 | 16.4 | 
240 | 89 | 24.0 | 15.2 | Interleaved Arrays 
———— AE A A 
Module Entry 
Arguments name point 
Fixed, real IHEP DF IHEPDFO 
Fixed, complex IHEP DX IHEPDXO 
a + T4, + R*E Short float 
real IHEPDS IHEPDSO 
c + T4, + R*d complex IHEPDW IHEPDWO 
Long float 
—————— a a e | real IHEPDL IHEPDLO 
40 | 50 | 65 | 75 | complex  IHEPDZ IHEPDZO 
------- +-------4-------1------ 
286 | 107 | 30.3 | 20.4 | Function: 
97.4 | 35.6 | 11.9 | 7.9 | 
320 | 119 | 33.9 | 23.5 | To produce a scalar with a value which is 
129 | 46.3 | 15.6 | 9.8 | the product of all the elements in the 
a nici Lar ad 


array argument. 
Method: 


The elements of the array are used in row 
major order to multiply the current pro- 
duct. 


For fixed-point arguments, each element 
is converted to floating-point by using 
the PL/I Library conversion package. 
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Error and Exceptional Conditions: IHEPSF 


I : OVERFLOW, UNDERFLOW Source Target 
H : IHEPSF, IHEPSX, IHEPDF, IHEPDX: binary short a + Ta + R*(e+Ta) 
ABS (elerent of the array) > 
7.2*10**75: SIZE condition caused in decimal short b + Ta + Rf (e+T2) 
conversion package 
binary long C + Ta + R*(f+T3) 
Implementation: decimal long. A + Ta + R*(f+T3) 
e Module sizes: f77-7T------ T====5== T77-7---—- T------- t------ 1 
| | 30 | 40 | 50 | 65 | 75 | 
Module Bytes ----+------+------- +------- I-]------ I----- -d 
la | 1208 | 416 | 163 | 43.8 | 29.3 | 
IHEPSF 160 | b | 1294 | 439 | 172 | 45.9 | 29.9 | 
IHEPSS 72 | c | 1257 | 434 | 169 | 45.0 | 30.2 | 
IHEPSL 72 | d | 1343 | 456 | 178 | 47.1 | 30.8 | 
IHEPSX 256 | e | 613 | 188 | 62.3 | 16.1 | 10.3 | 
IHEPSW 96 | £ | 1331 | 357 | 74.8 | 18.1 | 11.2 | 
IHEPSZ 96 Ladino l1l------— dm 1------- 1------ 4 
IHEPDF 144 
IHEPDS 88 
IHEPDL 88 
IHEPDX 272 IHEPSS 
IHEPDW 120 
IHEPDZ 120 a+ Tz + Råb 
p ES qoe qo ere ene que ce 1 
e Execution times: | | 30 | 40 | 50 | 65 | 75 | 
[---1------ ł------- ł------- ł------- ł------ 1 
Approximate execution times in micro- | a | 456 | 191 | 82.3 | 24.6 | 17.7 | 
seconds for the System/360 models given | b | 372 | 96.9 | 27.0 | 5.9 | 3.3 | 
below are oktained from the appropriate LL dro 1------- 1------- 1-l----- 4 
formula. 
Constants used in the formulas are: IHEPSL 
R = number of elements in the array a + T2°+ R*b 
T, = sum of times required to execute r77-T------ T------- T------- T------- T------ 1 
IHEJXI using IHEJXIY and IHEJXIA | | 30 | 40 | 50 | 65 | 75 | 
ļ---4------ ł------- ł------- ł------- ł------ 1 
T. = time to execute IHEJXS via | a | -252 | 22.5 | 41.5 | 16.0 | 12.9 | 
IHEJXSY | b | 1112 | 276 | 43.5 | 9.1 | 5.3 | 
ls leo y EA fisher ener dacccLe Mesas J 


Ta = time for the appropriate conver- 
sion using IHEDMA 


The binary and decimal source data is 


always fixed-point; target data is 
short or long floating-foint. 
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IHEPSX IHEPDS 


Source Target a + T, + R*b 
binary short a + Ta + R+ (e+2*T3) rn —yoo--—- TT77-----— TERR o TERR eee 1 
| | 30 | 40 | 50 | 65 | 75 | 
decimal short b + T2 + R*(e*2*T4) ----+------ I------- I------- I------- I------ 4 
| a | 418 | 178 | 80.0 | 23.7 | 17.5 | 
binary long C + Ta + R*(f£+2*T3) | b | 492 | 143 | 43.3 | 11.4 | 7.2 | 
ii hada Lio Lia disease J 
decimal long d + Ta + R*(£*2*T4) 
Sal rese qoc qeu pee qc 1 
| 30 | 40 | 50 | 65 | 75 | IHEPDL 
[---+------ ł------- ł------- ł------- ł------ 4 
a | 1405 | 494 | 192 | 52.1 | 35.5 | a + T, + R*b 
k | 1491 | 517 | 200 | 54.1 | 36.0 | 
c | 1481 | 520 | 199 | 53.2 | 36.2 | f---1T------ T------- T7-7----- T------- T-7----- 4 
a | 1567 | 542 | 208 | 55.2 | 36.8 | | | 30 | 40 | 50 | 65 | 75 | 
e | 1993 | 553 | 3173 | 41.6 | 24.1 | p---+------ I------- I------- I------- I------ 4 
£ | 5043 | 1285 | 241 | 54.1 | 31.0 | | a | -290 | 8.8 | 67.5 | 20.7 | 15.5 | 
---1------ don mo d-~- 4-~----- dono 4 | b | 1264 | 331 | 63.8 | 14.8 | 9.2 | 
A AS Las iuris A es dolce J 
IHEPSW 
IHEPDX 
a + Ta * R*b 
Source Target 
ot qu qucm qoc i —- 1 
| 30 | 40 | 50 | 65 | 75 | binary short a + T4, + R*(e*2*T4) 
F==-Ff-===== eres qe pesa qe 1 
a | 7614 | -69.8 | 7:7 | 8.9 | 10.6 | decimal short b + T, + R*(e*2*T4) 
b | 1508 | 382 | 108 | 24.1 | 12.1 | 
Lol dL_____-- domo d______- i-2----- J binary long C + T, + R*(f*2*T4) 
decimal long d + T4, + R*(£*2*T4) 
IHEPSZ ‘ion moe em E T------- = 1 
| | 30 | 40 | 50 | 65 | 75 | 
a + Ta + R*b E ł------- ł------- ł------- ł------ 1 
| a | 1382 | 497 | 208 | 56.8 | 41.2 | 
poor mm T-7-7----- T------- T-=-----:- r7-7---- 1 | b | 1468 | 520 | 216 | 58.9 | 41.8 | 
| 30 | 40 | 50 | 65 | 75 | | c | 1374 | 501 | 205 | 56.7 | 40.8 | 
----+------ I------- +------- I------- 4------ 1 | d | 1460 | 523 | 214 | 58.8 | 41.4 | 
a |-3622 | -791 | -56.2 | -4.3 | 2.6 | | e | 2047 | 574 | 182 | 44.4 | 26.3 | 
b | 4580 | 1124 | 180 | 37.7 | 20.1 | | £ | 5097 | 1306 | 250 | 56.9 | 33.2 | 
A IO Lara AAA P oie a decas J bs E asco ss A A È AE E EN J 
IHEPDF IHEPDW 
Target a + T4, + R*b 
short a + T, + R*(e + T3) f^7-7T-7^----- T^77---- T^777-7--- PR q---e-- 1 
A | | 30 | 40 | 50 | 65 | 75 | 
long C + T4, + R*(f + Ta) ----+~------ I------- +------- +------- +------ 4 
| a | -814 | -96 | 12.1 | 8.9 | 12.5 | 
r---r------ q------- q------- q------- T------ 1 | b | 1694 | 452 | 132 | 32.0 | 17.7 | 
| 30 | 40 | 50 | 65 | 75 | Lordi L1-------— dro mo d------- dr J 
ļ---4------ ł------- ł------- ł------- ł------4 
a | 1075 | 365 | 141 | 37.8 | 25.0 | 
c | 1124 | 382 | 147 | 39.0 | 25.9 | 
e | 645 | 201 | 69.0 | 18.2 | 12.1 | 
£ | 1363 | 371 | 81.5 | 20.2 | 13.0 | 
A AA dois dz diccllocez Lis J 
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IHEPDZ 


a + Ta + R*b 


ae pos ese ee ype eee oe ome ee ace WT 1 

| 30 | 40 | 50 | 65 | 75 | 
--4------ ł------- ł------- 4------- +------ 4 
a |-3841 | -852 | -61.5 | -4.7 | 4.2 | 
E | 4830 | 1214 | 211 | 46.1 | 25.7 | 
ss tl ee as n'a Lom we duca Acl d. leue J 


POLY (A,X) 


Module Names and Entry Points: 


Module Entry 
Arguments name point 

Fixed, real 

vector X IHEYGF IHEYGFV 

scalar X IHEYGF IHEYGFS 
Fixed, complex 

vector X IHEYGX IHEYGXV 

scalar X IHEYGX IHEYGXS 
Short float, real 

vector X IHEYGS IHEYGSV 

scalar X IHEYGS IHEYGSS 
Short float, complex 

vector X IHEYGW IHEYGWV 

scalar X IHEYGW IHEYGWS 
Long float, real 

vector X IHEYGL IHEYGLV 

scalar X IHEYGL IHEYGLS 
Long float, complex 

vector X IHEYGZ IHEYGZV 

scalar X IHEYGZ IHEYGZS 

Function: 

Vector X: 
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Let the arguments be arrays declared as 


A(m:in) and X(p:q). Then the function 
computed is: 
nem j-1 
Alm) + Y) Alm + p«r X(r + i) 
j=1 i=0 


unless n = m, when result is A(m). 


If qep<n-m- 1, for 


E + i > q, X(f + i) = X(q). 


then, 


Scalar X: 


This may ke interpreted as a special case 


of vector X, that is, a vector with 
element, X(1), which is equal to X. 
the function computed is: 


one 
Then 


nem 
>> Adm + jsx**j 
j=0 


A floating-point result is obtained in 


both cases. 


Method: 


1. Vector X, (q - p2 n-m- 1): 


POLY(A,X) is evaluated by 
plication and addition, 


nested  multi- 
i.e., 


(... (A(N)*X(k) + A(n-1))*X(k-1) + 


where k = p* n- m- 1. 


2. Vector X, (q - p« n-m- 1): 

In the expression above, the terms in X 
with subscript ranging from k down to q 
are all made equal to X(q). The evalua- 
tion is treated as for scalar X until 
sufficient terms in X have been made 
equal to X(q), when the computation  con- 
tinues as in (1.). 


3. Scalar X: 


Terms in X with subscript ranging from k 
to p are equal to X. 


For fixed-point arguments each element is 
converted to floating-point, by using the 
PL/I Library conversion package. 


Error and Exceptional Conditions: 


I : OVERFLOW, UNDERFLOW 


H : IHEYGF, IHEYGX: 

ABS (element of the array) > 
7.2*10**75: SIZE condition caused in 
conversion package 


Implementation: 


e Module sizes: 


Module Bytes 
IHEYGF 432 
IHEYGS 240 
IHEYGL 240 
IHEYGX 688 
IHEYGW 280 
IHEYGZ 280 


e Execution times: 


Let the arguments be declared as A(m:n) 
and X(p:q), or X, and T be the time for 
one conversion using the arithmetic 
conversion director  IHEDMA. Then the 
approximate execution times in micro- 
seconds for the System/360 models given 
are obtained from the appropriate for- 
mula.  'Short' or 'long' refers to the 
floating-point result. 


MAU. OQ HO Oo try 


IHEYGF 


Scalar X: 
short a + 2*T + (n-m)*(b+T) 
long C + 2*T + (n-m)*(d+T) 


Vector X, (q-p2n-m- 1): 


short e + T + (n-m)*(£+2*T) 
long g + T + (n-m)*(h+2*T) 
Vector X, (q- p « n - m - 1): 
short i + 2*T + (n-m)*(b+T) + 
(q-p*t1)*(j*T) 
long k + 2*T + (n-m)*(d+T) + 
(q-p*1)* (147) 
piene quem Tror to T 
| 30 | 40 | 50 | 65 | 
ł------ ł------- ł------- ł------- + 
| 2297 | 834 | 338 | 99.0 | 
| 904 | 387 | 100 | 26.0 | 
| 2408 | 870 |. 351 | 102 | 
| 1706 | 484 | 124 | 30.4 | 
| 2629 | 910 | 370 | 97.1 | 
| 1480 | 459 | 155 | 41.2 | 
| 2740 | 946 | 383 | 99.4 | 
| 3044 | 844 | 200 | 50.0 | 
| 3370 | 1197 | 491 | 140 | 
| 258 | 91 | 35 | 9.5 | 
| 3581 | 1233 | 504 | 143 | 
| 280 | 101 | 39 | 10.7 | 
dazi doccecac Ll doses L 
IHEYGX 
Scalar X: 
short a + 4*T + (n-m)*(b+2*T) 
long C + 4*T + (n-m)*(d+2*T) 
Vector X, (q- p2n-m- 1): 
short e + 2*T + (n-m)*(f+4*T) 
long g + 2*T + (n-m)*(h+4*T) 
Vector X, (q - p«n-m- 1): 
short i + 4*T + (n-m)*(b*2*T) 
(q-p*t1)*(j42*T) 
long k + 4*T + (n-m)*(d+2*T) 


(q-p*t1)*(142*T) 


T T quem ai y ia 1 
| 30 | 40 | 50 | 65 | 75 | 
--i------ +------- 4------- ł------- ł------ 4 
a | 3245 | 1174 | 478 | 142 | 95.4 | 
b | 2345 | 664 | 220 | 58.3 | 33.5 | 
c | 3471 | 1221 | 496 | 145 | 97.6 | 
d | 5519 | 1433 | 301 | 73.2 | 42.6 | 
e | 3187 | 1114 | 447 | 127 | 83.5 | 
£ | 4368 | 1226 | 399 | 103 | 58.1 | 
g | 3533 | 1161 | 465 | 129 | 85.7 | 
h {10636 | 2746 | 556 | 133 | 76.4 | 
i | 4087 | 1459 | 592 | 171 1114 | 
j | 545 | 187 | 72.8 | 21.1 | 12.7 | 
k | 4243 | 1506 | 610 | 176 |116 | 
1 | 567 | 196 | 76.8 | 22.3 | 13.6 | 
TRE AA bolle di Liria Lar J 
IHEYGS, IHEYGL, IHEYGW, IHEYGZ 
Scalar X: a + (n-m)*b 
Vector X: 
(q-p2n-m-1): c + (n-m) *d 
(q-p<n-m-1): e + (n-m)*b + (q-p+1)*f 
IHEYGS 
aaa eae aR ae Irene re 3 TEA 1 
| 30 | 40 | 50 | 65 | 75 | 
frem fm mmm ae j------ ł------ 1 
a | 1232 | 430 | 182 | 49.5 | 33.3 | 
b | 461 | 121 | 37.6 | 9.8 | 5.3 | 
c | 1871 | 623 | 259 | 69.6 | 45.4 | 
d | 490 | 128 | 40.9 | 10.5 | 5.7 | 
e | 2140 | 733 | 304 | 83.6 | 54.7 | 
£ | 29 | 7.5 | 3.3 | 0.7 | O.4 | 
ss Le===2s E rie AA Le J 
IHEYGL 
—P— quem qoem qoem quee 1 
| 30 | 40 | 50 | 65 | 75 | 
--------- ł------- ł------- ł------- ł------ 1 
a | 1320 | 451 | 189 | 49.6 | 33.3 | 
b | 1241 | 308 | 56.9 | 13.1 | 7.3 | 
c | 1959 | 644 | 266 | 69.7 | 45.4 | 
d | 1270 | 316 | 60.2 | 13.7 | 7.7 | 
e | 2228 | 755 | 311 | 83.6 | 54.7 | 
f| 29 | 7.5 | 3.3 | 0.7 | 0.4 | 
m OD Lor Loi P2232 lecca J 
IHEYGW 
Eu: BEE a cas quem qc EE qoos 1 
| 30 | 40 | 50 | 65 | 75 | 
---4------ ł------- ł------- ł------- ł------ 1 
a | 1396 | 475 | 198 | 54.9 | 36.1 | 
b | 1672 | 425 | 126 | 30.5 | 15.0 | 
c | 2035 | 667 | 275 | 75.0 | 48.1 | 
d | 1701 | 432 | 129 | 31.2 | 15.4 | 
e | 2304 | 775 | 320 | 88.9 | 57.4 | 
f | 29 | 7.5 | 3.3 | 0.7 | 0.4 | 
IM A Le ds Loa i J 
Chapter 4. Arrays: Array Functions 79 


IHEYGZ 


75 


1 65 | 


50 


40 


30 


| | ao | so | so | es | 75 | 
ee ee ee 

| 

| 

| 

| 

| 

L 
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ABS 
complex fixed-point 34-35 
complex floating-point 35-36 
ADD 
complex arguments 31 
real arguments 29 


ALL 72-73 
'and' operator 8-9 
ANY 72-73 


array indexers 
interleaved arrays 71-72 
simple arrays 70-71 
assignment operations 
bit string 12-13 
character string 16-17 
ATAN 
complex arguments 67-68 
real arguments 49-51 
ATAND (real arguments) 49-51. 
ATANH 
complex arguments 67-68 
real arguments 54-56 


bit string operators 8-10 
BOOL 14 


comparison operator 
bit string 
kyte-aligned 11 
general 11-12 
character string 15-16 
concatenation operator 
bit string 10-11 
character string 14-15 
COS 
complex arguments 62-65 
real arguments 44-46 
COSD (real arguments) 44-46 
COSH 
complex arguments 62-65 
real arguments 51-53 


DIVIDE (complex fixed-point) 32-34 
division operator 

complex fixed-point 24 

complex floating-roint 26 


ERF (real arguments) 56-58 
ERFC (real arguments) 56-58 
EXP 
complex arguments 60-61 
real arguments 40-42 
exponentiation operator 
complex operations 
floating-point exponents 28-29 
integer exponents 26-28 
real operations 
floating-point exponents 22-23 
integer exponents 20-22 


INDEX 


fill operations 
bit string 12-13 
character string 16-17 


HIGH 16-17 


IHEABU 
see: ABS (complex fixed-point) 
IHEABV 
see: ABS (complex fixed-point) 
IHEABW 
see: ABS (complex floating-point) 
IHEABZ 
see: ABS (complex floating-point) 
IHEADD 
see: ADD (real arguments) 
IHEADV 
see: ADD (complex arguments) 
IHEAPD 
see: shift-and-assign, shift-and-load 
(real operations) 
IHEATL 
see: ATAN (real arguments); ATAND (real 
arguments) 
IHEATS 
see: ATAN (real arguments); ATAND (real 
arguments) 
IHEATW 
see: ATAN (complex arguments); ATANH 
(complex arguments) 
IHEATZ 
see: ATAN (complex arguments); ATANH 
(complex arguments) 
IHEBSA 
see: ‘and' operator 
IHEBSC 
see: comparison operator (bit string, 
byte-aligned) 
IHEBSD 
see: comparison operator (bit string, 
general) 
IHEBSF 
see: BOOL 
IHEBSI 
| see: INDEX (bit string) 
IHEBSK 
see: concatenation operator (bit string); 
REPEAT (bit string) 
IHEBSM 
see: assignment operations (bit string); 
fill operations (bit string) 
IHEBSN 
see: 'not' operator 
IHEBSO 
see: 'or' operator 
IHEBSS 
see: SUBSTR (bit string) 
IHECSC 
see: comparison operator (character 
string) 
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IHECSI 
see: INDEX (character string) 
IHECSK 
see: concatenation operator (character 
string); REPEAT (character string) 
IHECSM 
see: assignrent operations (character 
string); fill operations (character 
string); HIGH; LOW 
IHECSS 
see: SUBSTR (character string) 
IHEDVU 
see: DIVIDE (complex fixed-point) 
IHEDVV 
see: DIVIDE (complex fixed-point) 
IHEDZW 
see: division operator (complex 
floating-point) 
IHEDZZ 
see: division operator (complex 
floating-point) 
IHEEFL 
see: ERF (real arguments); ERFC (real 
arguments) 
IHEEFS 
see: ERF (real arguments); ERFC (real 
arguments) 
IHEEXL 
see: EXP (real arguments) 
IHEEXS 
see: EXP (real arguments) 
IHEEXW 
See: EXP (complex arguments) 
IHEEXZ 
see: EXP (complex arguments) 
IHEHTI 
see: ATANH (real arguments) 
IHEHTS 
see: ATANH (real arguments) 
IHEJXI 


See: array indexers (interleaved arrays) 


IHEJXS 
see: array indexers (simple arrays) 
IHELNI 
see: LOG (real arguments); LOG2; LOG10 
IHELNS 
see: LOG (real arguments); LOG2; LOG10 
IHELNW 
see: LOG (complex arguments) 
IHELNZ 
see: LOG (complex arguments) 
IHEMPU 
see: MULTIPLY (complex fixed-point) 
IHEMPV 
see: MULTIPLY (complex fixed-point) 
IHEMXB 
see: MAX (real arguments); MIN (real 
arguments) 
IHEMXD 
see: MAX (real arguments); MIN (real 
arguments) 
IHEMXI 
See: MAX (real arguments); MIN (real 
arguments) 
IHENXS 
see: MAX (real arguments); MIN (real 
arguments) 
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IHEMZU 
see: multiplication operator (complex 
fixed-point); division operator 
(complex fixed-point) 
IHEMZV 
see: multiplication operator (complex 
fixed-point); division operator 
(complex fixed-point) 
IHEMZW 
see: multiplication operator (complex 
floating-point) 
IHEMZZ 
see: multiplication operator (complex 
floating-point) 
IHENL1 
See: ALL; ANY 
IHENL2 
See: ALL, ANY 
IHEPDF 
see: PROD 
IHEPDL 
see: PROD 
IHEPDS 
see: PROD 
IHEPDW 
see: PROD 
IHEPDX 
see: PROD 
IHEPDZ 
see: PROD 
IHEPSF 
see: PROD 
IHEPSL 
see: PROD 
IHEPSS 
see: PROD 
IHEPSW 
see: PROD 
IHEPSX 
see: PROD 
IHEPSZ 
see: PROD 
IHESHL 
see: SINH (real arguments); COSH (real 
arguments) 
IHESHS 
see: SINH (real arguments); COSH (real 
arguments) 
IHESMF 
see: SUM 
IHESMG 
See: SUM 
IHESMH 
See: SUM 
IHESMX 
see: SUM 
IHESNL ' 
see: SIN (real arguments); SIND (real 
arguments); COS (real arguments); COSD 
(real arguments) ; 
IHESNS 
see: SIN (real arguments); SIND (real 


arguments); COS (real arguments); COSD 


(real arguments) 
IHESNW 
see: SIN (complex arguments); SINH 
(complex arguments); COS (complex 
arguments); COSH (complex arguments) 


IHESNZ 
see: SIN (complex arguments); SINH 
(complex arguments); COS (complex 
arguments); COSH (complex arguments) 
IHESQL 
see: 
IHESQS 
see: 
IHESQW 
see: 
IHESQZ 
See: 
IHESSF 
see: 
IHESSG 
see: 
IHESSH 
See: 
IHESSX 
See: 
IHETHI 
see: 
IBETHS 
See: 
IHETNI 
See: TAN (real arguments); TAND (real 
arguments) 
IHETNS 
See: TAN (real arguments); TAND (real 
arguments) 
IHETNW 
see: TAN (complex arguments); TANH 
(complex arguments) 
IHETNZ 
see: TAN (complex arguments); TANH 
(complex arguments) 
IHEXIB 
See: exponentiation operator (real 
tions, integer exponents) 
IHEXID 
See: exponentiation operator (real 
tions, integer exponents) 
IHEXII 
see: exponentiation operator (real 
tions, integer exponents) 
IHEXIS 
see: exponentiation operator (real 
tions, integer exponents) 
IHEXIU 
See: exponentiation operator (complex 
operations, integer exponents) 
IHEXIV 
See: exponentiation operator (complex 
operations, integer exponents) 
IHEXIW 
See: exponentiation operator (complex 
opefations, integer exponents) 
iHEXIZ 
See: exponentiation operator (complex 
Operations, integer exponents) 
IHEXXL 
see: 
tions, 
IHEXXS 
see: exponentiation operator (real opera- 
tions, floating-point exponents) 


SQRT (real arguments) 


SQRT (real arguments) 
SQRT (complex arguments) 
SQRT (complex arguments) 
SUM 

SUM 

SUM 

SUM 

TANH (real arguments) 


TANH (real arguments) 


opera- 


opera- 


opera- 


opera- 


exponentiation operator (real opera- 
floating-point exponents) 


IHEXXW 
see: exponentiation operator (complex 
operations, floating-point exponents) 
IHEXXZ 
See: exponentiation operator (complex 
operations, floating-point exponents) 
IHEYGF 
See: 
IHEYGL 
see: 
IHEYGS 
See: 
IHEYGW 
see: 
IHEYGX 
See: 
IHEYGZ 
see: 
INDEX 
bit string 13-14 
character string 17-18 
indexers 
see: array indexers 


POLY 
POLY 
POLY 
POLY 
POLY 


POLY 


LOG 
complex arguments 61-62 
real arguments 42-44 


LOG2 (real arguments) 42-44 


LOG10 (real arguments) 42-44 
LOW 16-17 

MAX (real arguments) 30-31 
MIN (real arguments) 30-31 


multiplication operator 
complex fixed-point 24-25 
complex floating-point 25-26 


MULTIPLY (complex fixed-point) 31-32 


"not" operator 10 
‘or’ 


operator 9 


POLY 78-80 
PROD 75-78 


REPEAT 
bit string 10-11 
character string 14-15 


shift-and-assign, shift-and-load (real 


operations) 23-24 
SIN 
complex arguments 62-65 
real arguments 44-46 
SIND (real arguments) 44-46 
SINH 
complex arguments 62-65 
real arguments 51-53 


Index 
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SQRT 
complex arguments 58-60 
real arguments 39-40 
SUBSTR 
bit string 13 
character string 17 
SUM 73-75 
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TAN 
complex arguments 65-67 
real arguments 47-49 
TAND (real arguments) 47-49 
TANH 
complex arguments 65-67 
real arguments 53-54 


READER'S COMMENTS 


Title: IBM System/360 Operating System Form: C28-6590-0 
PL/I Subroutine Library 
Computational Subroutines 


Is the material: Yes No 


Easy to read? a= = 
Well organized? == — 


Complete? == Ss sea 
Well illustrated? = c 
Accurate? se =D 
Written for your technical level? -—- -—— 


How did you use this publication? 


----- As an introduction to the subject —----For additional 
knowledge 


Please check the items that describe your position: 


oo Customer personnel ----Operator -----Sales Representative 
-----IBM personnel ----Programmer =  -—--- Systems Engineer 
--—-—- Manager ----Customer Engineer —----Trainee 

-----Systems Analyst --—-Instructor |. | |  ----- Other--------------- 


Please check specific criticisms, give page numbers, and explain below: 


----- Clarification on pages 
e Addition on pages 
-----Deletion on pages 
----- Error on pages 


Explanation: 


If you wish a reply, be sure to include your name and address. 


C28-6590-0 


FIRST CLASS 
PERMIT NO. 33504 
NEW YORK, N.Y. 


BUSINESS REPLY MAIL 


NO POSTAGE NECESSARY IF MAILED IN THE UNITED STATES 





POSTAGE WILL BE PAID BY... 


IBM CORPORATION 


1271 AVENUE OF THE AMERICAS 
NEW YORK, N.Y. 10020 


ATTENTION: PUBLICATIONS, DEPT. D39 


® 0 a 0 0000 66 è 0 è 0 0 0 0 è è è e e 0 è è e e è 0 0 e è e e è è e e a e è e è a e e e e 0 e rn... 0 è è e 0 0 e 0 0 è e 0 S 0 è è € è è è S è 0 0 è e 0 0 0 e e 0 0 è è è è è è e 6 è sè è è è se 0 è 0 0 è è 0 è è 0 0 0 è 0 HOD * 


BIMI 


International Business Machines Corporation 
Data Processing Division 
112 East Post Road, White Plains, N.Y. 10601 


t paund | 
germ 


vsqui | 


0-0659-879 - 


C28-6590-0 — 


SIM 


® 
International Business Machines Corporation 
Data Processing Division 
112 East Post Road, White Plains, N.Y. 10601 


VSN ul paund 


0-06S9-8Z9 


